求10^8内的半素数

来源：互联网发布：冰川网络股票怎么样编辑：程序博客网时间：2024/05/20 15:10

每个合数都可以写成几个质数相乘的形式。例如
15=3*5; 9=3*3;可以写成两个质数相乘。
12=2*2*3可以写成三个质数相乘。
1~30之间只能写成两个质数相乘的数有4, 6, 9, 10, 14, 15, 21, 22, 25, 26
请问1~99999999之间有多少个合数只能被写成两个质数相乘的形式。《输出17427258》

VS2010 c++控制台程序

我chuo，内存需求太大，另外算乘法是有负数或溢位干扰，不得不用long long

前几次算得why 34853287 44832018，怎么差别那么大的，谷歌答案http://www.r-bloggers.com/project-euler-%E2%80%94-problem-187/，

根据答案修代码

// ceshiyongde.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include <iostream>
using namespace std;

int _tmain(int argc, _TCHAR* argv[])
{ /*
每个合数都可以写成几个质数相乘的形式。例如
15=3*5; 9=3*3;可以写成两个质数相乘。
12=2*2*3可以写成三个质数相乘。
1~30之间只能写成两个质数相乘的数有4, 6, 9, 10, 14, 15, 21, 22, 25, 26
请问1~99999999之间有多少个合数只能被写成两个质数相乘的形式。<17427258> */

long i,j,k,cnt=0,sum,sumSuShuShuHe,temp,count=0,t;
bool *pBoolSuShu=new bool[99999999];
long long *parr=new long long[5761456];
long long chenji;

pBoolSuShu[0]=false;
pBoolSuShu[1]=false;

for(i=2; i<99999999; i++)
{
pBoolSuShu[i]=true;
}

for(i=2; i<99999999; i++)
{
if(pBoolSuShu[i])
{
for(j=i+i; j<99999999; j+=i)
{
pBoolSuShu[j]=false;
}
}
}

for(i=2; i<99999999; i++)
{
if(pBoolSuShu[i])
{
parr[cnt]=i;
cnt++;
}
}

for(i=0; i<cnt; i++)
{
for(j=i; j<cnt; j++)
{
chenji=parr[i] * parr[j];

if( chenji<= 99999999)
count++;
else
break;
}
}

cout<<count<<endl;
delete []pBoolSuShu;
delete []parr;

getchar();
return 0;
}