1. 程式人生 > >[SDOI2010]古代豬文

[SDOI2010]古代豬文

putc form algo stdin program main 圖書 很多 質數

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;
}

[SDOI2010]古代豬文