題意:給定a,b,n,p,計算G(a,b,n,p)的值。
其中,
且1 ≤a, b, n, p≤2*109 ,p爲奇素數,a,b<p.
解題思路:G(a,b,n,p)分爲三個部分相乘,於是可將三部分分開來分別取餘數,最後再計算G(a,b,n,p)的值。
對於前面兩個括號部分,可以分別應用二分快速冪取模算法得到結果。
對於後面的中括號的內容,可以通過推導證明括號裏面的兩式之和一定爲整數,能夠取模。
推導過程如下:
因此,後面中括號取模的值即爲2Xn對P取模的值。
根據公式,a^b%c = a^(b%phi(c)+pphi(c))%c (如果b>phi(c),phi(c)爲c的歐拉函數值)
可將括號中兩式的指數F(n)模p-1後變爲power,這個過程用矩陣快速冪實現。
然後可推出Xn與Xn-1的遞推關係,Yn與Yn-1的遞推關係,用矩陣快速冪即可求解。
代碼如下:
#include <stdio.h>
typedef struct
{
__int64 matrix[2][2];
}Matrix;
__int64 powermod(__int64 x,__int64 n,__int64 p)
{
__int64 res = 1;
for(;n;n>>=1)
{
if(n&1)
res = (res*x)%p;
x = (x*x)%p;
}
return res;
}
Matrix multi(Matrix x,Matrix y,__int64 p)
{
Matrix res;
int i,j,k;
__int64 sum;
for(i = 0;i<2;i++)
for(j = 0;j<2;j++)
{
sum = 0;
for(k = 0;k<2;k++)
sum+=(x.matrix[i][k]*y.matrix[k][j])%p;
res.matrix[i][j] = sum%p;
}
return res;
}
Matrix Mat_powermod(Matrix x,__int64 n,__int64 p)
{
Matrix res;
int i,j;
for(i = 0;i<2;i++)
for(j = 0;j<2;j++)
{
if(i==j)
res.matrix[i][j] = 1;
else
res.matrix[i][j] = 0;
}
for(;n;n>>=1)
{
if(n&1)
res = multi(res,x,p);
x = multi(x,x,p);
}
return res;
}
int main()
{
__int64 p,T,a,b,n;
__int64 keep1,keep2,keep3;
__int64 power,result;
Matrix res;
scanf("%I64d",&T);
while(T--)
{
scanf("%I64d%I64d%I64d%I64d",&a,&b,&n,&p);
keep1 = (powermod(a,(p-1)/2,p)+1)%p;
keep2 = (powermod(b,(p-1)/2,p)+1)%p;
if(n==0)
power = 1;
else
{
res.matrix[0][0] = 1;
res.matrix[0][1] = 1;
res.matrix[1][0] = 1;
res.matrix[1][1] = 0;
res = Mat_powermod(res,n-1,p-1);
power = (res.matrix[0][0]+res.matrix[0][1])%(p-1);
}
power+=p-1;
res.matrix[0][0] = (a+b)%p;
res.matrix[0][1] = (2*a*b)%p;
res.matrix[1][0] = 2%p;
res.matrix[1][1] = (a+b)%p;
res = Mat_powermod(res,power-1,p);
keep3 = (2*(((res.matrix[0][0]*(a+b))%p+(res.matrix[0][1]*2)%p)%p))%p;
result = 1;
result = (result*keep1)%p;
result = (result*keep2)%p;
result = (result*keep3)%p;
printf("%I64d\n",result);
}
return 0;
}