poj 3590 The shuffle Problem——DP+置换

题目：http://poj.org/problem?id=3590

bzoj 1025 的弱化版。大概一样的 dp 。

输出方案的时候小的环靠前。不用担心 dp 时用 > 还是 >= 来转移，因为答案的那个 lcm 有唯一分解。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=,M=;
int T,n,dp[M][N],dy[M][N],pri[M],cnt,sta[M],top,ans,sm;
bool vis[N];
int Mx(int a,int b){return a>b?a:b;}
void init()
{
  for(int i=;i<=;i++)
    {
      if(!vis[i])pri[++cnt]=i;
      for(int j=;j<=cnt&&i*pri[j]<=;j++)
    {
      vis[i*pri[j]]=;
      if(i%pri[j]==)break;
    }
    }
}
void print()
{
  printf("%d ",ans);
  sort(sta+,sta+top+);
  for(int i=;i<=n-sm;i++)printf("%d ",i);
  for(int i=,cr=n-sm+;i<=top;i++)
    {
      for(int j=cr+,k=;k<sta[i];j++,k++)
    printf("%d ",j);
      printf("%d ",cr); cr+=sta[i];
    }
  puts("");
}
int main()
{
  init();
  scanf("%d",&T);
  while(T--)
    {
      scanf("%d",&n);
      memset(dp,,sizeof dp); memset(dy,,sizeof dy);
      dp[][]=; int i;
      for(i=;i<=cnt&&pri[i]<=n;i++)
    {
      for(int j=;j<=n;j++)dp[i][j]=dp[i-][j];
      for(int k=,p=pri[i];k<=&&p<=n;k++,p*=pri[i])//&&p<=n!! or exceed int
        for(int j=p,d;j<=n;j++)
          if((d=dp[i-][j-p]*p)>dp[i][j])//>? >=? no influence for 'only'
        {
          dp[i][j]=d; dy[i][j]=p;
        }
    }
      i--; top=; ans=;
      for(int j=;j<=n;j++)
    if(dp[i][j]>ans)ans=dp[i][j],sm=j;
      for(int t=i,cr=sm;t;t--)
    {
      if(!dy[t][cr])continue;
      sta[++top]=dy[t][cr];
      cr-=dy[t][cr];
    }
      print();
    }
  return ;
}