日期:2014-05-16 浏览次数:20453 次
迟到了一天的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; }