【网络流与线性规划24题】【机器人路径规划问题】【IDA*】【题解】

传送门:http://www.wikioi.com/problem/1236/

丧心病狂的一道题,蒟蒻的IDA*只能写到这份上了，dep在24以内都可以承受……估计国外大爹们的O(n^6)的算法还没有我骗分高吧……

国外大爹的研究成果:http://wenku.baidu.com/link?url=Ia3BIZISzgXFNL8v6d2dS9tS5Vrg2heCZkUW3R3HU8kM2X-sI2s73BbsMyzZur1TKqrPC8lRIc8IgKJ0uAVEV-z-GtqK0vP_DWWeh48MOBK

我的Code:

#include<queue>
#include<stack>
#include<cstdio>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Nosolution ((n==19)? 20: 30)
using namespace std;
int w[1001][1001];
void add(int u,int v){
    w[u][++w[u][0]]=v;
    w[v][++w[v][0]]=u;
}
int Brk[1001];
int b[1001];
int lca[1001][1001];
int s,t,n;
int dep;
int now;
int d[1001];
int vis[1001];
int p[1001];
int hash[1001];
void spfa(int s){
    int vis[1001];
    memset(d,0x7f,sizeof(d));
    memset(vis,0,sizeof(vis));
    d[s]=0;
    queue<int>q;
    q.push(s);
    vis[s]=1;
    while(!q.empty()){
        int u=q.front();q.pop();vis[u]=0;
        for(int j=1;j<=w[u][0];j++){
            int v=w[u][j];
            if(d[v]>d[u]+1){
                d[v]=d[u]+1;
                p[v]=u;
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
}
int cost=0;
bool dfs(int deep){
    if(b[t]==2)return true;
    if(deep+lca[now][t]+cost>dep)return false;
    for(int i=1;i<=Brk[0];i++){
        int u=Brk[i];
        for(int j=1;j<=w[u][0];j++){
            int v=w[u][j],O;
            if(!b[v]){
                Brk[i]=v;
                if(u==now)now=v;
                if(hash[u]&&!hash[v]&&b[u]!=2)cost--;
                if(hash[v]&&!hash[u]&&b[u]!=2)cost++;
                vis[v]=u;
                swap(b[v],b[u]);
                if(dfs(deep+1))
                    return true;
                if(hash[u]&&!hash[v]&&b[v]!=2)cost++;
                if(hash[v]&&!hash[u]&&b[v]!=2)cost--;
                if(now==v)now=u;
                vis[v]=O;
                Brk[i]=u;
                swap(b[v],b[u]);
            }
        }
    }
    return false;
}
int main(){
    cin>>n>>s>>t;now=s;Brk[++Brk[0]]=s;
    //G.resize(n+1);
    for(int i=0;i<n;i++){
        cin>>b[i];
        b[i]^=1;
        if(b[i])Brk[++Brk[0]]=i;
        int s,u,v;
        cin>>s;
        for(int j=0;j<s;j++){
            cin>>v;
            if(!lca[i][v]&&lca[v][i]){
                add(i,v);
            }
            lca[i][v]=1;
        }    
    }

    b[s]=2;
//    for(int i=0;i<Brk.size();i++)cout<<Brk[i]<<" ";
//    for(int i=0;i<edges.size();i++){
//        printf("#%d  %d -> %d\n",
//        i,edges[i].u,edges[i].v);
//    }
    for(int i=0;i<n;i++){
        spfa(i);
        for(int j=0;j<n;j++)
            lca[i][j]=d[j];
        if(i==s){
            int u=t;
            while(u!=s){
                hash[u]=1;
                if(b[u]==1)cost++;
                u=p[u];
            }
            hash[s]=1;
        }
    }memset(vis,-1,sizeof(vis));
    for(dep=1;dep<=6;dep++){
        
        if(dfs(0)){
            cout<<dep<<endl;
            return 0;
        }
    
    }
    cout<<Nosolution<<endl;
    return 0;
}