题目链接:HDU 4291
题意:g(0)=0,g(1)=1,g(n)=3g(n-1)+g(n-2),mod=1000000007。求g(g(g(n)))%mod。
分析:
看了题解才发现是循环节....... 完全不会啊.
这是别人博客的解释:https://blog.csdn.net/weixin_34233421/article/details/86059283
对于ans=g(g(n))%mod,如果g(g(n)%mod)%mod这样显然是不对的。g(g(n)%L1)%mod才是正确的。为什么呢?因为对于本题的g(0)=0,g(1)=1来说,对于mod每L1次就会重新回来,所以你在计算的过程中用L1去模掉中间值跟最后模掉mod是一样的!!!
这是循环节代码:
void get_mod(int mod)
{
long long a,b,c;
a=0;b=1;
int i=0;
while(1)
{
i++;
c=3*b+a;
c%=mod;
a=b;
b=c;
if(a==0&&b==1)
{
printf("%d\n",i);
return ;
}
}
}
上面代码跑出来的3个循环节
const int mod1=1e9+7;
const int mod2=222222224;
const int mod3=183120;
剩下就简单了,套个快速幂板子完事.
以下是代码:
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=3e5+10;
const int mod1=1e9+7;
const int mod2=222222224;
const int mod3=183120;
const int maxn=4;
int Mod;
ll tmp[maxn][maxn];
void multi(ll a[][maxn],ll b[][maxn],int n)//n是矩阵大小
{
int i,j,k;
memset(tmp,0,sizeof tmp);
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
for(k=1;k<=n;k++)
{
tmp[i][j]+=a[i][k]*b[k][j];
tmp[i][j]%=Mod;
}
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
a[i][j]=tmp[i][j];
}
ll res[maxn][maxn];
void Pow(ll a[][maxn],ll m,int n,int md)// a 是初始矩阵 res是答案数组 m是幂,n是矩阵大小
{
if(m==-1)
{
a[1][1]=0;
return ;
}
Mod=md;
for(int i=1;i<maxn;i++)
for(int j=1;j<maxn;j++)
res[i][j]=(i==j);
while(m)
{
if(m&1) multi(res,a,n); //res=res*a;
multi(a,a,n); //a=a*a
m>>=1;
}
}
void Show(int n,ll res[][maxn])
{
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
cout<<res[i][j]<<" ";
}
cout<<endl;
}
}
ll a[maxn][maxn];
int main()
{
ll n;
while(~scanf("%lld",&n))
{
a[1][1]=3;a[1][2]=1;
a[2][1]=1;a[2][2]=0;
Pow(a,n-1,2,mod3);
a[1][1]=3;a[1][2]=1;
a[2][1]=1;a[2][2]=0;
Pow(a,res[1][1]-1,2,mod2);
//Show(2,res);
a[1][1]=3;a[1][2]=1;
a[2][1]=1;a[2][2]=0;
Pow(a,res[1][1]-1,2,mod1);
printf("%lld\n",res[1][1]);
}
return 0;
}
/*
0 0
1 1
3 2
10 3
33 4
*/