题目: Ipad,IPhone
思路:这题。。。好坑反正,WA了刚好20次。。。
题目给的式子,分成两部分,对于第一部分,直接快速幂即可,对于后面的部分,后面括号里是整数的证明 出题人 AekdyCoin 用二次剩余进行了解说,请点击 here
对于后面的部分,是不是整数我给直接忽略了。。。不然怎么进行取余呢。。事实证明,前半部分是有存在感的,根据欧拉准则,即二次剩余的欧拉判别法 ,我们知道a是模p的二次剩余时a^((p-1)/2)=1(mod p),非二次剩余时,a^((p-1)/2)=-1(mod p) ,所以二次剩余时,前半部分的解为4,而非二次剩余的时候,前面的解为0,同时后半部分也不会有整数解,因此前半部分的存在时很有必要的。
好吧,继续说后半部分,既然是整数可以进行取余了,指数部分是菲波那契数列数列,可以构造矩阵这个不难,同时要用指数循环节进行优化,不然很容易溢出,然后就跪了。这里对素数p进行取余,所以指数部分,我们对其欧拉值进行取余,也就是对p-1取余得到斐波那契取余项。
然后再构造矩阵:
然后。。。因为矩阵初始化的时候忘了对2*(a+b),-(a+b)^2等进行取余,WA了差不多20次吧。。。
#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
using namespace std;
typedef __int64 LL;
LL mod;
struct Matrix
{
LL m[3][3];
}E,D;
LL a,b,n;
void init()
{
for(int i=1;i<=2;i++)
for(int j=1;j<=2;j++)
E.m[i][j]=(i==j);
D.m[1][1]=1;
D.m[1][2]=1;
D.m[2][1]=1;
D.m[2][2]=0;
}
Matrix Multi(Matrix A,Matrix B)
{
Matrix ans;
for(int i=1;i<=2;i++)
for(int j=1;j<=2;j++)
{
ans.m[i][j]=0;
for(int k=1;k<=2;k++)
ans.m[i][j]=(ans.m[i][j]+A.m[i][k]*B.m[k][j])%mod;
}
return ans;
}
Matrix Pow(Matrix A,LL k)
{
Matrix ans=E;
while(k)
{
if(k&1)
{
k--;
ans=Multi(ans,A);
}
else
{
k/=2;
A=Multi(A,A);
}
}
return ans;
}
LL Fib(LL n)
{
if(n==0)
return 1;
Matrix ans=Pow(D,n);
return ans.m[1][1]%mod;
}
void Print(Matrix A)
{
for(int i=1;i<=2;i++)
{
for(int j=1;j<=2;j++)
cout<<A.m[i][j]<<" ";
cout<<endl;
}
}
LL get(Matrix A,Matrix B,LL fib)
{
Matrix ans=Pow(A,fib-1);
ans=Multi(ans,B);
return ans.m[1][1];
}
LL Pow(LL a,LL b)
{
LL ans=1;
while(b)
{
if(b&1)
{
b--;
ans=(ans*a)%mod;
}
else
{
b/=2;
a=(a*a)%mod;
}
}
return ans;
}
int main()
{
init();
int t;
cin>>t;
init();
while(t--)
{
cin>>a>>b>>n>>mod;
LL aa=(1+Pow(a,(mod-1)/2))%mod;
LL bb=(1+Pow(b,(mod-1)/2))%mod;
if(aa==0||bb==0)
{
cout<<0<<endl;
continue;
}
Matrix T,cnt;
T.m[1][1]=2*(a+b)%mod;
T.m[1][2]=-(a-b)*(a-b)%mod;
T.m[2][1]=1;
T.m[2][2]=0;
cnt.m[1][1]=2*(a+b)%mod;
cnt.m[2][1]=2;
LL tmp;
mod--;
tmp=Fib(n);
mod++;
LL ans=get(T,cnt,tmp);
ans=4*ans%mod;
if(ans<0)
ans+=mod;
cout<<ans<<endl;
}
return 0;
}