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。

Source

Sdoi2010 Contest2

 

传说中的数论板子大合集，题面是真的好笑；

费马小定理+CRT+lucas+快速幂+逆元+组合数，大概也就这么点吧；

题目要求mod p；

考虑幂很大我们用欧拉降幂：
a^b mod p =a^(b%phi(p)) mod p;

使用的条件是a,p互质且p为质数，这个是用费马小定理搞的；

当a,p互质且p为质数，a^p mod p=a，所以a^(p-1) mod p=1;

所以a^b mod p= a^(b mod (p-1)) mod p;

用欧拉降幂的话phi(p)=p-1,所以是一样的；

然后问题的关键在于求：

n有1e9,所以要用lucas求大组合数：

lucas定理：

然后模数不是质数，那么我们考虑常用的方法对每个质因子求一遍然后最后用CRT合并；

因为只有一次询问，那么我们用sqrt(n)的复杂度枚举n的每个因子，加到四个模线性方程组中，最后用CRT算出来；

记得注意各个部分的模数不同，lucas和算组合数的时候用prime，CRT的时候用p-1,最后快速幂的时候用p;

//MADE BY QT666
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=100050;
const int Mod=999911659;
int prime[]={2,3,4679,35617};
int n,g;
ll Ans,sum[N];
ll qpow(ll x,ll y,ll mo){
    ll ret=1;
    while(y){
	if(y&1) (ret*=x)%=mo;
	(x*=x)%=mo;y>>=1;
    }
    return ret;
}
ll jc[10][N];
ll C(int n,int m,int x){
    if(n<m) return 0;
    return jc[x][n]*qpow(jc[x][m],prime[x]-2,prime[x])%prime[x]*qpow(jc[x][n-m],prime[x]-2,prime[x])%prime[x];
}
ll lucas(int n,int m,int x){
    if(m==0) return 1;
    return C(n%prime[x],m%prime[x],x)*lucas(n/prime[x],m/prime[x],x)%prime[x];
}
void Crt(){
    for(int i=0;i<4;i++){
	ll M=(Mod-1)/prime[i];
	(Ans+=sum[i]*M%(Mod-1)*qpow(M,prime[i]-2,prime[i])%(Mod-1))%=(Mod-1);
    }
}
void work(int d){
    for(int i=0;i<4;i++) (sum[i]+=lucas(n,d,i))%=prime[i];
}
int main(){
    scanf("%d%d",&n,&g);if(g==999911659) {puts("0");return 0;}
    for(int i=0;i<4;i++){
	jc[i][0]=1;
	for(int j=1;j<=50000;j++) jc[i][j]=(jc[i][j-1]*j)%prime[i];
    }
    g%=Mod;
    for(int i=1;i*i<=n;i++){
	if(n%i==0){
	    if(i*i==n) work(i);
	    else work(i),work(n/i);
	}
    }
    Crt();
    printf("%lld\n",qpow(g,Ans,Mod));
    return 0;
}