洛谷 P6280 [USACO20OPEN] Exercise G

题目来源于：洛谷

题目本质：dp，素数筛法，质数

本题与P4161基本一模一样

首先，分析题目发现，某个排列的需要进行恰好 K 步变回原样，这个时候K的值就是这个排序中各个环的长的的最小公倍数(lcm)。然后需要计算所有K的和，就是计算所有符合条件的质数的乘积的和。其次根据唯一分解，一些数的最小公倍数就是除了1以外其它各个数分解质因数后各个质数的最大次幂相乘。最后，枚举每个质数取几个的情况，设dp(i,j)为用了前i 个质数组成的环的长度的和为j 的时候的所有满足条件的K的值。

完整代码如下：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e4+5;
int p[N],c[N];		//p[i]表示第i个质数，c[i]表示最多能用几个
ll dp[N][N];		//dp[i][j]前i个物品拼成面值为j 
int isp(int x){
	if(x==1){
		return 0;
	}
	for(int i=2;i*i<=x;i++){
		if(x%i==0){
			return 0;
		}
	}
	return 1; 
}
int main(){
	int n,m;
	int cnt=1;
	cin>>n>>m;
	for(int i=2;i<=n;i++){
		if(isp(i)){
			p[cnt]=i;
			c[cnt]=0;
			int x=i;
			while(x<=n){
				c[cnt]++;
				x*=i;	
			}
			cnt++; 
		}
	}
	memset(dp,0,sizeof(dp));
	dp[0][0]=1;		//将数组初始化 
	for(int i=1;i<cnt;i++){
		for(int k=0;k<=n;k++){
			dp[i][k]=dp[i-1][k];
		}
		for(int k=p[i];k<=n;k++){
			int x=p[i];
			for(int j=1;j<=c[i]&&x<=k;j++){
				dp[i][k]=(dp[i][k]+(dp[i-1][k-x]*x)%m)%m;		//第i个质用j次,x=p[i]^j 
				x*=p[i];				
			}
		}	
	}
	ll ans=0;
	for(int i=2;i<=n;i++){
		ans=(ans+dp[cnt-1][i])%m;
	}
	printf("%lld\n",ans+1);
	return 0;
}