FZU 1683 纪念SlingShot

Problem 1683 纪念SlingShot

Accept: 682 Submit: 2368
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

236

Sample Output

Case 1: 37

Case 2: 313

思路：构造矩阵。矩阵的快速幂和矩阵的乘法，（注意两个矩阵相乘不满足交换律）

由此我们可以得到两个矩阵 {1 3 2 7

0 3 2 7

0 1 0 0

0 0 1 0 }

{ 9 0 0 0

5 0 0 0

3 0 0 0

1 0 0 0 }

代码：

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#define mod 2009using namespace std;struct node
{int m[5][5];
};node a,I;
int len=4;void init()
{I.m[1][1]=9;I.m[2][1]=5;I.m[3][1]=3;I.m[4][1]=1;a.m[1][1]=1;a.m[1][2]=3;a.m[1][3]=2;a.m[1][4]=7;a.m[2][1]=0;a.m[2][2]=3;a.m[2][3]=2;a.m[2][4]=7;a.m[3][1]=0;a.m[3][2]=1;a.m[3][3]=0;a.m[3][4]=0;a.m[4][1]=0;a.m[4][2]=0;a.m[4][3]=1;a.m[4][4]=0;return ;
}node mut(node c,node b)
{node ans;for(int i=1;i<=4;i++){for(int j=1;j<=4;j++){ans.m[i][j]=0;for(int k=1;k<=4;k++){ans.m[i][j]=(ans.m[i][j]+b.m[i][k]*c.m[k][j])%mod;}}}return ans;
}node solve(int kk)
{node ans=I;node aa=a;while(kk){if(kk&1){ans=mut(ans,aa);}aa=mut(aa,aa);kk/=2;}return ans;
}int main()
{init();int t;scanf("%d",&t);int n;int cas=0;while(t--){scanf("%d",&n);printf("Case %d: ",++cas);if(n<=2){ if(n==0) printf("1\n");else if(n==1) printf("4\n");else printf("10\n");continue;}else{node fin = solve(n-2);printf("%d\n",fin.m[1][1]%mod);}       }return 0;}

FZU 1683 纪念SlingShot相关推荐

FZU - 1683 纪念SlingShot
Problem 1683 纪念SlingShot Accept: 719 Submit: 2492 Time Limit: 1000 mSec Memory Limit : 32768 K ...
FZU 1683 纪念SlingShot (简单的矩阵快速幂)
纪念SlingShot 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,对于给 ...
fzu 1683 纪念SlingShot 矩阵
题意: 对于给出的递推关系,求前n项和思路: 数据大,有规律.矩阵快速幂吧.然而为啥这题long long 就TLE.智障题. 构造: Sn=Sn-1+F(n)= Sn-1 +3F(n-1)+2F( ...
【FOJ】Problem 1683 纪念SlingShot
Problem 1683 纪念SlingShot. 思路构造矩阵T S(n) = Tn-2[0][0] × S(2) + Tn-2[0][1] × f(2) + Tn-2[0][2] × f(1) ...
FOJ 1683 纪念SlingShot
来源:http://acm.fzu.edu.cn/problem.php?pid=1683 概述:已知 F(n)=3 * F(n-1)+2 * F(n-2)+7 * F(n-3),n>=3,其中 ...
纪念SlingShot 矩阵快速幂构造矩阵
纪念SlingShot Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u Submit ...
纪念SlingShot
很多时候题目里面没有给你n的范围,所以把所有的情况都列出,没有列全的话也会出现time limit 的错误!!! 及时给了不明确,也要列出来题意:f(0)=1,f(1)=3,f(2)=5,f(n)= ...
纪念SlingShot 矩阵
第一次独立做矩阵,只WA了一次,好高兴~~ 题意:已知 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,对于给 ...
第一篇博客，纪念一下，终于开通啦！
终于有网了,可以把自己积累的内容都放到博客上来了. 纪念一下. 转载于:https://www.cnblogs.com/tinazzz/p/7067482.html

FZU 1683 纪念SlingShot

Accept: 682 Submit: 2368
Time Limit: 1000 mSec Memory Limit : 32768 KB

Problem Description

Input

Output

Sample Input

Sample Output

思路：构造矩阵。矩阵的快速幂和矩阵的乘法，（注意两个矩阵相乘不满足交换律）

FZU 1683 纪念SlingShot相关推荐

最新文章

热门文章

FZU 1683 纪念SlingShot

Accept: 682 Submit: 2368 Time Limit: 1000 mSec Memory Limit : 32768 KB

Problem Description

Input

Output

Sample Input

Sample Output

思路： 构造矩阵。矩阵的快速幂和矩阵的乘法， （注意两个矩阵相乘不满足交换律）

FZU 1683 纪念SlingShot相关推荐

最新文章

热门文章

Accept: 682 Submit: 2368
Time Limit: 1000 mSec Memory Limit : 32768 KB

思路：构造矩阵。矩阵的快速幂和矩阵的乘法，（注意两个矩阵相乘不满足交换律）