POJ 3590 The shuffle Problem [置换群 DP]

传送门

---

$1A$太爽了

从此$Candy?$完全理解了这种$DP$做法

和bzoj1025类似，不过是求最大的公倍数，并输出一个字典序最小的方案

依旧枚举质因子和次数，不足的划分成1

输出方案从循环长度小的到大的输出

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=;
typedef long long ll;
inline int read(){
    char c=getchar();int x=,f=;
    while(c<''||c>''){if(c=='-')f=-; c=getchar();}
    while(c>=''&&c<=''){x=x*+c-''; c=getchar();}
    return x*f;
}
int n,m;
int p[N];
bool notp[N];
void sieve(int n){
    for(int i=;i<=n;i++){
        if(!notp[i]) p[++m]=i;
        for(int j=;j<=m&&i*p[j]<=n;j++){
            notp[i*p[j]]=;
            if(i%p[j]==) break;
        }
    }
}
int f[N][N],g[N][N];
int a[N],tot;
void getPath(int i,int j){
    if(i==) return;
    if(g[i][j]!=) a[++tot]=g[i][j];
    getPath(i-,j-g[i][j]);
}
void dp(){
    f[][]=;
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++){
            f[i][j]=f[i-][j];
            g[i][j]=;
            for(int k=p[i];k<=j;k*=p[i])
                if(f[i-][j-k]*k>f[i][j]){
                    f[i][j]=f[i-][j-k]*k;
                    g[i][j]=k;
                }
        }
    int ans=,now=;
    for(int i=;i<=n;i++) if(f[m][i]>f[m][ans]) ans=i;
    printf("%d ",f[m][ans]);
    for(;now<=n-ans;now++) printf("%d ",now);
    //puts("totot");
    tot=;
    getPath(m,ans);
    sort(a+,a++tot);
    //for(int i=1;i<=tot;i++) printf("a %d %d\n",i,a[i]);
    for(int i=;i<=tot;i++){
        int _=now;
        for(int j=;j<a[i];j++) printf("%d ",++now);
        printf("%d ",_);
        now++;
    }
    puts("");
}
int main(){
    freopen("in","r",stdin);
    sieve();
    int T=read();
    while(T--){
        n=read();
        dp();
    }
}