Длинка

#include <cstdio>
#include <iostream>
#include <algorithm>
#define ll long long
#define N 11111
#define P 1000000000
using namespace std;

class bigint{
    public:
    int sz;
    ll a[N];
    bigint(ll x=0){
        memset(a,0,sizeof(a));
        sz=0;
        while (x){
            a[++sz]=x%P;
            x/=P;
        }
        if (!sz) ++sz;
    }
    void print(){
        printf("%I64d",a[sz]);
        for (int i=sz-1;i>=1;i--) printf("%.9I64d",a[i]);
    }
    friend ostream& operator << (ostream& output, bigint & q){
        q.print();
        return output;
    }
    ll & operator [] (int x){
        return a[x];
    }
    friend bool operator < (bigint & p,bigint &q){
        if (p.sz<q.sz) return 1;
        if (p.sz>q.sz) return 0;
        for (int i=p.sz;i;i--){
            if (p[i]>q[i]) return 0;
            if (p[i]<q[i]) return 1;
        }
        return 0;
    }
    friend bool operator == (bigint & p,bigint & q){
            return !(p<q) && !(q<p);
    }
    friend bool operator != (bigint & p, bigint & q){
        return p<q || q<p;
    }
    friend bool operator >= (bigint & p,bigint & q){
        return !(p<q);
    }
    friend bool operator <= (bigint &p, bigint & q){
        return !(q<p);
    }
    friend bool operator > (bigint &p, bigint & q){
        return q>p;
    }
    friend bigint operator + (bigint &a,bigint &b){
        bigint res;
        int n=max(a.sz,b.sz);
        for (int i=1;i<=n;i++){
            res[i]=a[i]+b[i]+res[i];
            res[i+1]=res[i]/P;
            res[i]%=P;
        }
        if (res[n+1]) ++n;
        res.sz=n;
        return res;
    }
    friend bigint operator * (bigint &a,bigint &b){
        bigint res;
        int n=a.sz,m=b.sz,k=n+m-1;
        for (int i=1;i<=n;i++)
            for (int j=1;j<=m;j++){
                res[i+j-1]=a[i]*b[j]+res[i+j-1];
                res[i+j]+=res[i+j-1]/P;
                res[i+j-1]%=P;
            }
        if (res[k+1]) ++k;
        res.sz=k;
        return res;
    }
    friend bigint operator - (bigint &a,bigint &b){
        bigint res;
        int k=a.sz;
        for (int i=1;i<=k;i++){
            res[i]+=a[i]-b[i];
            if (res[i]<0){
                res[i]+=P;
                res[i+1]--;
            }
        }
        while (res[k]==0) --k;
        res.sz=k;
        return res;
    }
    friend void operator -=(bigint &p,bigint & q){
        p=p-q;
    }
    friend void operator +=(bigint &p,bigint & q){
        p=p+q;
    }
    friend void operator *=(bigint &p, bigint & q){
        p=p*q;
    }
    friend bigint operator ^ (bigint & p, int y){
        bigint x=p,res=1;
        while (y){
            if (y&1) res*=x;
            y>>=1;
            x*=x;
        }
        return res;
    }
    friend bigint operator ! (bigint & p){
        bigint res=1,I;
        int n=p[1];
        for (int i=1;i<=n;i++){
            I=i;
            res=res*I;
        }
        return res;
    }

};

int main(){
    int x,y;
    cin>>x>>y;
    bigint f1=x;
    bigint f=f1^y;
    cout<<f1<<"^"<<y<<"="<<f;
}