状压DP假装学习

• 直接学状压题型 = =
• 哇，好难
• 哇，好像看懂了吧 = =
• 打不出来 = =
• 忘了 = =
• yi 这位运算什么意思？

先学一波位运算吧

1. 1 << n   即1*2^n，一般用来表示状态, Eg: n = 2,1 << n = 4 即可表示二进
    制四位的状态 Eg: 1111
2. i & (i<<1) 
         Eg: i = 2;
         i << 1 = 8;
         2 & 8 =
             1000
           & 0001
        ______________
             0000
/*code*/
        tot = 0;    //tot -> 合法的状态标号 用can储存状态 Eg:can[3] = 5 = (101)_2
        for(int i = 0;i < (1<<n);i++){
            if((i & (i<<1)) == 0){         //判断有没有相邻的
                can[tot++] = i;            //初始化状态
            }
        }
3. cur[i] = (cur[i] | (1<<j))
    Eg:cur[i] = 5 // (101)_2
       j = 1;
       1 << j = 2; //(010)_2
       cur[i] = 
          101
       |  010
       _________
          111
       即在第二个位置插入 '1'
       用来修改地图 

/*code*/
        for(int i = 1;i <= m;i++){
            for(int j = 0;j < n;j++){
                int num;
                scanf("%d",&num);
                if(num == 0) cur[i] = (cur[i] | (1<<j));
            }
        }

前方高能膜法门:https://blog.csdn.net/accry/article/details/6607703

旅行商问题(TSP)

int dp[1 << MAX_N][MAX_N];  //DP记忆数组
                            //用二进制的每一位代表该点是否访问过
                            //0011代表已走过第一个和第二个

void solve(){
    //初始化
    for(int S = 0;S < 1 << n;S++){
        fill(do[S],dp[S]+n,INF);
    }
    //初始化从零点访问
    dp[(1 << n) - 1][0] = 0;
    //核心
    for(int S = (1 << n)-2; S >= 0;S--){    //从全访问 -> 0访问
        for(int v = 0;v < n;v++){           //从0点到最后一点
            for(int u = 0;u < n;u++){       
                if(!(S >> u & 1)){          //判断每位(点)是否访问过，若没访问过则递推式成立
                    dp[S][v] = min(dp[S][v],dp[S | 1 << u][u] + d[v][u]);   //不访问d[v][u]和访问d[v][u[]]
                }
            }
        }
    }
    cout<<dp[0][0]
}

Corn Fields

• 这是第二天了啊 = =
• 靠着前一天题解的记忆打了好久终于成功了
• 哼哼 有牌面

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>

#define MAX_N 13
#define mod 100000000
#define mm(a,b) memset(a,b,sizeof(a))

int dp[MAX_N][1 << MAX_N];      //记忆数组
int map[MAX_N];                 //field地图
int can[1 << MAX_N];            //第几种状态的二进制表示
int cnt;                        //第几种状态
int N,M;                        //大小
int ans;                        //答案

void TestMap();             //输出地图
void TestInitCan();         //输出初始化Can
void TestInitDpOne();       //输出初始化第一行Dp

using namespace std;
int main(){
    while(~scanf("%d %d",&N,&M)){
        cnt = 0;
        for(int i = 0;i < (1<<M);i++){
            if((i & (i<<1)) == 0){          //判断有没有相邻的
                can[cnt] = i;
                cnt++;
            }
        }
        //初始化
        mm(map, 0);
        mm(dp, 0);
        ans = 0;
        //读取地图,0 -> 肥沃,1 -> 贫瘠
        for(int i = 1;i <= N;i++){
            for(int j = 0;j < M;j++){
                int tmp;
                scanf("%d",&tmp);
                if(tmp == 0){
                    map[i] = (map[i] |(1<<j));
                }
            }
        }
        //对第一行进行初始化
        for(int i = 0;i < cnt;i++){
            if((map[1] & can[i]) == 0){
                dp[1][can[i]] = 1;
            }
        }

        for(int i = 1;i < N;i++){                       //枚举每一行
            for(int j = 0;j < cnt;j++){                 //枚举每一行的所有可行状态
                if((can[j] & map[i]) == 0){             //满足条件
                    for(int k = 0;k < cnt;k++){         //枚举下一行的所有可行状态
                        if(((can[k] & map[i+1]) == 0) && ((can[k] & can[j]) == 0)){
                            dp[i+1][can[k]] = dp[i+1][can[k]] + dp[i][can[j]];
                        }
                    }
                }
            }
        }

        for(int i = 0;i < cnt;i++){  //各个状态ans之和
            ans += dp[N][can[i]];
            ans = ans % mod;
        }
        cout<<ans<<endl;
        /* 发现结果错了，输出了三个中间值，发现是输入地图的时候 0 和 1 取反了*/
        //TestMap();
        //TestInitCan();
        //TestInitDpOne();
    }
    return 0;
}

void TestMap(){
    cout<<"Map:"<<endl;
    for(int i = 1;i <= N;i++){
        cout<<map[i]<<endl;
    }
    return;
}

void TestInitCan(){
    cout<<"Can:"<<endl;
    for(int i = 0;i < cnt;i++){
        cout<<can[i]<<endl;
    }
    return;
}

void TestInitDpOne(){
    cout<<"DpOne:"<<endl;
    for(int i = 0;i < cnt;i++){
        cout<<"map[1] = "<<map[1]<<"   can[i] = "<<can[i]<<endl;
        cout<<"ans: "<<dp[1][can[i]]<<endl;
    }
    return;
}