hdu 6395 Sequence (17多校) //矩阵快速幂+二分+分块

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;
}