杭州网络预赛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;

}

请见:积性函数