bzoj1951 [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

题解

这题是我见过的最恶心的数论题，没有之一

HINT：
Lucas定理，扩展欧几里得定理，中国剩余定理，费马小定理，欧拉定理

题目让求g^{$\sum C_n^{k(k|n)}$}%99911659的答案，令p=999911659，由费马小定理，可将原始化简为g^{$\sum C_n^{k(k|n)}$%(p-1)} %p。
其中，p-1=999911658，是一个巨大的合数，但是我们可以将其拆为2、3、4679、35617四个质数的乘积，然后对四个质数分别暴力枚举n所有的因子，用Lucas定理计算组合数的值，最后用中国剩余定理合并答案。

技巧：
对于取模数为质数的情况，我们可以根据费马小定理，用快速幂求出逆元，应当比用扩展欧几里得快一些
暴力枚举n的因子的时候，只需枚举到$\sqrt n$即可
然而蒟蒻并不懂if(j*j<n)这一句的重要性，虽然去掉后提交仍然正确，但是那样的话连样例都过不了QAQ
还请各位神犇指出原因，本蒟蒻将不胜感激QwQ
UPD：防止$\sqrt n$为整数时算重。

ps：感觉这道题自己的时间效率不怎么样诶，耗时是许多人的n倍，不开森

CODE：

#include<cmath>
#include<cstdio>
typedef long long ll;
const ll M=999911659;
ll prime[5]={0,2,3,4679,35617};
ll fac[40000];
ll inv[40000];
ll a[5];
ll n,g,ans,x,y;
inline ll pow(ll a,ll b,ll p)
{
    ll ans=1;
    a%=p;
    for(;b;b>>=1,a=1ll*a*a%p)
      if(b&1) ans=1ll*ans*a%p;
    return ans;
}
void exgcd(ll a,ll b,ll &x,ll &y)
{
    if(!b) {x=1,y=0;return;}
    exgcd(b,a%b,y,x);
    y-=(a/b)*x;
}
inline ll lucas(ll n,ll m,ll p)
{
    ll ans=1;
    while(n&&m)
    {
        ll a=n%p,b=m%p;
        if(a<b) return 0;
        ans=(1ll*ans*fac[a]*inv[fac[b]]*inv[fac[a-b]])%p;
        n/=p,m/=p;
    }
    return ans;
}
int main()
{
    scanf("%lld%lld",&n,&g);
    inv[0]=1;fac[0]=1;fac[1]=1;
    for(ll i=1;i<=4;i++)
    {
        for(ll j=1;j<=prime[i];j++)
        {
            fac[j]=(fac[j-1]*j)%prime[i];
            inv[j]=pow(j,prime[i]-2,prime[i]);
        }
        for(ll j=1;j<=sqrt(n);j++)
          if(!(n%j))
          {
            a[i]=(1ll*a[i]+lucas(n,j,prime[i]))%prime[i];
            if(j*j<n) a[i]=(1ll*a[i]+lucas(n,n/j,prime[i]))%prime[i];
          }
    }
    for(ll i=1;i<=4;i++)
    {
        exgcd((M-1)/prime[i],prime[i],x,y);
        ans=(1ll*ans+1ll*a[i]*x*((M-1)/prime[i]))%(M-1);
    }
    ans=(ans%(M-1)+M-1)%(M-1);
    if(!ans) ans=M-1;
    printf("%lld",pow(g,ans,M));
    return 0;
}