NBUT 1580 素数分解dp

这道题目以前就看到过，现在想起来了，补一下题，动态规划还真不好想啊，转移方程如果能想起来的话，基本上就

能A了。。

思路：

先对素数打表，然后设一个数组d[i][j+1]用来表示前j个素数和为i的分解个数，然后我们对第j个素数进行分析，

有两种情况，用或者不用，即d[i][j]和d[i-primer[j]][j+1](注意这里可以用多次，所以用了一次后还可以在

前j个素数里面选择),然后把所有的d[0][j+1]初始化为1，即i=0的时候有解。其它初始化为0。

贴代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<vector>
#include<set>
#include<string>
#include<algorithm>
int primer[205];
using namespace std;
int init()
{
	int i,j,cnt=1;
	int visit[205];
	memset(visit,0,sizeof(visit));
	for(i=2; i<=200; i++)
	{	
		if(!visit[i])
			primer[cnt++] = i;
		for(j=i+i; j<=200; j+=i)
		{
			visit[j] = 1;
		}
	}
	return cnt; 
}
int main()
{	
	int i,j,n;
	int d[205][205];
	int m = init();
	memset(d,0,sizeof(d));
	for(j=0; j<m; j++)
			d[0][j+1] = 1;
	for(i=2; i<=200; i++)
		for(j=1; j<m; j++)
			{
				if(i >= primer[j])
					d[i][j+1] += d[i][j]+d[i-primer[j]][j+1];
				else
				{
					d[i][j+1] = d[i][j];
				}
			}
	while(cin >> n)
	{
		cout << d[n][m] << endl;
	} 
	return 0;
}