子集个数 题解

/* 
生成函数, f(x)=(1+x^2)(1+x^3)...(1+x^9973) 
暴力展开生成式,取所有质数系数,和即为答案 
模数为1e9+7时，答案为243425655
*/
#include<bits/stdc++.h>
#define N 6000005
#define Mod 1000000007
#define LL long long
using namespace std;
int num[N], prime[N], c[N], d[N];  // c[i]存放x^i的系数 , d[i]为临时数组 
int pn; 
void getPrime(int n)	
{
	pn = 0;
    memset(num, -1, sizeof(num));
    for(int i = 2; i <= n;i++)
	{
        if(num[i]) prime[++pn] = i;
        for(int j = 1;j <= pn && 1LL*i*prime[j] <= n;j++)
		{
            num[i * prime[j]] = 0;
            if(i % prime[j] == 0) break;
        }
    }
}
void calculateCoefficient()
{
	c[0] = 1;
	int maxSum = 0;
	for(int i = 1; i <= pn; i++)
	{
		maxSum += prime[i];			// maxSum为当前展开到的最大系数 
		for(int j = prime[i]; j <= maxSum; j++)
			d[j] = c[j - prime[i]];
		for(int j =