ICode9

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

[NOI2013]矩阵游戏(数列通项+费马小定理)

2020-09-22 22:31:22  阅读:308  来源: 互联网

标签:long 费马 int ret static 通项 MAXL NOI2013 mod


题目链接

Analysis

先把单独一行拿出来看,设 \(f_1\) 是这一行的第一个元素,有 \(f_i=f_{i-1}*a+b\)。所以 \(f_m=f_1a^{m-1}+\frac{a^{i-1}-1}{a-1}b\)。如果不会的可以再去补一下高中数学。

然后设 \(g_i\) 是第 \(i\) 行的 \(f_m\),有 \(g_i=(g_{i-1}c+d)a^{m-1}+\frac{a^{i-1}-1}{a-1}b\),然后换个元又变成上面的式子,搞一搞就出来了。

但是 \(n,m\) 太大怎么搞?我们有一个费马小定理,\(a^{p-1}=1\pmod p,a<p\)。然后就可以降到 \(p\) 以下了。

坑点:注意 \(a=1\) 时等比数列求和公式不存在,需要特判,而次时又需要模 \(p\) 的 \(n,m\),所以 \(n,m\) 两个都要模。

namespace Solve{
	const long long mod = 1e9 + 7;
	const int MAXL = 1000010;
	static char n[MAXL], m[MAXL];
	static int lenn, lenm;
	static long long a, b, c, d;
	long long ksm(long long x, long long y) {
		long long ret = 1;
		while (y) {
			if (y & 1) ret = (ret * x) % mod;
			x = (x * x) % mod;
			y >>= 1;
		}
		return ret;
	}
	void BruteForcePlus() {
		long long MOD = mod - 1;
		long long nn = 0, mm = 0;
		long long nnn = 0, mmm = 0;
		for (int i = 1; i <= lenn; i++) nn = (nn * 10 + n[i] - '0') % MOD, nnn = (nnn * 10 + n[i] - '0') % mod;
		for (int i = 1; i <= lenm; i++) mm = (mm * 10 + m[i] - '0') % MOD, mmm = (mmm * 10 + m[i] - '0') % mod;
		if (nn == 0) nn = MOD;
		if (mm == 0) mm = MOD;
		long long y = c * ksm(a, mm - 1) % mod;
		long long x = ((ksm(a, mm - 1) - 1 + mod) * ksm(a - 1, mod - 2) % mod * b % mod + d * ksm(a, mm - 1) % mod) % mod;
		long long ans = (ksm(a, mm - 1) + (ksm(a, mm - 1) - 1 + mod) * ksm(a - 1, mod - 2) % mod * b % mod) % mod;
		if (a == 1) {
			x = ((mmm - 1) * b % mod + d + mod) % mod;
			ans = (1 + (mmm - 1) * b % mod + mod) % mod;
		}
		if (y == 1) ans = (ans + x * (nnn - 1) % mod + mod) % mod;
		else ans = (ans * ksm(y, nn - 1) % mod + (ksm(y, nn - 1) - 1 + mod) * ksm(y - 1, mod - 2) % mod * x % mod) % mod;
		print(ans);
	}
	void MAIN() {
		scanf("%s%s", n + 1, m + 1);
		scanf("%lld%lld%lld%lld", &a, &b, &c, &d);
		lenn = strlen(n + 1);
		lenm = strlen(m + 1);
		BruteForcePlus();
	}
} using namespace Solve;

标签:long,费马,int,ret,static,通项,MAXL,NOI2013,mod
来源: https://www.cnblogs.com/zcr-blog/p/13715290.html

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

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

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

ICode9版权所有