P1951: [Sdoi2010]古代猪文

呜啊啊啊啊，选错了题，原以为很简单的优化+剪枝就能过结果牵扯到了一堆数论知识。我的错，贴上我的代码（已经尽量优化了）

 const maxn=;
 var n,g,i,j,ans:longint;
 tem:int64;
 function pow(g,x:longint):longint;
 var
 now,tem,i,t:int64;
 begin
   tem:=; now:=g;
   if g= then exit();
   while now<maxn do
     begin
       now:=now*now;
       tem:=tem*;
     end;
   now:=now mod maxn;
   if (now=) and (x div tem>=) then exit();
   if x div tem>= then t:=pow(now,x div tem)
     else t:=;
   now:=g; x:=x-tem*(x div tem);
   while x<> do
     begin
       if (x mod )= then begin
         t:=(t*now) mod maxn;
       end;
       now:=(now*now) mod maxn;
       x:=x div ;
     end;
   exit(t mod maxn);
 end;
 function c(x:longint):longint;
 var i:longint;
 tem:int64;
 begin
   tem:=g;
   for i:= to x do
     tem:=pow(tem,n-i+) mod maxn;
   for i:= to x do tem:=trunc(exp(ln(tem)/i));
   exit(tem);
 end;
 begin
   readln(n,g);
   tem:=;
   for i:= to trunc(sqrt(n)) do
     if n mod i= then
       begin
         if i*i<>n then tem:=tem*((c(i) mod maxn)*(c(n div i) mod maxn) mod maxn) mod maxn
           else tem:=(tem*c(i)) mod maxn;
       end;
   writeln(tem mod maxn);
 end.

然而正确的是。。

#include <cstring>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstdio>
using namespace std;
#define maxn 35620
typedef long long LL;
#define MOD 999911659
#define M 999911658
LL w[]={,,,};
LL a[];
LL fac[][maxn];
LL Pow(LL a,LL b,LL mod)//快速幂
{
    LL ans=;
    while(b)
    {
        if(b&)  ans=(ans*a)%mod;
        a=(a*a)%mod;
        b >>= ;
    }
    return ans;
}
void init()//预处理阶乘
{
    for(int i=;i<;i++)
    {
        fac[i][]=;
        for(int j=;j<=w[i];j++)
        {
            fac[i][j]=(fac[i][j-]*j)%w[i];
        }
    }
}
/********************************
*  组合数取模用费马小定理
*********************************/
LL C(LL n,LL m,int x)//组合数取模
{
    if(n < m) return ;
    return (fac[x][n] * Pow((fac[x][n-m]*fac[x][m]),w[x]-,w[x]))%w[x];
}
/*******************************
*    lucas 处理大组合数取模
********************************/
LL Lucas(LL n,LL m,int x)//lucas定理
{
    if(m==) return ;
    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){x = ;y = ;return a;}

    LL ans = exgcd(b,a%b,x,y);

    LL t = x;x = y;y = t - a/b*y;

    return ans;
}
/***************************************************************************************
*      中国剩余定理：
*       x = b1 % m1
*       x = b2 % m2
*       x = b3 % m3
*       .
*       gcd(m1,m2,m3,...) = 1;
*       M = m1 * m2 * m3 *.....;
*       M1 = m2 * m3 * ...., M2 = m1 * m3 * ...., M3 = m1 * m2 * m4 *....., ......;
*       M1 * M(-1) = 1 % M ,M2 * M2(-1) = 1 % M;
*       res = (M1(-1)*b1 + M2(-1)*b2+.....)%M;res即为所求值
*       注：如果取模的值相同：都是m1  那么 bn的值可以相加计算；
*       略屌。。。。。。。。。。。。。。
*
*****************************************************************************************/
LL CRT()//计算组合数和取模之后的值
{
    LL i,d,x0,y0,ans=;
    for(i = ;i < ;i++)//中国剩余定理
    {
        d=M/w[i];
        exgcd(d,w[i],x0,y0);
        ans=(ans+d*x0*a[i])%M;
    }
    if(ans <= ) ans += M;
    return ans;
}
int main()
{
    init();
    LL g,n;
    while(cin>>n>>g)
    {
        memset(a,,sizeof(a));
        for(int i=;i*i<=n;i++)
        {
            if(n%i==)
            {
                LL tmp=n/i;
                for(int j=;j<;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%MOD,CRT(),MOD)<<endl;
    }
    return ;
}

我就不想多说了=-=

(转载请注明出处：http://www.cnblogs.com/Kalenda/)