HDU6395 Sequence 矩阵快速幂+整除分块

题意：

         给你式子，让你求F[n]是多少，一看就知道要用到矩阵快速幂，但是这道题恶心就在于p/n这个东西，是动态变化的，不是一个常数，但是我们可以想办法令他变成常数，怎么搞呢？用到一个叫整除分块的小玩意来搞这个，可以达到根号n复杂度。然后相同p/n的区间就作为一个常数去搞就可以了，还有一点应该注意的就是p/(p/l)有可能比n还要大，导致出现错误，所以要写成r=min(n,p/(p/l)).就这么多了

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define LL long long int
const LL MOD=1e9+7;
struct node
{
	LL m[3][3];
	node()
	{
		memset(m,0,sizeof(m));
	}
};
LL A,B,C,D,P,n,E,ans;
node cla(node a,node b)//这里可以 node cla(node a,node b,int l) 第三个l表示的是你矩阵的大小，但是这里求斐波那契，只需要2就好了，就不用写了
{
/*	for(int i=0;i<3;i++)
	{
		for(int j=0;j<3;j++)
		printf("%lld ",b.m[i][j]);
		printf("\n");
	}
    printf("a%lld b%lld  a%lld b%lld %lld %lld\n",a.m[1][0],b.m[0][0],a.m[1][1],b.m[1][0],a.m[1][2],b.m[2][0]);*/
	node c;
	for(int i=0;i<3;i++)
		for(int j=0;j<3;j++)
			for(int k=0;k<3;k++)
			if(a.m[i][k]&&b.m[k][j])//剪枝（添条件，设门槛），提高效率，有一个是0，相乘肯定是0
			{
				c.m[i][j]+=a.m[i][k]*b.m[k][j];
				c.m[i][j]%=MOD;
			}
	return c;
}
LL POW(int k)
{
	node a,c;
	for(int i=0;i<3;i++) c.m[i][i]=1;
	a.m[0][0]=D,a.m[0][1]=C,a.m[0][2]=1;
	a.m[1][0]=1,a.m[1][1]=0,a.m[1][2]=0;
	a.m[2][0]=0,a.m[2][1]=0,a.m[2][2]=1;
	while(k)
	{
		if(k%2) c=cla(c,a);
		a=cla(a,a);
		k/=2;
	}
	LL temp1=(c.m[1][0]*B%MOD+c.m[1][1]*A%MOD+c.m[1][2]*E%MOD)%MOD;
//	printf("@@%lld %lld %lld\n",c.m[1][0]*B,c.m[1][1]*A,c.m[1][2]*E);
	LL temp2=(c.m[0][0]*B%MOD+c.m[0][1]*A%MOD+c.m[0][2]*E%MOD)%MOD;
	A=temp1,B=temp2;
//	printf("%lld %lld\n",A,B);
	ans=B;
}
int main()
{
/*	freopen("in.txt","r",stdin); //输入重定向，输入数据将从in.txt文件中读取
	freopen("out.txt","w",stdout); //输出重定向，输出数据将保存在out.txt文件中*/
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld%lld%lld%lld%lld%lld",&A,&B,&C,&D,&P,&n);
        if(n<=2)
        {
            if(n==1)
                printf("%lld\n",A);
            if(n==2)
                printf("%lld\n",B);
            continue;
        }
        int r;
        for(int l=3;l<=n;l=r+1)
        {
            if(P/l==0){
            	E=0;
                POW(n-l+1);
                break;
            }
            else{
                r=min(n,P/(P/l));
                E=P/l;
            //    printf("%d %d %d %d\n",l,r,r-l+1,E);
                POW(r-l+1);
            }
        }
        printf("%lld\n",ans);
    }
}