[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\leqslant N\leqslant 50\)
20%的數據中，\(1\leqslant N\leqslant 10^3\)
40%的數據中，\(1\leqslant N\leqslant 10^5\)
100%的數據中，\(1\leqslant G\leqslant 10^9,1\leqslant N\leqslant 10^9\)

---

題目要求\(G^{\sum\limits_{d|N}\binom{N}{d}}\)

可是直接求顯然很差求，注意到答案%999911659，因此指數上須要%999911658

\(999911658=2\times3\times4679\times35617\)

而後咱們就能夠分別對每一個質數算出答案，而後用crt合併便可

/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
typedef long long ll;
typedef long double ld;
typedef unsigned int ui;
typedef unsigned long long ull;
template<typename T>inline T min(T x,T y){return x<y?x:y;}
template<typename T>inline T max(T x,T y){return x>y?x:y;}
inline char gc(){
    static char buf[1000000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
template<typename T>inline T frd(T x){
    int f=1; char ch=gc();
    for (;ch<'0'||ch>'9';ch=gc())   if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=gc()) x=(x<<1)+(x<<3)+ch-'0';
    return x*f;
}
template<typename T>inline T read(T x){
    int f=1;char ch=getchar();
    for (;ch<'0'||ch>'9';ch=getchar())  if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=getchar())    x=(x<<1)+(x<<3)+ch-'0';
    return x*f;
}
inline void print(int x){
    if (x<0)    putchar('-'),x=-x;
    if (x>9)    print(x/10);
    putchar(x%10+'0');
}
namespace Math{
    int mlt(int a,int b,int p){
        int res=1;
        for (;b;b>>=1,a=1ll*a*a%p)  if (b&1)    res=1ll*res*a%p;
        return res;
    }
    int gcd(int a,int b){return !b?a:gcd(b,a%b);}
    void exgcd(int a,int b,int &x,int &y){
        if (!b){x=1,y=0;return;}
        exgcd(b,a%b,x,y);
        int t=x; x=y,y=t-a/b*y;
    }
    int Ex_GCD(int a,int b,int c){
        int d=gcd(a,b),x,y;
        if (c%d)    return -1;
        a/=d,b/=d,c/=d;
        exgcd(a,b,x,y);
        x=(1ll*x*c%b+b)%b;
        return x;
    }
}
const int N=4e4,Mod=999911659;
int fac[N+10],inv[N+10];
void prepare(int p){
    fac[0]=inv[0]=inv[1]=1;
    for (int i=1;i<p;i++)   fac[i]=1ll*fac[i-1]*i%p;
    for (int i=2;i<p;i++)   inv[i]=1ll*(p-p/i)*inv[p%i]%p;
    for (int i=1;i<p;i++)   inv[i]=1ll*inv[i-1]*inv[i]%p;
}
int C(int n,int m,int p){
    if (n<m)    return 0;
    if (n<p&&m<p)   return 1ll*fac[n]*inv[m]%p*inv[n-m]%p;
    return 1ll*C(n/p,m/p,p)*C(n%p,m%p,p)%p;
}
int work(int n,int p){
    prepare(p); int res=0;
    for (int i=1;i*i<=n;i++){
        if (n%i)    continue;
        res=(res+C(n,i,p))%p;
        if (i*i!=n) res=(res+C(n,n/i,p))%p;
    }
    return res;
}  
int main(){
    int N=read(0),G=read(0),_P=Mod-1,Ans=0;
    static int P[4]={2,3,4679,35617};
    for (int i=0;i<4;i++)   Ans=(Ans+1ll*Math::Ex_GCD(_P/P[i],P[i],1)*work(N,P[i])%_P*(_P/P[i])%_P)%_P;
    printf("%d\n",Ans?Math::mlt(G,Ans,Mod):0);
    return 0;
}