http://blog.csdn.net/z309241990/article/details/11917495
题目大意:给一个数n, 最多由三个素数通过加、乘运算得到n 有多少种组合。
解题思路:
首先通过筛选法找出所有的素数,然后依次算出2个素数的和与积的个数,对于每次输入的n,采取如下措施:
1、若n为素数,则答案+1;
2、直接加上n由2个素数相加和相乘的所有情况数、
3、对于3个数,我们只需要枚举第一个素数,因为根据前面我们已经得到所有2个素数相加和相乘的所有的情况数(分开算),然后剩下的就比较简单了。分3种情况
对于加法:
只可能出现a+b+c, a+b+b,a+a+a的情况,由于前面算的时候默认第一个数比第二个数小,所以2个数不存在重复与交换的情况。仔细想想会发现对于a+b+c我们依次枚举前面的a那么会被重复算3次(分别为a+b+c,b+a+c,c+a+b),同样的对于a+b+b的情况我们算了2次(a+b+b,b+a+b)然而对于a+a+a我们只算了一次,所以我的想法就是对于a+b+b的情况我们人为的加一次,而对于a+a+a的情况我们人为的加2次,所以最终的方案数就等于总的计数除以3
对于纯乘法跟纯加法,补加后面除以3.
最后就是对于混合计算~
这个不需要去重,直接枚举素数,加上后面的2个素数乘积的方案数即可
代码实现:
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
using namespace std;
#define N 80100
#define MOD 1000000007
int is[N], prm[N];
int na[N], nb[N];
int res = 0, ans;
int Getprm()
{
int i, j, k = 0;
int s, e = (int)(sqrt(N * 1.0) + 1);
memset(is, 1, sizeof(is));
prm[k++] = 2;
is[0] = is[1] = 0;
for(i = 4; i < N; i += 2) is[i] = 0;
for(i = 3; i < e; i += 2)
if(is[i])
{
prm[k++] = i;
for(s = i *2, j = i * i; j < N; j += s)
is[j] = 0;
}
for(; i < N; i += 2) if(is[i]) prm[k++] = i;
return k;
}
void init()
{
memset(na, 0, sizeof(na));
memset(nb, 0, sizeof(nb));
for(int i = 0; i < ans; i ++)
{
for(int j = i; j < ans; j ++)
{
int temp = prm[i] + prm[j];
if(temp > N) break;
na[temp] ++;
}
}
for(int i = 0 ; i < ans; i ++)
{
for(int j = i; j < ans; j ++)
{
long long temp = (long long)prm[i] * (long long)prm[j];
if(temp > N) break;
nb[temp] ++;
}
}
}
int main()
{
ans = Getprm();
init();
int n;
while(~scanf("%d", &n))
{
res = 0;
if(is[n]) res ++;
res += na[n] + nb[n];
int cn = 0;
for(int i = 0; prm[i] < n; i ++)
{
cn += na[n - prm[i]];
}
for(int i = 0; prm[i] * 2 < n; i ++)
{
if(is[n - prm[i] * 2])
{
cn ++;
if(prm[i] * 3 == n)
cn ++;
}
}
res = (res + cn / 3) % MOD;
cn = 0;
for(int i = 0; prm[i] < n; i ++)
{
if(n % prm[i]) continue;
cn += nb[n / prm[i]];
}
for(int i = 0; prm[i] * prm[i] <= n; i ++)
{
if(n % (prm[i] * prm[i])) continue;
if(is[n/prm[i]/prm[i]])
{
cn ++;
if(prm[i] * prm[i] * prm[i] == n)
cn ++;
}
}
res = (res + cn / 3) % MOD;
for(int i = 0; prm[i] < n; i ++)
{
res = res + nb[n-prm[i]];
}
printf("%d\n", res % MOD);
}
return 0;
}