【SDOI2010】古代猪文（数论）

最新推荐文章于 2021-10-31 19:27:58 发布

xyyxyyx

最新推荐文章于 2021-10-31 19:27:58 发布

阅读量149

点赞数

分类专栏：题解

本文链接：https://blog.csdn.net/xyyxyyx/article/details/102482785

版权

题解专栏收录该内容

107 篇文章 3 订阅

订阅专栏

文章目录

题意
思路
注意
代码

题意

给出 $N,G(1\le N,G \le 10^9)$ ，求 $G^{\sum_{k|N}C_N^k}\;mod\;999911659$

思路

根据欧拉定理（或者费马小定理），我们需要求 $\sum_{k|N}C_N^k\;mod\;999911658$ ，而模数是个合数。

一上来就想套Lucas？

慢着。把模数分解质因数： $999911658 = 2 * 3 * 4679 * 35617$ ，和一般拓展卢卡斯不同的是，所有质因数的指数都是1，所以可以直接用卢卡斯定理代替拓展卢卡斯麻烦的阶乘取模。这算是一个小技巧吧，要是下次碰到给定而不是输入的模数，那就尝试这样化简题目。

最后注意一下特判，欧拉定理成立的条件是 $G$ 与模数互质。因为保证N大于0，所以不可能有指数为0的情况，所以只要 $999911659\mid G$ 就输出0。

有这样一组数据999911657 999911659，使得指数在取模之后为0，不特判答案将会是1。

注意

欧拉定理的条件

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int Mod = 999911659, phi = 999911658;
const int mod[4] = {2, 3, 4679, 35617}, N = 35618;
int n, g, ans, fac[4][N], ifac[4][N];

int fpow(int x, int y, int p){
	int ret = 1;
	while (y){
		if (y & 1) ret = 1LL * ret * x % p;
		x = 1LL * x * x % p;
		y >>= 1;
	}
	return ret;
}

int inv(int x, int p){
	return fpow(x, p-2, p);
}

void init(int id, int p)
{
	fac[id][0] = 1;
	for (int i = 1; i < p; ++ i)
		fac[id][i] = 1LL * fac[id][i-1] * i % p;
	ifac[id][p-1] = inv(fac[id][p-1], p);
	for (int i = p-2; i >= 0; -- i)
		ifac[id][i] = 1LL * ifac[id][i+1] * (i+1) % p;
}

int Lucas(int n, int m, int id, int p){
	if (n < p && m < p) return 1LL * fac[id][n] * ifac[id][m] % p * ifac[id][n-m] % p;
	return 1LL * Lucas(n / p, m / p, id, p) * Lucas(n % p, m % p, id, p) % p;
}

int CRT(int a, int b, int p){
	return 1LL * inv(p / b, b) * (p / b) % p * a % p;
}

int solve(int n, int m, int p)
{
	int ret = 0;
	for (int i = 0; i < 4; ++ i)
		(ret += CRT(Lucas(n, m, i, mod[i]), mod[i], phi)) %= p;
//	cout << n << " " << m << " " << ret << endl;
	return ret;
}

int main()
{
	for (int i = 0; i < 4; ++ i)
		init(i, mod[i]);
	cin >> n >> g;
	if (g % Mod == 0){
		puts("0");
		return 0;
	}
	ans = 0;
	for (int i = 1; i * i <= n; ++ i)
		if (n % i == 0){
			(ans += solve(n, i, phi)) %= phi;
			if (i * i != n)
				(ans += solve(n, n / i, phi)) %= phi;
		}
	ans = fpow(g, ans, Mod);
	cout << ans << endl;
	return 0;
}