hdu 4291 A Short problem

数学题，找循环节！！

首先g(g(g(n)))=g(x) mod 1e9+7 则可知x有循环节1e9+7;

之后x=g(g(n)),则可算出g(n)的循环节，在算出n的循环节就可以了！！

代码如下：

 #include<iostream>
 #include<stdio.h>
 #include<algorithm>
 #include<iomanip>
 #include<cmath>
 #include<cstring>
 #include<vector>
 #define ll __int64
 #define mod1 1000000007
 #define mod2 222222224
 #define mod3 183120
 #define mod4 240
 using namespace std;
 struct ma
 {
     ll a[][];
 };
 ma mul(ma aa,ma ab,ll mod)
 {
     ma ans;
     for(int i=;i<;i++)
     for(int j=;j<;j++){
         ans.a[i][j]=;
         for(int k=;k<;k++)
             ans.a[i][j]=(ans.a[i][j]+aa.a[i][k]*ab.a[k][j])%mod;
     }
     return ans;
 }
 ma pows(ma aa,ll b,ll mod)
 {
     ma ans;
     ans.a[][]=ans.a[][]=;
     ans.a[][]=ans.a[][]=;
     while(b){
         if(b&) ans=mul(ans,aa,mod);
         b>>=;
         aa=mul(aa,aa,mod);
     }
     return ans;
 }
 int main(){
     ll n;
     ma p;
     p.a[][]=;p.a[][]=p.a[][]=;p.a[][]=;
     while(cin>>n){
         ma ans;
         if(n==){
             printf("0\n");
             continue;
         }
         ans=pows(p,n%mod4,mod3);
         ans=pows(p,ans.a[][],mod2);
         ans=pows(p,ans.a[][],mod1);
         printf("%I64d\n",ans.a[][]);
     }
     return ;
 }