知天易逆天难(未解决)

知天易，逆天难！（未解决）
#include<iostream>
#include<math.h>
using namespace std;
char s[100011];
int phi(int n)
{
int i,t=(int)sqrt((double)n),p=n; for(i=2;i<=t;i++)
{
if(p%i==0)
{
while(p%i==0)
p=p/i;
n=n-n/i;
}
}
if(p>1)
return n-n/p;
return n;
}
void sub(int p)
{
int a,i,len=strlen(s),t=0; for(i=len-1;i>=0;i--)
{
if(t==0 && p==0)
break;
a=p%10;
s[i]=s[i]-'0'+t-a;
if(s[i]<0)
{
s[i]+=10;
t=-1;
}
else
t=0;
s[i]=char(s[i]+48);
p=p/10;
}
}
int mod(int p,int x)
{
int i,len=strlen(s),sum=0; if(x>0)
{
int t=0,temp;
temp=s[len-1]-'0'-1;
if(temp<0)
{
temp=temp+10;
t=-1;
}
s[len-1]=char(temp+48); for(i=len-2;i>=0;i--)
{
temp=s[i]-'0'+t;
if(temp<0)
{
temp=temp+10;
t=-1;
}
else
t=0;
s[i]=char(temp+48);
}
}
//cout<<s<<endl;
for(i=0;i<len;i++)
{
sum=sum*10+s[i]-'0';
if(sum>p)
sum=sum%p;
}
return sum;
}
__int64 epow(int n,int p) {
__int64 s=1,t=2;
while(n)
{
if(n%2==1)
s=(s*t)%p;
t=(t*t)%p;
n=n>>1;
}
return s;
}
int main()
{
int i,c,m,r,t,p,len,sum; __int64 ans,temp;
while(scanf("%s",s)!=EOF) {
m=0;
p=1;
sum=0;
scanf("%d",&c);
len=strlen(s);
if(c==0)
{
printf("0\n");
continue;
}
if(len<=8)
{
int q=0;
for(i=0;i<len;i++)
sum=sum*10+s[i]-'0';
//cout<<sum<<endl;
ans=epow(sum-1,c);
//printf("%I64d\n",ans); ans=ans*sum%c;
printf("%I64d\n",ans); continue;
}
while(c%2==0)
{
m++;
p=p*2;
c=c/2;
}
//cout<<m<<endl;
r=phi(c);
t=mod(c,0);
sub(m);
temp=mod(r,1);
//cout<<s<<endl;
//printf("%I64d\n",temp); temp=epow(temp,c);
ans=t*temp%c;
ans=ans*p;
printf("%I64d\n",ans); }
return 0;
}。