前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >【HDU - 4345 】Permutation(DP)

【HDU - 4345 】Permutation(DP)

作者头像
饶文津
发布2020-06-02 15:08:28
2510
发布2020-06-02 15:08:28
举报
文章被收录于专栏:饶文津的专栏饶文津的专栏

BUPT2017 wintertraining(15) #8F

题意

1到n的排列,经过几次置换(也是一个排列)回到原来的排列,就是循环了。 现在给n(<=1000),求循环周期的所有可能数。

题解

问题等价于几个正整数加起来等于n,求最小公倍数的可能数。 因为1不影响最小公倍数,所以等价于求几个正整数加起来小于等于n,最小公倍数的可能数。

最小公倍数与每个质因子在正整数里最大出现次数有关,所以枚举质因子的幂,进行dp。

dp[i][j]表示前i个质数,和为j时,最小公倍数的可能数。

dp[0][0]=1 转移就是

(dp[i][j]=sum_{k=1}^{p[i]^kle j}{dp[i-1][j-p[i]^k]})

最后把dp[cnt][j]累加起来,答案就是dp[cnt][n]了。 可以写成一维的。

代码

代码语言:javascript
复制
#include <cstdio>
#define N 1001
typedef long long ll;
ll dp[N];
bool com[N];
int cnt;
int prime[N];
int main(){
	for(int i=2;i<N;++i)
	if(!com[i]){
		prime[++cnt]=i;
		for(int j=i+i;j<N;j+=i)
			com[j]=true;
	}
	dp[0]=1;
	for(int i=1;i<=cnt;++i)
		for(int j=N-1;j>=prime[i];--j)
			for(int k=prime[i];k<=j&&k<N;k*=prime[i])
				dp[j]+=dp[j-k];
	for(int i=1;i<N;++i)
		dp[i]+=dp[i-1];
	int n;
	while(~scanf("%d",&n))
		printf("%lld\n", dp[n]);
	return 0;
}
本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
原始发表:2017-07-11 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 题意
  • 题解
  • 代码
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档