线性推,矩阵乘法+快速幂求通项。
传送门:点击打开链接
#include <cstdio>
#include <cstring>
#include <cmath>
#include <iostream>
using namespace std;
#define LL long long
struct Mat{
LL f[2][2];
};
LL MOD;
Mat mul(Mat a,Mat b)
{
LL i,j,k;
Mat c;
memset(c.f,0,sizeof(c.f));
for(i=0;i<2;i++)
for(j=0;j<2;j++)
for(k=0;k<2;k++)
c.f[i][j]=(c.f[i][j]+a.f[i][k]*b.f[k][j])%MOD;//可以改为不%MOD
return c;
}
Mat pow_mod(Mat e,LL b)
{
Mat s;
s.f[0][0]=s.f[1][1]=1;
s.f[0][1]=s.f[1][0]=0;
while(b)
{
if(b&1)
s=mul(s,e);
e=mul(e,e);
b=b>>1;
}
return s;
}
int main()
{
LL a,b,n,m;
//while(~scanf("%lld%lld%lld%lld",&a,&b,&n,&m))
while(~scanf("%I64d%I64d%I64d%I64d",&a,&b,&n,&m))
{
MOD=m;
LL p,q,ans;
Mat e;
p=2*a;
q=a*a-b;
e.f[0][0]=p;e.f[0][1]=1;
e.f[1][0]=-q;e.f[1][1]=0;
e=pow_mod(e,n-1);
ans=((p*e.f[0][0]+2*e.f[1][0])%MOD+MOD)%MOD; //可能负数结果
//printf("%lld\n",ans);
printf("%I64d\n",ans);
}
return 0;
}