【团体程序设计天梯赛-练习集L3-015】球队“食物链”【状压dp+路径记录】

题意：

    找一个环，满足相邻的前一个胜过后一个，输出答案序列，如果有多个答案，输出字典序最小的方案

题解：

    首先，如果有答案的话，方案中1肯定是在第一个(因为答案是个环

    最多有20个点，状压dp很好想到，这是最简单的了差不多，dp[s][i]代表以i为结尾的状态s能否被达到，把dp[1][1]初始化为1

    难点是路径的记录，我是这样处理的，在记录正向边的时候，用pre[i]记录那些点能够到达i的点(也就是方向边)，然后从后面的状态开始，把满足条件的dp[1<<n-1][i]的状态mark为1,然后枚举pre[i]中的每一个值j,如果dp[((1<<n)-1)^(1<<(i-1)][j]，那把这个状态的mark也标记为1，这样从前往后标记，这样所有答案的路径都会被mark为1，然后就从1号节点开始遍历一遍就好了

#include<cstdio>
#include<algorithm>
using namespace std;
#define PB push_back
const int N=25;
const int M=2e6+10;
int f[M];
int n;
char mp[N][N];
vector<int>G[N],pre[N];
int dp[M][N];
bool vis[N][N],mark[M][N],ok[N][N][N];
int main()
{
    for(int i=0;i<M;i++)
        f[i]=f[i/2]+(i&1);
    while(~scanf("%d",&n)){
        for(int i=1;i<=n;i++)scanf("%s",mp[i]+1);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(mp[i][j]=='W'){
                    if(!vis[i][j]){
                        G[i].PB(j);
                        vis[i][j]=1;
                        pre[j].PB(i);
                    }
                }
                else if(mp[i][j]=='L'){
                    if(!vis[j][i]){
                        G[j].PB(i);
                        vis[j][i]=1;
                        pre[i].PB(j);
                    }
                }
            }
        }
        dp[1][1]=1;
        int mx=1<<n;
        for(int s=1;s<mx;s++){
            for(int i=1;i<=n;i++)if(dp[s][i])
                for(int j=0;j<G[i].size();j++){
                    int v = G[i][j];
                    if(s&(1<<(v-1)))continue;
                    dp[s|(1<<(v-1))][v]=1;
                }
        }
        bool flag=0;
        for(int i=1;i<=n;i++){
            if((dp[mx-1][i])&&(find(G[i].begin(),G[i].end(),1)!=G[i].end()))mark[mx-1][i]=1,flag=1;
        }
        if(!flag){puts("No Solution");continue;}
        for(int i=mx-1;i>1;i--){
            for(int j=1;j<=n;j++)if(mark[i][j]){
                int d=i^(1<<(j-1));
                for(int k=0;k<pre[j].size();k++)if(dp[d][pre[j][k]])
                    mark[d][pre[j][k]]=1;
            }
        }
        printf("1");
        int now=1,s=1;
        for(int i=2;i<=n;i++){
            sort(G[now].begin(),G[now].end());
            for(int j=0;j<G[now].size();j++){
                int v=G[now][j];
                if(s&(1<<(v-1)))continue;
                if(mark[s^(1<<(v-1))][v]){
                    s|=(1<<(v-1));
                    now=v;
                    printf(" %d",v);
                    break;
                }
            }
        }
        puts("");
    }
    return 0;
}