bzoj1951 [Sdoi2010]古代猪文

http://www.elijahqi.win/archives/3185
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

nmd 可能今天学傻了非得写扩展lucas T飞了浪费我时间tmd

普通的CRT博主总是记不住在这里推导一下尝试记住
有同余方程

$$\begin{array}{c} x\equiv c_1\pmod {m_1}\\ x\equiv c_2\pmod {m_2} \\ x\equiv c_3\pmod {m_3}\\ ...\\ x\equiv c_k\pmod {m_k} \end{array}$$
设 $M=m1\times m2\times m3...$
设 $n_i\equiv0(mod\ \frac{M}{m_i})$
设 $n_i\equiv1(mod\ m_i)$
则 $\sum\limits_{i=1}^{k}c_i \times n_i$是一个解 把上面的方程带入即可发现这确实是一个解
设 $n_i=\frac{M}{m_i}x$ $n_i=m_i y+1$
那么即可用上式联立方程求解
$\frac{M}{m_i}x-m_i y=1$ 使用exgcd求出x即我们想要的一个 $x$ 然后再乘上相应的系数就是想要的$n_i$ 然后把所有加起来即可

#include<cstdio>
#include<cctype>
#include<algorithm>
#define ll long long
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f; 
}
int n,g,w[3],mi,prime[5]={0,2,3,4679,35617};
inline int ksm(ll b,int t,int md){static int tmp;
    for (tmp=1;t;t>>=1,b=b*b%md) if (t&1) tmp=tmp*b%md;return tmp;
}
inline 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;
}
inline int inv(int a,int b){static int t1,t2;
    exgcd(a,b,t1,t2);return (t1%b+b)%b;
}
const int N=55617;
int jc[5][N],invg[5][N];
const int md=999911658;
inline int get(int id,int n,int m){
    if (m>n) return 0;
    return (ll)jc[id][n]*invg[id][m]%md*invg[id][n-m]%md;
}
inline int c(int id,int n,int m){static int mod,tmp;mod=prime[id];
    for (tmp=1;n&&m;n/=mod,m/=mod) tmp=(ll)tmp*get(id,n%mod,m%mod)%mod;return tmp;
}
inline void init(int id,int md){
    jc[id][0]=1;for (int i=1;i<md;++i) jc[id][i]=jc[id][i-1]*i%md;invg[id][md-1]=inv(jc[id][md-1],md);
    for (int i=md-2;~i;--i) invg[id][i]=(ll)invg[id][i+1]*(i+1)%md;
}
inline void inc(int &x,int v){x=x+v>=md?x+v-md:x+v;}
int main(){
    freopen("bzoj1951.in","r",stdin);
    n=read();g=read();if (g==md+1) {puts("0");return 0;}
    for (int i=1;i<=4;++i) init(i,prime[i]);int mi=0;
    for (int i=1;i*i<=n;++i){static int tmp;
        if (n%i) continue;tmp=0;
        for (int j=1;j<=4;++j) {
            //printf("%d\n",c(j,n,n/i));
            inc(tmp,(ll)c(j,n,n/i)*(md/prime[j])%md*inv(md/prime[j],prime[j])%md);
        }inc(mi,tmp);
        if (i*i!=n){tmp=0;
            for (int j=1;j<=4;++j) {
            //  printf("%d\n",c(j,n,i));
                inc(tmp,(ll)c(j,n,i)*(md/prime[j])%md*inv(md/prime[j],prime[j])%md);
            }inc(mi,tmp);
        }//printf("%d\n",mi);
    } 
    printf("%d\n",ksm(g,mi,md+1));
    return 0;
}

放博主sb 做法

#include<cstdio>
#include<cctype>
#include<algorithm>
#define ll long long
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f; 
}
//999911658=2*3*4679*35617
int n,g,w[3],mi,prime[5]={0,2,3,4679,35617};
inline int ksm(ll b,int t,int md){static int tmp;
    for (tmp=1;t;t>>=1,b=b*b%md) if (t&1) tmp=tmp*b%md;return tmp;
}
inline 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;
}
inline int inv(int a,int b){static int t1,t2;
    exgcd(a,b,t1,t2);return (t1%b+b)%b;
}
inline int jc(int n,int p){
    if (!n) return 1;ll t=1,tt=n%p;
    for (int i=2;i<p;++i) if (i%p) t=t*i%p;
    t=ksm(t,n/p,p);
    for (int i=2;i<=tt;++i) if (i%p) t=t*i%p;
    t=t*jc(n/p,p)%p;return t;
}
inline int calc(int n,int p){
    static ll t1,t2,nm1;
    t1=jc(n,p);nm1=0;t2=1;
    for (ll i=p;i<=n;i=i*p) nm1+=n/i;
    for (int i=1;i<=2;++i){
        if (!w[i])continue;
        for (ll j=p;j<=w[i];j=j*p) nm1-=w[i]/j;
        t2=t2*jc(w[i],p)%p;
    }t1=t1*inv(t2,p)%p;t1=t1*ksm(p,nm1,p);return t1;
}
const int md=999911658;
inline void inc(int &x,int v){x=x+v>=md?x+v-md:x+v;}

int main(){
    freopen("bzoj1951.in","r",stdin);
    n=read();g=read();if (g==md+1) {puts("0");return 0;}
    for (int i=1;i*i<=n;++i){static int t,tmp;
        if (n%i) continue;
        t=n/i;w[1]=t;w[2]=n-t;tmp=0;
        for (int j=1;j<=4;++j) {
            inc(tmp,(ll)calc(n,prime[j])*(md/prime[j])%md*inv(md/prime[j],prime[j])%md);
        }inc(mi,tmp);
        if (i*i!=n){
            t=i;w[1]=t;w[2]=n-t;tmp=0;
            for (int j=1;j<=4;++j) {
                inc(tmp,(ll)calc(n,prime[j])*(md/prime[j])%md*inv(md/prime[j],prime[j])%md);
            }inc(mi,tmp);
        }
    } 
    printf("%d\n",ksm(g,mi,md+1));
    return 0;
}