poj3311

poj3311
我写的第一篇状压dp，调试了好久23333，其实状压就是二进制结合位运算，把整体的状态用二进制来表示，一般是0/1，用没用过或者走没走没走过的问题，其实这种思想在搜索中也有涉及过，与dp结合灵活度较高。状态占1维，还有其他维。
这个题就是f[i][j]表示，i状态下，当前在j位置的最短距离，转移过程与floyed相似。第一次尝试就掉坑（痛苦的体验啊

#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
#include<cmath>
#include<ctime>
#include<set>
#include<map>
#include<stack>
#include<cstring>
#define inf 2147483647
#define ls rt<<1
#define rs rt<<1|1
#define lson ls,nl,mid,l,r
#define rson rs,mid+1,nr,l,r
#define N 100010
#define For(i,a,b) for(long long i=a;i<=b;i++)
#define p(a) putchar(a)
#define g() getchar()

using namespace std;
long long n;
long long dis[20][20],m[20][20];
long long f[1<<12][20];
void in(long long &x){
    long long y=1;
    char c=g();x=0;
    while(c<'0'||c>'9'){
        if(c=='-')y=-1;
        c=g();
    }
    while(c<='9'&&c>='0'){
        x=(x<<1)+(x<<3)+c-'0';c=g();
    }
    x*=y;
}
void o(long long x){
    if(x<0){
        p('-');
        x=-x;
    }
    if(x>9)o(x/10);
    p(x%10+'0');
}
int main(){
    while(cin>>n&&n){

    For(i,0,n)
      For(j,0,n){
          in(m[i][j]);
          dis[i][j]=m[i][j];
      }
          
      
    For(i,0,n)
        For(j,0,n)
            For(k,0,n)
                if(dis[i][j]>dis[i][k]+m[k][j])
                    dis[i][j]=dis[i][k]+m[k][j];

    memset(f,-1,sizeof(f));
    f[1][0]=0;

    For(i,1,(1<<(n+1))-1){
        i=i|1;
        For(j,0,n){
            if(f[i][j]!=-1)
            For(k,0,n)
                if((j!=k)&&(f[(1<<k)|i][k]==-1||f[(1<<k)|i][k]>f[i][j]+dis[j][k]))
                    f[(1<<k)|i][k]=f[i][j]+dis[j][k];
                
            
        }

    }
    o(f[(1<<(n+1))-1][0]);
    p('
');
}
    return 0;
}