neu1250矩阵快速幂哪~~

其实呢，自己乱测试数据也是会发现BUG的

一直都会是自己的代码错误！

折腾了我一大晚上到头来竟然矩阵数组的存储也是long long类型

不过终归是解决了稍微学到一点教训吧

一个是写完代码自己可以代几组数据试一下有时候会出现那种很明显的错误比如这个题出现负数神马啦

再一个经验是Long long 真的比int 好用？？

自己的程序几乎每个地方都可能超出范围

再最后是简单的程序也要多练现在自己敲一遍矩阵的乘法

#include<stdio.h>
#include<string.h>
#define da 100000007
struct M
{long long s[3][3];
};
struct M multiply(struct M a,struct M b){struct M c;memset(c.s,0,sizeof(c.s));for(int i=0;i<3;i++)for(int j=0;j<3;j++)for(int k=0;k<3;k++)c.s[i][j]=(c.s[i][j]+(a.s[i][k]*b.s[k][j])%da)%da;return c;
}
struct M paw(struct M a,long long t){if(t==1)return a;else{struct M b=paw(a,t/2);if(t&1){return multiply(multiply(b,b),a);}elsereturn multiply(b,b);}
}
int main()
{struct M a,c;a.s[0][0]=3;a.s[0][1]=2;a.s[0][2]=1;a.s[1][0]=1;a.s[1][1]=0;a.s[1][2]=0;a.s[2][0]=0;a.s[2][1]=1;a.s[2][2]=0;int a1,a2,a3;long long n;while(scanf("%d%d%d%lld",&a1,&a2,&a3,&n)!=EOF){if(n==0) { printf("%d\n",a1%da);}else if(n==1) {printf("%d\n",a2%da);}else if(n==2) { printf("%d\n",a3%da);}else{c=paw(a,n-2);long long x=(c.s[0][0]*(a3%da)+c.s[0][1]*(a2%da)+c.s[0][2]*(a1%da))%da;printf("%lld\n",x);}}return 0;
}

　　按位与“&”功能强大？

for(int i=0;i<N;i++)
        for(int j=0;j<N;j++)
           for(int k=0;k<N;k++)
              {c.s[i][j]+=(a.s[i][k]%da)*(b.s[k][j]%da);
               c.s[i][j]%=da;//因为各对应行列为i,j
              }//也对

以后矩阵快速幂用这个自己写的吧....

转载于:https://www.cnblogs.com/cgf1993/archive/2013/04/12/3017226.html

neu1250矩阵快速幂哪~~相关推荐

矩阵快速幂+构造方法
与快速幂一样,可以将递推式通过二进制的方式来进行优化,这个学了快速幂就是十分容易理解大概的板子如下: struct mat///自己定义大小的矩阵 {ll m[11][11]; }; mat mul ...
【做题】SRM701 Div1 Hard - FibonacciStringSum——数学和式＆矩阵快速幂
原文链接 https://www.cnblogs.com/cly-none/p/SRM701Div1C.html 题意:定义"Fibonacci string"为没有连续1的01串 ...
快速幂 + 矩阵快速幂
快速幂 1 #include<iostream> 2 #include<algorithm> 3 #include<cstring> 4 #define LL lo ...
HDU4549（矩阵快速幂+快速幂）
f(n)=a^f(n-1) + b^f(n-2):计算矩阵部分用矩阵快速幂:计算a的幂次和b的幂次用快速幂. #include<iostream> #include<algorith ...
[HNOI2008]GT考试[矩阵快速幂+kmp优化的dp]
解题思路:假如说我们用f[i]表示长度为i的串能组合成无不吉利数字的组合的个数的话我们无法找到f[i]和f[i+1]的关系,就是我们下一位填某个数字会不会出现不吉利串,这就和你前面的串末尾于不吉利串重 ...
I-Matrix Power Series POJ - 3233 矩阵快速幂+分治
I-Matrix Power Series POJ - 3233 矩阵快速幂+分治 Problem Description Given a n × n matrix A and a positive ...
H - Fibonacci POJ - 3070 （矩阵快速幂）
H - Fibonacci POJ - 3070 (矩阵快速幂) Description In the Fibonacci integer sequence, F0 = 0, F1 = 1, and ...
HDU 6185 Covering 矩阵快速幂递推
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=6185 题目描述: 一个4*n的矩形, 你用1*2的矩形覆盖有多少种方案, n <= 1e18 ...
bzoj 1409 Password 矩阵快速幂+欧拉函数
可以发现,该数组的mi就是斐波那契数列所以要矩阵快速幂搞出第n位但是斐波那契数列上涨的很快,这就需要欧拉定理了 p^phi(q)%q=1(gcd(p,q)==1) p是素数,所以可以用然后需要5 ...

neu1250矩阵快速幂哪~~

neu1250矩阵快速幂哪~~相关推荐

最新文章

热门文章