UVALive-3722 留个坑，为什么费马小定理求逆元不对？？

#include <iostream>
#include <cstdlib>
#include <queue>
#include <algorithm>
#include <cmath>
#include <fstream>
#define inf 2147483647
#define N 1000010
#define p(a) putchar(a)
#define For(i,a,b) for(long long i=a;i<=b;++i)

using namespace std;
long long x,a,n,c,t,k;
void in(long long &x){
    long long y=;char c=getchar();x=;
    while(c<''||c>''){if(c=='-')y=-;c=getchar();}
    while(c<=''&&c>=''){ x=(x<<)+(x<<)+c-'';c=getchar();}
    x*=y;
}
void o(long long x){
    if(x<){p('-');x=-x;}
    if(x>)o(x/);
    p(x%+'');
}

long long ksm(long long a,long long b){
    long long r=;
    while(b>){
        if(b&)
            r=r*a%c;
        a=a*a%c;
        b>>=;
    }
    return r;
}

long long gcd(long long a,long long b){
    return (a%b==?b:gcd(b,a%b));
}

long long sum(long long a,long long b){
    if(b==)
        return a%c;
    long long half=ksm(a,b>>);
    long long temp=sum(a,b>>);
    if(b%==)
        return (temp%c+half*temp%c)%c;
    else
        return (temp%c+half*temp%c+half*half%c*a%c)%c;
}

int main(){
    while(cin>>x>>a>>n>>c&&x){
        t=ksm(a,n);
        //if(gcd(a-1,c)==1)
        //    o(((x%c*t%c-a%c*ksm(a-1,c-2)%c*((t-1)%c+c)%c)%c+c)%c);
        //else
            o(((x%c*t%c-sum(a,n))%c+c)%c);
        p('\n');
    }
    return ;
}

https://vjudge.net/problem/UVALive-3722