1951: [Sdoi2010]古代猪文
Time Limit: 1 Sec Memory Limit: 64 MB
Submit: 1525 Solved: 627
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-1取模的结果。分解质因数，可以发现999911658=2*3*4679*35617。所以我们就可以分别对每个质因数求解，然后用中国剩余定理合并。而对每个质因数求解的问题就是求组合数模一个小质数，预处理+Lucas定理。
P.S.如何O(n)预处理1-n关于p的逆元(n<p)？
我们从小到大依次计算。
对于i，令p=x*i-y(0<y<i)，
即 x*i=y(mod p)。
所以i的逆元inv[i]=x*inv[y]，而inv[y]已经计算出来了。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define maxn 100005
using namespace std;
const int p[5]={2,3,4679,35617,999911658};
int n,m,fac[maxn],inv[maxn];
inline int read()
{
int x=0,f=1;char ch=getchar();
while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
inline ll power(ll x,ll y,ll p)
{
ll ans=1;
for(;y;y>>=1,x=x*x%p) if (y&1) ans=ans*x%p;
return ans;
}
inline void exgcd(ll a,ll b,ll &x,ll &y)
{
if (!b){x=1;y=0;return;}
exgcd(b,a%b,x,y);
int t=x;x=y;y=t-a/b*x;
}
inline ll getinv(ll a,ll b)
{
ll x,y;
exgcd(a,b,x,y);
return (x%b+b)%b;
}
inline int c(int n,int m,int p)
{
if (n<m) return 0;
if (n<p&&m<p) return fac[n]*inv[m]%p*inv[n-m]%p;
return c(n/p,m/p,p)*c(n%p,m%p,p)%p;
}
int main()
{
n=read();m=read();
ll ans=0;
m%=(p[4]+1);
if (!m){printf("0\n");return 0;}
F(i,0,3)
{
fac[0]=1;
F(j,1,p[i]-1) fac[j]=fac[j-1]*j%p[i];
inv[0]=inv[1]=1;
F(j,2,p[i]-1) inv[j]=(p[i]/j+1)*inv[j-p[i]%j]%p[i];
F(j,2,p[i]-1) inv[j]=inv[j-1]*inv[j]%p[i];
ll t1=0,t2=getinv(p[4]/p[i],p[i]);
for(int j=1;j*j<=n;j++) if (n%j==0)
{
(t1+=c(n,j,p[i]))%=p[i];
if (j*j==n) continue;
(t1+=c(n,n/j,p[i]))%=p[i];
}
(ans+=(ll)p[4]/p[i]%p[4]*t2%p[4]*t1%p[4])%=p[4];
}
(ans+=p[4])%=p[4];
cout<<power(m,ans,p[4]+1)<<endl;
}