杭州网络预赛1004 Deciphering Password

 

真糟糕，被%耍了一晚多写两个就超时了……真不应该到现在才来把这题干掉

#include <iostream>
using namespace std;

#define MAXN 1001

__int64 a,b;
int prim[MAXN],pcnt;
bool is_prim[MAXN];

void prepare(){
    memset(is_prim,true,sizeof(is_prim));
    pcnt=0;
    __int64 i,j;
    for(i=2;i<MAXN;i++){
        if(is_prim[i]){
            for(j=i;i*j<MAXN;j++)
                is_prim[i*j]=false;
            prim[pcnt++]=i;
        }
    }
}



int main(){
    __int64 i,j,k,ans,tt,T=1;
    prepare();
    while(scanf("%d%d",&a,&b)!=EOF){
        ans=1;
        
        b%=10007;//( j * b + 1) %10007 = ( j % 10007 * b % 10007 % 10007 + 1 ) % 10007
        for(i=0;i<pcnt && prim[i]*prim[i]<=a && a!=1;i++){
            for(j=0;a%prim[i]==0;j++)
                a/=prim[i];
            tt=(j*b+1)*(j*b+2)/2%10007;//立方和公式
            tt*=tt;
            ans=ans*tt%10007;
        }
        if(a!=1){
            tt=(b+1)*(b+2)/2%10007;
            tt*=tt;
            ans=ans*tt%10007;
        }
        printf("Case %I64d: %I64d\n",T++,ans);
    }
    return 0;
}

 

请见:积性函数