2013 ACM/ICPC 长沙网络赛-Goldbach

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