fuzhou 1683 纪念SlingShot ***

Problem 1683 纪念SlingShot

Accept: 361    Submit: 1287
Time Limit: 1000 mSec    Memory Limit : 32768 KB

Problem Description

已知 F（n）=3 * F（n-1）+2 * F（n-2）+7 * F（n-3），n>=3，其中F（0）=1，F（1）=3，F（2）=5，对于给定的每个n，输出F（0）+ F（1）+ …… + F（n） mod 2009。

Input

第一行是一整数m，代表总共有m个cases。

Output

对于每个case，输出一行。格式见样例，冒号后有一个空格。

Sample Input

2 3 6

Sample Output

Case 1: 37 Case 2: 313

Source

FOJ月赛-2009年2月- Coral

 

Submit  Back  Status

 

题意不看，刚开始觉得好水。直接快速幂就可以做。

后来发现求和的，囧...

————————————————————————————

 

  1 /*
  2 矩阵构造：
  3 假设前N项和为S(n), 那么可以推出 S(n)=S(n-1)+f(n);
  4 得到一个矩阵:
  5 
  6 | 1 1 0 0 | | s(n-1) |   |  S(n)  |
  7 | 0 3 2 7 | |  f(n)  | = | f(n+1) |
  8 | 0 1 0 0 | | f(n-1) |   |  f(n)  |
  9 | 0 0 1 0 | | f(n-2) |   | f(n-1) |
 10 
 11 然后每一个N，一次快速幂，过了？额，下面就是这样的代码
 12 超时ing...
 13 
 14 那么就要优化了。优化：题目中矩阵是不变的，只是单纯的问
 15 n的值。保存2^i的结果。每一个数字都能由2^i组成。（二进制优化）
 16 
 17 这样就节省了很多的时间。
 18 484 ms
 19 
 20 特别注意：
 21 在最初的编写过程中，我只是用 (M_tom.mat[1][1]*4+M_tom.mat[1][2]*5)%mod;
 22 这样做的结果是答案偏小了。问题在于还要加上 (M_tom.mat[1][3]*3+M_tom.mat[1][4]*1);
 23 这就是数学的问题了.....囧
 24 动手试一试
 25 */
 26 
 27 /*
 28 #include<iostream>
 29 #include<cstdio>
 30 #include<cstring>
 31 #include<cstdlib>
 32 using namespace std;
 33 typedef __int64 LL;
 34 const __int64 mod=2009;
 35 
 36 struct Matrix
 37 {
 38     LL mat[6][6];
 39     void ini()
 40     {
 41         memset(mat,0,sizeof(mat));
 42     }
 43     void first_ini()
 44     {
 45         for(LL i=1;i<=4;i++)
 46         for(LL j=1;j<=4;j++)
 47         if(i==j) mat[i][j]=1;
 48         else mat[i][j]=0;
 49     }
 50 }M_hxl,M_tom;
 51 
 52 
 53 void make_first()
 54 {
 55     M_hxl.mat[1][1]=1;M_hxl.mat[1][2]=1;M_hxl.mat[1][3]=0;M_hxl.mat[1][4]=0;
 56     M_hxl.mat[2][1]=0;M_hxl.mat[2][2]=3;M_hxl.mat[2][3]=2;M_hxl.mat[2][4]=7;
 57     M_hxl.mat[3][1]=0;M_hxl.mat[3][2]=1;M_hxl.mat[3][3]=0;M_hxl.mat[3][4]=0;
 58     M_hxl.mat[4][1]=0;M_hxl.mat[4][2]=0;M_hxl.mat[4][3]=1;M_hxl.mat[4][4]=0;
 59 }
 60 
 61 Matrix Multiply(Matrix cur,Matrix now)
 62 {
 63     Matrix ww;
 64     ww.ini();
 65     for(LL i=1;i<=4;i++)
 66     for(LL k=1;k<=4;k++)
 67     if(cur.mat[i][k])
 68     {
 69         for(LL j=1;j<=4;j++)
 70         if(now.mat[k][j])
 71         {
 72             ww.mat[i][j]+=cur.mat[i][k]*now.mat[k][j];
 73             if(ww.mat[i][j]>=mod)
 74             ww.mat[i][j]%=mod;
 75         }
 76     }
 77     return ww;
 78 }
 79 
 80 void power_sum2(LL n)
 81 {
 82       M_tom.first_ini();
 83       while(n)
 84       {
 85           if(n&1)
 86           {
 87               M_tom=Multiply(M_hxl,M_tom);
 88           }
 89           n=n>>1;
 90           M_hxl=Multiply(M_hxl,M_hxl);
 91       }
 92       LL sum=0;
 93       sum=(sum+M_tom.mat[1][1]*4)%mod;
 94       sum=(sum+M_tom.mat[1][2]*5)%mod;
 95       sum=(sum+M_tom.mat[1][3]*3)%mod;//!!!开始时没有加
 96       sum=(sum+M_tom.mat[1][4]*1)%mod;//!!!开始时没有加
 97       printf("%I64d\n",sum);
 98 }
 99 
100 int main()
101 {
102     LL T,n,i;
103     while(scanf("%I64d",&T)>0)
104     {
105         for(i=1;i<=T;i++)
106         {
107             scanf("%I64d",&n);
108             printf("Case %I64d: ",i);
109             if(n==0)printf("1\n");
110             else if(n==1) printf("4\n");
111             else if(n==2) printf("9\n");
112             else
113             {
114                 make_first();
115                 power_sum2(n-1);
116             }
117         }
118     }
119     return 0;
120 }
121 
122 
123 */
124 
125 #include<iostream>
126 #include<cstdio>
127 #include<cstring>
128 #include<cstdlib>
129 using namespace std;
130 typedef __int64 LL;
131 const __int64 mod=2009;
132 
133 struct Matrix
134 {
135     LL mat[6][6];
136     void ini()
137     {
138         memset(mat,0,sizeof(mat));
139     }
140     void first_ini()
141     {
142         for(LL i=1;i<=4;i++)
143         for(LL j=1;j<=4;j++)
144         if(i==j) mat[i][j]=1;
145         else mat[i][j]=0;
146     }
147     void init()
148     {
149         mat[1][1]=1;mat[1][2]=1;mat[1][3]=0;mat[1][4]=0;
150         mat[2][1]=0;mat[2][2]=3;mat[2][3]=2;mat[2][4]=7;
151         mat[3][1]=0;mat[3][2]=1;mat[3][3]=0;mat[3][4]=0;
152         mat[4][1]=0;mat[4][2]=0;mat[4][3]=1;mat[4][4]=0;
153     }
154 }M_hxl[40],M_tom;
155 
156 Matrix Multiply(Matrix cur,Matrix now)
157 {
158     Matrix ww;
159     ww.ini();
160     for(LL i=1;i<=4;i++)
161     for(LL k=1;k<=4;k++)
162     if(cur.mat[i][k])
163     {
164         for(LL j=1;j<=4;j++)
165         if(now.mat[k][j])
166         {
167             ww.mat[i][j]+=cur.mat[i][k]*now.mat[k][j];
168             if(ww.mat[i][j]>=mod)
169             ww.mat[i][j]%=mod;
170         }
171     }
172     return ww;
173 }
174 
175 void power_sum2(LL n)
176 {
177       M_tom.first_ini();//！！
178       LL cnt=0;
179       while(n)
180       {
181           if(n&1)
182           {
183               M_tom=Multiply(M_hxl[cnt],M_tom);
184           }
185           n=n>>1;
186           cnt++;//模拟二进制.
187       }
188       LL sum=0;
189       sum=(sum+M_tom.mat[1][1]*4)%mod;
190       sum=(sum+M_tom.mat[1][2]*5)%mod;
191       sum=(sum+M_tom.mat[1][3]*3)%mod;//!!!开始时没有加
192       sum=(sum+M_tom.mat[1][4]*1)%mod;//!!!开始时没有加
193       printf("%I64d\n",sum);
194 }
195 
196 void pripare()
197 {
198     LL i;
199     M_hxl[0].init();
200     for(i=1;i<32;i++)
201     M_hxl[i]=Multiply(M_hxl[i-1],M_hxl[i-1]);
202 }
203 
204 int main()
205 {
206     LL T,n,i;
207     pripare();//打表。
208     while(scanf("%I64d",&T)>0)
209     {
210         for(i=1;i<=T;i++)
211         {
212             scanf("%I64d",&n);
213             printf("Case %I64d: ",i);
214             if(n==0)printf("1\n");
215             else if(n==1) printf("4\n");
216             else if(n==2) printf("9\n");
217             else
218             {
219                 power_sum2(n-1);
220             }
221         }
222     }
223     return 0;
224 }

 

转载于:https://www.cnblogs.com/tom987690183/p/3289399.html