题目就不废话了 传送门:http://acm.hdu.edu.cn/showproblem.php?pid=3463
首先我们肯定都能想到去打表
然后我们从3-20000打出2261个素数
对于偶数而言 O(2261*2261)的枚举还是可以接受的 但是对于奇数而言O(2261*2261*2261)就显然要超时了
所以我们换种思路 因为偶数由2个素数构成 而奇数是由3个构成 所以我们可以在构造偶数的同时构造奇数
设result[i] 表示数i的表示方法数
开始先初始化result为0
我们先枚举所有的素数 prime[i] 然后用已经得到的偶数的结果来构造一个奇数(初始的时候还没有得到任何偶数的结果)
然后在使用prime[i]前 我们已经得到的偶数的结果中是不含prime[i]的 因此有result[j+prime[i]]+=result[j](j为偶数)
也就是对j+prime[i]这个奇数,每一次把使用了这个prime[i]的情况加进去(由于j的结果还不含prime[i],所以不会得到重复)
然后我们再使用这个prime[i],从prime[0]-prime[i-1],把使用了这个prime[i]的可能性加到result[prime[t]+prime[i]](t从0到i-1)
这样就是对prime[t]+prime[i]这个偶数,把使用了prime[i]的情况加进去
这样我们可以确保每次构造时不会发生重复 并且时间复杂度只有O(2261*20000)
贴段AC代码 46ms
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <math.h>
using namespace std;
int result[20001];
int prime[2300];
int count;
inline bool isPrime(int n)
{
if(n%2==0) return false;
for(int i=0; i<count; i++)
if(n%prime[i]==0) return false;
return true;
}
int main(void)
{
memset(result,0,sizeof(result));
count = 0;
for(int i=3; i<=19999; i++)
if(isPrime(i)) prime[count++]=i;
for(int i=0; i<count; i++)
{
for(int j=0; j<=20000; j+=2)
{
if(prime[i]+j<=20000){
result[prime[i]+j]+=result[j];
}
else break;
}
for(int j=0;j<i;j++){
if(prime[i]+prime[j]<=20000)
result[prime[i]+prime[j]]++;
else break;
}
}
int m;
while(scanf("%d",&m)!=EOF) printf("%d\n",result[m]);
}