1951: [Sdoi2010]古代猪文

题目：1951: [Sdoi2010]古代猪文

tag：Lucas定理，快速幂，扩展欧几里得，中国剩余定理

第一次做这么恶心的数论题，真心写的想吐了。

题目还是粘一下吧，还是蛮好的一个题呀：

1951: [Sdoi2010]古代猪文

Time Limit: 1 Sec   Memory Limit: 64 MB
Submit: 681   Solved: 236
[ Submit][ Status]

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

[ Submit][ Status]

题意：题目给出 n，g，要求所有满足 x 整除 n 的情况下，C（n，x）的和对 999911659 取余的结果作为幂，然后求g的这么多次幂取余的结果。

思路：由于p=999911659 是个大素数，根据上篇写的指数循环节的相关知识可以得到 （g^y ）%p = （g^（y%phi（p）+phi（p）））%p 同时phi（p）=p-1 带进去得到 （g^y）%p = （g^(y%（p-1）)）% p

  关于由n和约数得到的组合数，用Lucas定理可以解决，跟前几篇写的一样，不过打阶乘表，只需要打到 35617 ，剩下的就不废话了。

  讨论一下p ,  由于当前的 mod = p-1 = 999911658 = 2 * 3 * 4679 * 35617 , 如果当前需要进行取余的数已经求出来了，那么我们直接对mod取余就很容易是溢出的结果来进行取余，所以要把这个取余分成四个同余方程，然后用中国剩余定理合并，就搞定了

#include <cstring>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstdio>
using namespace std;
#define maxn 35620
typedef long long LL;
LL w[4]={2,3,4679,35617};
LL a[4]={0,0,0,0};
LL fac[4][maxn];
LL Pow(LL a,LL b,LL mod)
{
    LL ans=1;
    while(b)
    {
        if(b&1)
        {
            b--;
            ans=(ans*a)%mod;
        }
        else
        {
            b/=2;
            a=(a*a)%mod;
        }
    }
    return ans;
}
void init()
{
    for(int i=0;i<4;i++)
    {
        fac[i][0]=1;
        for(int j=1;j<=w[i];j++)
        {
            fac[i][j]=(fac[i][j-1]*j)%w[i];
        }
    }
}
LL C(LL n,LL m,int x)
{
    if(n<m)
        return 0;
    return (fac[x][n]*Pow((fac[x][n-m]*fac[x][m]),w[x]-2,w[x]))%w[x];
}
LL Lucas(LL n,LL m,int x)
{
    if(m==0)
        return 1;
    return (Lucas(n/w[x],m/w[x],x)*C(n%w[x],m%w[x],x))%w[x];
}
LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    else
    {
        LL ans=exgcd(b,a%b,x,y);
        LL t=x;
        x=y;
        y=t-a/b*y;
        return ans;
    }
}
LL CRT(int r)
{
    LL M=1;
    LL i,d,x0,y0,ans=0;
    for(i=0;i<r;i++)
        M*=w[i];//好像很容易溢出的样子，CRT弱爆了
    for(i=0;i<r;i++)
    {
        d=M/w[i];
        exgcd(d,w[i],x0,y0);
        ans=(ans+d*x0*a[i])%M;
    }
    while(ans<= 0)
       ans+=M;
    return ans;
}
int main()
{
    init();
    LL g,n;
    while(cin>>n>>g)
    {
        g%=999911659;
        memset(a,0,sizeof(a));
        for(int i=1;i*i<=n;i++)
        {
            if(n%i==0)
            {
                LL tmp=n/i;
                for(int j=0;j<4;j++)
                {
                    if(tmp!=i)
                        a[j]=(a[j]+Lucas(n,i,j))%w[j];
                    a[j]=(a[j]+Lucas(n,tmp,j))%w[j];
                }
            }
        }
        cout<<Pow(g,CRT(4),999911659)<<endl;
    }
    return 0;
}