http://acm.hdu.edu.cn/showproblem.php?pid=6395
p/n 取值只有 sqrt(p) 级数 ,然后想同的p/n分成一块,因为连续递减,二分每个块的右边界(好吧,是我懒得想怎么快速求有边界),然后每一块矩阵快速幂。
至于为什么是sqrt的级数,,,,emmm,,,,不会证明,,,,但是感觉分的块数肯定爆减。
转移方程:
#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define FI first
#define SE second
const LL mod =1e9+7;
const int MX = 1e6+5;
struct mat{
LL m[3][3];
};
mat matMul(mat aa,mat bb,int sz){
mat re;memset(re.m,0,sizeof(re.m));
for(int i=0;i<sz;i++)
for(int j=0;j<sz;j++)
for(int k=0;k<sz;k++)
re.m[i][j]=(re.m[i][j]+aa.m[i][k]*bb.m[k][j])%mod;
return re;
}
mat matQuickMod(mat jin,LL x,int sz){
mat re;memset(re.m,0,sizeof(re.m));
for(int i=0;i<sz;i++) re.m[i][i]=1;
while(x){
if(x&1) re=matMul(re,jin,sz);
x=x>>1;jin=matMul(jin,jin,sz);
}
return re;
}
LL getlager(LL id,LL P,LL n){
LL l=id,r=n,re=-1;
while(l<=r){
LL mid=(l+r)/2;
if(P/id==P/mid)re=mid,l=mid+1;
else r=mid-1;
}
return re;
}
int main(){
int T;cin>>T;while(T--){
LL A,B,C,D,P,n;
cin>>A>>B>>C>>D>>P>>n;
mat y;
y.m[0][0]=D;y.m[0][1]=1;y.m[0][2]=0;
y.m[1][0]=C;y.m[1][1]=0;y.m[1][2]=0;
y.m[2][0]=0;y.m[2][1]=0;y.m[2][2]=1;
LL L=3,R;
mat x;memset(x.m,0,sizeof(x.m));
for(int i=0;i<3;i++)x.m[i][i]=1;
while(L<=n){
R=getlager(L,P,n);
y.m[2][0]=P/L;
x=matMul(x,matQuickMod(y,R-L+1,3),3);
L=R+1;
}
mat ans;memset(ans.m,0,sizeof(ans.m));
ans.m[0][0]=B;ans.m[0][1]=A;ans.m[0][2]=1;
ans=matMul(ans,x,3);
cout<<ans.m[0][0]<<endl;
}
return 0;
}