POJ 2992 Divisors (因子个数)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>

using namespace std;
typedef long long LL;

vector <int > vt;
int a[440][440];

int prime[500],nprime;
bool isprime[500];

void doprime()
{
    nprime=0;
    LL i,j;
    for(i=1;i<500;i+=2)
        isprime[i]=true;
    isprime[1]=false;
    isprime[2]=true;
    for(i=2;i<500;i++)
    {
        if(isprime[i])
        {
            prime[nprime++]=i;
            for(j=i*i;j<500;j+=i)
                isprime[j]=false;
        }
    }
}

void divd(int x)
{
    vt.clear();
    int i,j,t=x;
    for(i=0;i<nprime && prime[i]*prime[i]<=x;i++)
    {
        while(t%prime[i]==0)
        {
            vt.push_back(prime[i]);
            t/=prime[i];
        }
        if(t==1)
            break;
    }
    if(t>1)
        vt.push_back(t);
}

void work()
{
    int i,j;
    for(i=2;i<=431;i++)
    {
        memcpy(a[i],a[i-1],sizeof(a[i]));
        divd(i);
        for(j=0;j<vt.size();j++)
            a[i][vt[j]]++;
    }
}

int main()
{
    doprime();
    work();
    int n,k,i,j;
    LL ans;
    int use[500];
    while(scanf("%d%d",&n,&k)!=EOF)
    {
        memset(use,0,sizeof(use));
        for(i=2;i<=431;i++)
            use[i]+=a
[i]-a[k][i]-a[n-k][i];
        ans=1;
        for(i=2;i<=431;i++)
            ans*=(use[i]+1ll);
        printf("%lld\n",ans);
    }
    return 0;
}