NOI 2013 矩阵游戏

最新推荐文章于 2022-10-29 23:30:14 发布

雪融了就该花开

最新推荐文章于 2022-10-29 23:30:14 发布

阅读量120

点赞数 2

分类专栏：计算机算法文章标签：矩阵线性代数算法抽象代数

本文链接：https://blog.csdn.net/weixin_38839688/article/details/118212294

版权

计算机算法专栏收录该内容

15 篇文章 1 订阅

订阅专栏

原题传送门：https://uoj.ac/problem/124

大致题目意思：

其中

于是，网上便出现了名为矩阵费马小定理的说法，我一头雾水证明不出来，

只是花了一个多小时，才证明出了一个递推表达式符合费马小定理的情况，证明方法如下：

于是，默认可以认为矩阵乘法也符合欧拉定理（大嘘）

所以我便写了一手取模，矩阵快速幂，通过了所有的比赛时用例，皆大欢喜。

代码如下：

#include "bits/stdc++.h"
typedef long long LL;
using namespace std;
LL a, b, c, d, n1, m1;
string n,m;
LL MOD = 1e9 + 7;
struct M{
	LL v[2][2];
	M(){
		memset(v,0,sizeof(v));
	}
	friend M operator*(M a,M b){
		M c;
		for(int i=0;i<2;i++)
			for(int j=0;j<2;j++)
				for(int k=0;k<2;k++)
					c.v[i][j]=(c.v[i][j]+a.v[i][k]*b.v[k][j] % MOD)%MOD;
		return c;
	}
	friend M operator^(M a,LL b){
		M c;c.v[0][0]=c.v[1][1]=1;
		for(LL i=b;i;i>>=1,a=a*a)
			if(i&1)c=c*a;
		return c;
	}
}A,B,f;

LL divMod(string s, LL mod)
{
    LL res = 0;
    for (int i = 0; i < s.length(); i++) {
        res = (s[i] - '0') + res * 10;
        res = res % mod;
    }
    return res;
}

int main()
{
    cin >> n >> m >> a >> b >> c >> d;
    LL zq = MOD - 1;
    if (c == 1 && a == 1) n1 = (divMod(n, zq + 1) + zq) % (zq + 1);
        else n1 = (divMod(n, zq) + zq - 1) % zq;
    if (a == 1) m1 = (divMod(m, zq + 1) + zq) % (zq + 1);
        else m1 = (divMod(m, zq) + zq - 1) % zq;
	A.v[0][0]=1;A.v[0][1]=b;A.v[1][1]=a;
	B.v[0][0]=1;B.v[0][1]=d;B.v[1][1]=c;
	f.v[0][0]=f.v[0][1]=1;
	A=A^m1;B=A*B;B=B^n1;
	printf("%lld\n",(f*B*A).v[0][1]);

}

好的，实际上，这个代码过不了UOJ的附加Test数据。。。所以5 X 20 - 3 = 97分。

所以以上全部都是废话。

正确做法，直接对大整数做十进制的快速幂，没了。。。没了。。。蚌埠住了

贴一份其他大佬的AC代码，爬了

https://uoj.ac/submission/479717