ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

POJ - 3233 Matrix Power Series

2019-02-16 08:54:00  阅读:253  来源: 互联网

标签:node 3233 Matrix int mt POJ go include define


传送门

 

Solution:

1.矩阵分块 题解在这里

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cmath>
 4 #include<algorithm>
 5 #include<cstring>
 6 #define R register
 7 #define go(i,a,b) for(R int i=a;i<=b;i++)
 8 #define yes(i,a,b) for(R int i=a;i>=b;i--)
 9 #define ll long long
10 #define db double
11 using namespace std;
12 ll n,m;int mod;
13 struct node{ll mt[3][3];};
14 node calc(node x,node y,int a,int b,int c)
15 {
16     node z;memset(z.mt,0,sizeof(z.mt));
17     go(i,1,a)
18         go(j,1,b)
19         go(k,1,c)
20         z.mt[i][j]=(z.mt[i][j]+x.mt[i][k]*y.mt[k][j])%mod;
21     return z;
22 }
23 node a,b,c;
24 void sol()
25 {
26     while(m)
27     {
28         if(m&1) b=calc(b,a,2,2,2);
29         m>>=1;a=calc(a,a,2,2,2);
30     }
31 }
32 int main()
33 {
34     scanf("%lld",&n);m=n-2;mod=1000000007;
35     if(n<=2){printf("1");return 0;}
36     a.mt[1][1]=1,a.mt[1][2]=1,a.mt[2][1]=1,a.mt[2][2]=0;
37     b.mt[1][1]=1,b.mt[1][2]=0,b.mt[2][1]=0,b.mt[2][2]=1;
38     c.mt[1][1]=1;c.mt[2][1]=1;
39     sol();
40     c=calc(b,c,2,1,2);
41     printf("%lld",c.mt[1][1]);
42     return 0;
43 }
View Code

但是这个代码有点问题 就是k=1和k=2的时候会错,是因为递推式子的原因,详见题解

但还是AC了 因为没有k=1,2这么水的数据 这里说一下我就不改了

2.分治 题解

 

没啦qwq

 

标签:node,3233,Matrix,int,mt,POJ,go,include,define
来源: https://www.cnblogs.com/forward777/p/10386687.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有