阶乘质因子分解——lightoj1035

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define maxn 200

int primes[maxn],m;
bool vis[maxn];
void init(){
    for(int i=;i<maxn;i++){
        if(!vis[i])
            primes[++m]=i;
        for(int j=;j<=m;j++){
            if(i*primes[j]>=maxn)break;
            vis[i*primes[j]]=;
            if(i%primes[j]==)break;
        }
    }
} 

int c[maxn];
void divide(int n){
    memset(c,,sizeof c);
    for(int i=;i<=m && primes[i]<=n;i++){
        int p=primes[i];
        for(;p<=n;p*=primes[i])
            c[i]+=n/p;
    }
}

int main(){
    int t;cin>>t;init();
    for(int tt=;tt<=t;tt++){
        int n;cin>>n;
        divide(n);
        printf("Case %d: %d = %d (%d)",tt,n,primes[],c[]);
        for(int i=;i<=m;i++)
            if(c[i])
                printf(" * %d (%d)",primes[i],c[i]);
        puts("");
    }
}