bzoj1954 你是我孙子定理

Description

“在那山的那边海的那边有一群小肥猪。他们活泼又聪明，他们调皮又灵敏。他们自由自在生活在那绿色的大草坪，他们善良勇敢相互都关心……” ——选自猪王国民歌 很久很久以前，在山的那边海的那边的某片风水宝地曾经存在过一个猪王国。猪王国地理位置偏僻，实施的是适应当时社会的自给自足的庄园经济，很少与外界联系，商贸活动就更少了。因此也很少有其他动物知道这样一个王国。 猪王国虽然不大，但是土地肥沃，屋舍俨然。如果一定要拿什么与之相比的话，那就只能是东晋陶渊明笔下的大家想象中的桃花源了。猪王勤政爱民，猪民安居乐业，邻里和睦相处，国家秩序井然，经济欣欣向荣，社会和谐稳定。和谐的社会带给猪民们对工作火红的热情和对未来的粉色的憧憬。 小猪iPig是猪王国的一个很普通的公民。小猪今年10岁了，在大肥猪学校上小学三年级。和大多数猪一样，他不是很聪明，因此经常遇到很多或者稀奇古怪或者旁人看来轻而易举的事情令他大伤脑筋。小猪后来参加了全猪信息学奥林匹克竞赛(Pig Olympiad in Informatics, POI)，取得了不错的名次，最终保送进入了猪王国大学(Pig Kingdom University, PKU)深造。 现在的小猪已经能用计算机解决简单的问题了，比如能用P++语言编写程序计算出A + B的值。这个“成就”已经成为了他津津乐道的话题。当然，不明真相的同学们也开始对他刮目相看啦~ 小猪的故事就将从此展开，伴随大家两天时间，希望大家能够喜欢小猪。 题目描述 猪王国的文明源远流长，博大精深。 iPig在大肥猪学校图书馆中查阅资料，得知远古时期猪文文字总个数为N。当然，一种语言如果字数很多，字典也相应会很大。当时的猪王国国王考虑到如果修一本字典，规模有可能远远超过康熙字典，花费的猪力、物力将难以估量。故考虑再三没有进行这一项劳猪伤财之举。当然，猪王国的文字后来随着历史变迁逐渐进行了简化，去掉了一些不常用的字。 iPig打算研究古时某个朝代的猪文文字。根据相关文献记载，那个朝代流传的猪文文字恰好为远古时期的k分之一，其中k是N的一个正约数（可以是1和N）。不过具体是哪k分之一，以及k是多少，由于历史过于久远，已经无从考证了。 iPig觉得只要符合文献，每一种能整除N的k都是有可能的。他打算考虑到所有可能的k。显然当k等于某个定值时，该朝的猪文文字个数为N / k。然而从N个文字中保留下N / k个的情况也是相当多的。iPig预计，如果所有可能的k的所有情况数加起来为P的话，那么他研究古代文字的代价将会是G的P次方。 现在他想知道猪王国研究古代文字的代价是多少。由于iPig觉得这个数字可能是天文数字，所以你只需要告诉他答案除以999911659的余数就可以了。

Input

有且仅有一行：两个数N、G，用一个空格分开。

Output

有且仅有一行：一个数，表示答案除以999911659的余数。

Sample Input

4 2

Sample Output

2048

HINT

10%的数据中，1 <= N <= 50；
20%的数据中，1 <= N <= 1000；
40%的数据中，1 <= N <= 100000；
100%的数据中，1 <= G <= 1000000000，1 <= N <= 1000000000。

Source

记k=sigma（C（n,n/i）(其中i是n的因数)），则答案为G^k。

这个k可能很大啊，但是考虑mod数是一个质数，所以G^k=G^(k mod (P - 1))但是这里要注意特判G=P的情况（答案是0）

根据中国剩余定理，对P-1分解质因数，C（n,m）=n!/m!/(n-m)!，然后考虑到n可能特别大，所以可以用迭代的方法求出一个数的阶乘共有多少个因数p以及在mod p1（记p1=p^q）意义下除去p的倍数的数的积，这可以用二元组表示，求有多少个因数p很简单，假设当前阶乘是n，那么n/p数可以被p整除，但是考虑到一个数可能是p的若干次方，所以我们把n/p继续进行计算，第二个东西的话我们可以看出每p1个数在mod p1下是同余的，所以在迭代的时候用快速幂计算即可。

然后就可以用中国剩余定理合并，对于每个p1求出的答案f[i]，sigma（P/p1 *（P/p1）^-1(在mod p1的意义下) * f[i]）就是答案。

所以就枚举因数啦然后搞搞就吼了。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int P = 999911659;
const int MOD = P - 1;
int i, j, k, p[20], p1[20], G, n, up, bb[5][100005], tt;
struct sb{
	int x, y;
};
inline int ksm(int x, int y, int z)
{
	int b = 1;
	while (y)
	{
		if (y & 1) b = (long long)b * x % z;
		x = (long long)x * x % z;
		y >>= 1;
	}
	return b;
}
inline sb mul(sb a, sb b, int x)
{
	return (sb){a.x * (long long)b.x % p1[x], a.y + b.y};
}
inline int Inv(int a, int x)
{
	return ksm(a, p1[x] / p[x] * (p[x] - 1) - 1, p1[x]);
}
inline sb div(sb a, sb b, int x)
{
	return (sb){a.x * (long long)Inv(b.x, x) % p1[x], a.y - b.y};
}
inline int change(sb a, int x)
{
	return a.x * (long long)ksm(p[x], a.y, p1[x]) % p1[x];
}
inline int dfs(int n, int x)
{
	if (n < p[x]) return bb[x][n];
	int tot;
	if (n > p1[x])
	{
		int tim = n / p1[x];
		tot = ksm(bb[x][p1[x]], tim, p1[x]);
		tot = (long long)tot * bb[x][n - tim * p1[x]] % p1[x];
	}
	else tot = bb[x][n];
	return tot * (long long)dfs(n / p[x], x) % p1[x];
}
inline int dfs1(int n, int x)
{
	if (n < p[x]) return 0;
	else return n / p[x] + dfs1(n / p[x], x);
}
inline sb divide(int n, int x)
{
	sb h;
	h.x = dfs(n, x);
	h.y = dfs1(n, x);
	return h;
}
inline int crt(int n, int m, int x)
{
	sb a = divide(n, x), b = divide(m, x), c = divide(n - m, x);
	return change(div(a, mul(b, c, x), x), x);
}
inline int solve(int n, int m)
{
	int tot = 0;
	for(int i = 1; i <= tt; i ++)
		tot = (tot + (long long)P / p1[i] * Inv(P / p1[i], i) % MOD * crt(n, m, i)) % MOD;
	return tot;
}
int main()
{
	k = P - 1;
	for(i = 2; i * i <= P; i ++)
		if (k % i == 0)
		{
			p[++tt] = i;
			p1[tt] = 1;
			while (k % i == 0)
			{
				k /= i;
				p1[tt] *= i;
			}
		}
	if (k > 1) p1[++tt] = p[tt] = k;
	cin >> n >> G;
	if (G == P) {cout << 0; return 0;}
	for(j = 1; j <= tt; j ++)
	{
		bb[j][0] = 1;
		for(i = 1; i <= p1[j]; i ++)
			if (i % p[j] == 0) bb[j][i] = bb[j][i - 1];
			else bb[j][i] = bb[j][i - 1] * (long long)i % p1[j];
	}
	for(i = 1; i * i <= n; i ++)
		if (n % i == 0)
		{
			j = n / i;
			up += solve(n, n / i);
			if (up >= MOD) up -= MOD;
			if (i != j)
			{
				up += solve(n, n / j);
				if (up >= MOD) up -= MOD;
			}
		}
	cout << ksm(G, up, P) << endl;
}