1951: [Sdoi2010]古代猪文
Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 2988 Solved: 1275
[Submit][Status][Discuss]
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
分析：题目很长，大意就是求
.
对式子整体进行观察，求出次数后用快速幂处理一下就好了，关键就是次数非常大，而且不好处理.对于乘方而言，缩底数就用快速幂取模，缩指数就用欧拉定理.φ(999911659) = 999911658 = 2*3*4679*35617.设F(n)为次数那一大坨东西，那么最后的答案就是：
G^(F(N) % φ(999911659) + φ(999911659)) % 999911659.
指数的范围是缩小了，但是还是不能快速地求解，我们可以在根号时间内枚举约数，但是直接一个个地计算组合数会因为n,k非常大而T掉，需要缩小n,k的范围.这就要用到Lucas定理.但是模数φ(999911659)也很大，那就把它拆成2*3*4679*35617.运用Lucas定理算出对每一个质数取模得到的答案，就得到了一个同余方程组，利用中国剩余定理合并一下就能得到F(N)了.
实现要先预处理一下模质数p下的阶乘和逆元和逆元的阶乘.
主要的思想还是将数变小，每个部分都有对应的方法，还有就是要注意细节.
#include <cmath> #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; typedef long long ll; const ll mod = 999911659; ll prime[5] = { 0, 2, 3, 4679, 35617 }, jie[5][40000], ni[5][40000], nijie[5][40000], n, g, ans, anss[5]; void init() { for (int i = 1; i <= 4; i++) { ll p = prime[i]; jie[i][0] = 1; jie[i][1] = 1; ni[i][1] = 1; nijie[i][0] = 1; nijie[i][1] = 1; for (ll j = 2; j < 40000; j++) { jie[i][j] = (jie[i][j - 1] * j) % p; ni[i][j] = (p - p / j) * ni[i][p % j] % p; nijie[i][j] = (nijie[i][j - 1] * ni[i][j]) % p; } } } ll qpow(ll a, ll b) { a %= mod; ll res = 1; while (b) { if (b & 1) res = (res * a) % mod; a = (a * a) % mod; b >>= 1; } return res % mod; } ll Lucas(ll pos,ll a, ll b, ll p) { if (a < b) return 0; if (a < p && b < p) return jie[pos][a] * nijie[pos][b] % p * nijie[pos][a - b] % p; return Lucas(pos, a % p, b % p, p) * Lucas(pos, a / p, b / p, p) % p; } void solve(ll x) { for (int i = 1; i <= 4; i++) anss[i] = (anss[i] + Lucas(i,n, x, prime[i])) % prime[i]; } ll exgcd(ll a, ll b, ll &x, ll &y) { if (!b) { x = 1; y = 0; return a; } ll temp = exgcd(b, a % b, x, y), t = x; x = y; y = t - (a / b) * y; return temp; } ll getni(ll x, ll y) { ll tx, ty; exgcd(x, y, tx, ty); return tx % y; } void China() { ll temp = 1; for (int i = 1; i <= 4; i++) temp *= prime[i]; for (ll i = 1; i <= 4; i++) { ll mi = temp / prime[i]; ll niyuan = getni(mi, prime[i]); ans += anss[i] * mi % temp * niyuan % temp; ans %= temp; } ans = (ans + temp) % temp; } int main() { scanf("%lld%lld", &n, &g); init(); for (ll i = 1; i * i <= n; i++) { if (n % i == 0) { if (i * i == n) solve(i); else { solve(i); solve(n / i); } } } China(); printf("%lld ", qpow(g, prime[1] * prime[2] * prime[3] * prime[4] + ans)); return 0; }