ZOJ Goldbach 2013年长沙赛区网络赛

迟到了一天的AC。。。。

思路：

先把单个素数 或着 两个素数能组成的情况预处理一下，然后对于给出的 n，拿第三个素数去和两个素数的情况匹配，最后要注意去重。

详情见代码。

因为手残少敲了一个 else ，Debug了一晚上。。。

 

#include <iostream>

#include <cstring>

#include <cstdlib>

#include <cstdio>





using namespace std;





const long long int INF = 1000000007;





bool MarkVisit[80010] = {false};

long long int pri[80010];

long long int  s[80010] = {0};

long long int sm[80010] = {0};//a*a

long long int sa[80010] = {0};//a+a

long long int dm[80010] = {0};//a*b

long long int da[80010] = {0};//a+b





int main()

{

    int i,j;





    int top = 0;





    for(i = 2; i <= 80000; ++i)//挫到荼靡的素数筛。。。

    {

        if(MarkVisit[i] == false)

        {

            pri[top++] = i;





            for(j = i+i; j <= 80000; j += i)

            {

                MarkVisit[j] = true;

            }

        }

    }









    for(i = 0; i < top; ++i)//预处理

    {

        s[pri[i]]++;//一个素数的情况





        if(pri[i]*2 <= 80000)

        {

            sa[pri[i]*2]++;//两个相同的素数相加的情况

        }

        if(pri[i]*pri[i] <= 80000)

        {

            sm[pri[i]*pri[i]]++;//两个相同的素数相乘的情况

        }





        for(j = i+1; j < top; ++j)

        {

            if(pri[i] + pri[j] <= 80000)

            {

                da[pri[i] + pri[j]]++;两个不相同的素数相加的情况

            }





            if(pri[i]*pri[j] <= 80000)

            {

                dm[pri[i]*pri[j]]++;//两个不相同的素数相乘的情况

            }

        }

    }





    int n;





    long long int tempsum, sum;





    while(scanf("%d",&n) != EOF)

    {

        sum = s[n] + dm[n] + da[n] + sa[n] + sm[n];//先把一个或两个素数的情况统计出来





        for(i = 0,tempsum = 0; i < top; ++i)

        {

            if(n%pri[i] == 0)

            {

                if(n == pri[i]*pri[i]*pri[i])//如果三个数相同 则此种情况重复计算一次 *3 补成三次

                {

                    tempsum += sm[n/pri[i]]*3 + dm[n/pri[i]];

                }

                else tempsum += sm[n/pri[i]]*2 + dm[n/pri[i]];//如果只有来两个相同 则此种情况重复计算两次。

                                                              //因为另一种单独计算 所以只需对其中一次*2 补成三次





            }

            if(n >= pri[i])

            {

                sum += sm[n-pri[i]];//+ && *都出现的不会重复

                sum += dm[n-pri[i]];

                sum %= INF;





                if(n == pri[i]*3)//此种情况 和 * 形同  不再赘述

                {

                    tempsum += da[n-pri[i]] + sa[n-pri[i]]*3;





                }





                else tempsum += da[n-pri[i]] + sa[n-pri[i]]*2;

            //就是因为少了这个else...险些把这一晚上都废了。。。

            }

            else break;

        }

        sum %= INF;

        sum += tempsum/3;//因为每一种都出现了3次

        sum %= INF;





        printf("%lld\n",sum);

    }

    return 0;

}