bzoj 1951 [Sdoi2010]古代猪文（数论知识）

标签：

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。

【思路】

一道优(e)秀(xin)的数论题。

首先我们要求的是(G^sigma{ C(n,n/i),i|n })%P，即G^M %P，根据费马小定理G^(M-1) 1(mod P),我们要求的就是G^(M%(P-1)) %P。

考虑C(n,i)%(P-1)，由于n i P都比较大所以不好求组合数。发现P-1可以分解质因数为2,3,4679,35617，将C(n,i)对每一个质因子取模，会得到一个形为x ai(mod pi)的模线性方程组，可以用中国剩余定理确定x。对于C(n,i)%p，此时p比较小，我们可以用lucas定理求解。

总的来说就是先用O(sqrt(n))的时间枚举约数，然后用lucas定理求出不同模数下的ai，最后联立方程组，中国剩余定理解。

注意当G=P的时候费马小定理不成立，此时答案为0。

关于lucas的写法，a^p-2 %p是a在模p下的逆，因为a^(p-2) *a=a^(p-1)，由费马小定理得a^(p-1) %p=1。

【定理】

1— 费马小定理

假如p是质数，且(a,p)=1，那么 a^(p-1)≡1（mod p）

2— Lucas定理

C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)

【代码】

 1 #include<cmath>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<iostream>
 5 using namespace std; 
 6 
 7 typedef long long LL;
 8 const int N = 36000;
 9 const LL mod[5]={2,3,4679,35617,999911659}; 
10 
11 LL fac[4][N],a[4],n,G;
12 
13 void gcd(LL a,LL b,LL &x,LL& y) {
14     if(!b) { x=1;y=0; }
15     else gcd(b,a%b,y,x) , y-=x*(a/b);
16 }
17 LL pow(LL x,LL p,LL MOD) {
18     LL ans=1;
19     while(p) {
20         if(p&1) ans=(ans*x)%MOD;
21         x=(x*x)%MOD; p>>=1;
22     }
23     return ans;
24 }
25 LL C(LL n,LL m,int x) {
26     if(n<m) return 0;
27     return (fac[x][n]*pow(fac[x][n-m]*fac[x][m],mod[x]-2,mod[x]))%mod[x];
28 }
29 LL lucas(LL n,LL m,int x) {
30     if(!m) return 1;
31     return lucas(n/mod[x],m/mod[x],x)*C(n%mod[x],m%mod[x],x)%mod[x];
32 }
33 LL china() {
34     LL M=1,d,x,y,ans=0;
35     for(int i=0;i<4;i++) M*=mod[i];
36     for(int i=0;i<4;i++) {
37         d=M/mod[i];
38         gcd(d,mod[i],x,y);
39         ans=(ans+d*x*a[i])%M;
40     }
41     while(ans<=0) ans+=M;
42     ans=(ans+M)%M;
43     return ans;
44 }
45 void get_fac() {
46     for(int i=0;i<4;i++) {
47         fac[i][0]=1;
48         for(int j=1;j<=mod[i];j++)
49             fac[i][j]=(fac[i][j-1]*j)%mod[i];
50     }
51 }
52 int main() {
53     get_fac();
54     scanf("%lld%lld",&n,&G);
55     G%=mod[4];
56     if(!G) { puts("0"); return 0; }
57     for(int i=1;i*i<=n;i++) if(n%i==0) {
58         LL tmp=n/i;
59         for(int j=0;j<4;j++) {
60             a[j]=(a[j]+lucas(n,tmp,j))%mod[j];
61             if(tmp!=i) a[j]=(a[j]+lucas(n,i,j))%mod[j];
62         }
63     }
64     printf("%lld",pow(G,china(),mod[4]));
65     
66     return 0;
67 }

bzoj 1951 [Sdoi2010]古代猪文（数论知识）

标签：

原文地址：http://www.cnblogs.com/lidaxin/p/5223344.html