数学-HDU 4565 so easy

一道数学味很浓的题，如果不是看了题解真忘了那些数学知识了····

这是解析

http://blog.csdn.net/crazy______/article/details/9021169

看明白了直接矩阵快速幂就可以了（原谅我英语不好···）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include <algorithm>
#define ll long long
#define eps 1e-8
#define ms(x,y) (memset(x,y,sizeof(x)))
#define fr(i,x,y) for(ll i=x;i<=y;i++)
using namespace std;
struct juzhen
{
ll p[3][3];
juzhen(){};
};
juzhen mul(juzhen a,juzhen b)
{
juzhen c;
c.p[1][1]=a.p[1][1]*b.p[1][1]+a.p[1][2]*b.p[2][1];
c.p[1][2]=a.p[1][1]*b.p[1][2]+a.p[1][2]*b.p[2][2];
c.p[2][1]=a.p[2][1]*b.p[1][1]+a.p[2][2]*b.p[2][1];
c.p[2][2]=a.p[2][1]*b.p[1][2]+a.p[2][2]*b.p[2][2];
return c;
}
juzhen fast_mod(juzhen a,ll n,ll m)
{
juzhen ret;
ret.p[1][1]=ret.p[2][2]=1;
ret.p[1][2]=ret.p[2][1]=0;
while(n)
{
if(n&1)
{
ret=mul(ret,a);
fr(i,1,2)fr(j,1,2)
{
ret.p[i][j]%=m;
if(ret.p[i][j]<0)ret.p[i][j]+=m;
}
}
a=mul(a,a);
fr(i,1,2)fr(j,1,2)
{
a.p[i][j]%=m;
if(a.p[i][j]<0)a.p[i][j]+=m;
}
n/=2;
}
return ret;
}
int main()
{
ll a,b,n,m;
while(~scanf("%I64d%I64d%I64d%I64d",&a,&b,&n,&m))
{
juzhen ans;
ans.p[1][1]=2*a;
ans.p[1][2]=-(a*a-b);
ans.p[2][1]=1;
ans.p[2][2]=0;
ans=fast_mod(ans,n,m);
ll t=ans.p[2][1]*2*a+2*ans.p[2][2];
cout<<t%m<<endl;
}
}