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;
}