POJ 3421 X-factor Chains

线型素数筛+质因素分解+组合数。

AC后发现这样做效率有点低。。766ms。

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<time.h>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;

typedef long long LL;
int tol;
LL factor[];
const int N = ;
LL prime[N] = {},num_prime = ;
int isNotPrime[N] = {, };

void findfac(long long n)
{
    for(int i=; i<num_prime; i++)
    {
        if(n%prime[i]!=) continue;
        while()
        {
            if(n==||n%prime[i]!=) break;
            factor[tol++]=prime[i];
            n=n/prime[i];
        }
        if(n==) break;
    }
}

int num[],tot;
long long c[][];

void init()
{
    c[][]=;
    for(int i=; i<=; i++) c[i][]=;
    for(int i=; i<=; i++)
    {
        for(int j=; j<=; j++)
        {
            c[i][j]=c[i-][j-]+c[i-][j];
        }
    }

    for(long i =  ; i < N ; i ++)
    {
        if(! isNotPrime[i])
            prime[num_prime ++]=i;
        for(long j =  ; j < num_prime && i * prime[j] <  N ; j ++)
        {
            isNotPrime[i * prime[j]] = ;
            if( !(i % prime[j] ) )
                break;
        }
    }
}

int main()
{
    init();
    long long n;

    while(~scanf("%lld",&n))
    {
        if(n==) {
            printf("%d %lld\n",,);
            continue;
        }
        tol=;
        findfac(n);
        sort(factor,factor+tol);

        tot=,num[tot]=;
        for(int i=;i<tol;i++)
        {
            if(factor[i]==factor[i-]) num[tot]++;
            else
            {
                tot++;
                num[tot]=;
            }
        }
        tot++;
        int sum=tol;
        int ans1=tol;
        long long ans2=;

        for(int i=;i<tot;i++)
        {
            ans2=ans2*c[sum][num[i]];
            sum=sum-num[i];
        }
        printf("%d %lld\n",ans1,ans2);
    }

    return ;
}