BZOJ 1025: [SCOI2009]游戏|动态规划
2016-01-13 15:12
120 查看
这个题需要转化一下
原题中的转化必须要构成一些环
要求的就是环的大小的最小公倍数的个数
满足环的大小的和<=n
令F[i][j]表示前i个质数 和为j的方案数
然后就相当于背包问题了 令prime[i]^k(k=1,2,3......)做物品 往里装
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<algorithm>
#include<iostream>
#define T 1002
#define lowbit(x) (x&(-x))
using namespace std;
int a[T],prime[T],n,tot;
long long f[T][T],ans;
void get_prime()
{
for(int i=2;i<=n;i++)
if(!a[i])
{
prime[++tot]=i;
for(int j=i*2;j<=n;j+=i)
a[j]=1;
}
}
int main()
{
scanf("%d",&n);
get_prime();
f[0][0]=1;
for(int i=1;i<=tot;i++)
{
for(int j=0;j<=n;j++) f[i][j]=f[i-1][j];
for(int j=0;j<=n;j++)
for(int k=prime[i];k<=j;k*=prime[i])
f[i][j]+=f[i-1][j-k];
}
for(int i=0;i<=n;i++)ans+=f[tot][i];
cout<<ans;
return 0;
}
原题中的转化必须要构成一些环
要求的就是环的大小的最小公倍数的个数
满足环的大小的和<=n
令F[i][j]表示前i个质数 和为j的方案数
然后就相当于背包问题了 令prime[i]^k(k=1,2,3......)做物品 往里装
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<algorithm>
#include<iostream>
#define T 1002
#define lowbit(x) (x&(-x))
using namespace std;
int a[T],prime[T],n,tot;
long long f[T][T],ans;
void get_prime()
{
for(int i=2;i<=n;i++)
if(!a[i])
{
prime[++tot]=i;
for(int j=i*2;j<=n;j+=i)
a[j]=1;
}
}
int main()
{
scanf("%d",&n);
get_prime();
f[0][0]=1;
for(int i=1;i<=tot;i++)
{
for(int j=0;j<=n;j++) f[i][j]=f[i-1][j];
for(int j=0;j<=n;j++)
for(int k=prime[i];k<=j;k*=prime[i])
f[i][j]+=f[i-1][j-k];
}
for(int i=0;i<=n;i++)ans+=f[tot][i];
cout<<ans;
return 0;
}
相关文章推荐
- C++动态规划之最长公子序列实例
- C++动态规划之背包问题解决方法
- C#使用动态规划解决0-1背包问题实例分析
- 动态规划
- C++ 动态规划
- DP(动态规划) 解游轮费用问题
- 动态规划的用法——01背包问题
- 动态规划的用法——01背包问题
- 《收集苹果》 动态规划入门
- 《DNA比对》蓝桥杯复赛试题
- 《背包问题》 动态规划
- 初学ACM - 半数集(Half Set)问题 NOJ 1010 / FOJ 1207
- 关于爬楼梯的动态规划算法
- 动态规划 --- hdu 1003 **
- DP问题各种模型的状态转移方程
- 0-1背包解题过程
- 背包问题
- USACO 3.2.2:Stringsobits
- 字符串编辑距离
- HDU ACM Step 2.2.2 Joseph(约瑟夫环问题)