poj 2191 Mersenne Composite Numbers

#include<iostream>
using namespace std;
bool is_prime(int a)
{
    if(a==1)
        return false;
    for(int i=2;i*i<=a;++i)
        if(a%i==0)
            return false;
    return true;
}
long long num[100]={1},ans[10];
int main()
{
    int n;
    cin>>n;
    for(int i=1;i<64;++i)        
        num[i]=num[i-1]*2;    //无法使用pow(2.0,i)或1<<i
    for(int i=1;i<n;++i)
    {
        int ok=0;
        if(is_prime(i))    //根据题意,i必须是素数
        {
            long long s=num[i]-1;        
            int rear=0;
            for(long long p=2;p*p<=s;++p)    // 筛选s的素因子
            {
                while(s%p==0)        
                {
                    ans[rear++]=p;
                    s=s/p;
                    ok=1;
                }
            }
            if(ok)
            {
                ans[rear]=s;    
                cout<<ans[0];
                for(int j=1;j<=rear;++j)
                    cout<<" * "<<ans[j];
                cout<<" = "<<num[i]-1<<" = ( 2 ^ "<<i<<" ) - 1\n";
            }    
        }
    }
    return 0;
}