hdu5823(反演dp)

听说3^n也能水过去。。

其实应该是个经典题，求图染色这个np问题。

把问题拆成独立集来进行dp可以在3^n之内水过去。

拆成独立集的时候就发现，等价与一个经典的反演dp问题

然后复杂度就变成了 n*n*2^n

另外，偷到一套头文件宏定义。

#include <math.h>
#include <time.h>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <set>
#include <map>
#include <string>
#include <stack>
#include <queue>
#include <vector>
#include <bitset>
#include <iostream>
#include <algorithm>
#define pb push_back
#define fi first
#define se second
#define icc(x) (1<<(x))
#define lcc(x) (1ll<<(x))
#define lowbit(x) (x&-x)
#define debug(x) cout<<#x<<"="<<x<<endl
#define rep(i,s,t) for(int i=s;i<t;++i)
#define per(i,s,t) for(int i=t-1;i>=s;--i)
#define mset(g, x) memset(g, x, sizeof(g))
using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef unsigned int ui;
typedef double db;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef vector<int> veci;
const int mod=(int)1e9+,inf=0x3fffffff,rx[]={-,,,},ry[]={,,,-};
const ll INF=1ll<<;
const db pi=acos(-),eps=1e-;

template<class T> void rd(T &res){
    res = ; int ch,sign=;
    while( (ch=getchar())!='-' && !(ch>=''&&ch<=''));
    if(ch == '-') sign = ; else res = ch-'';
    while((ch=getchar())>=''&&ch<='') res = (res<<)+(res<<)+ch-'';
    res = sign?-res:res;
}
template<class T>void rec_pt(T x){
    if(!x)return;
    rec_pt(x/);
    putchar(x%^);
}
template<class T>void pt(T x){
    if(x<) putchar('-'),x=-x;
    if(!x)putchar('');
    else rec_pt(x);
}
template<class T>inline void ptn(T x){ pt(x),putchar('\n'); }
template<class T>inline void Max(T &a,T b){ if(b>a)a=b; }
template<class T>inline void Min(T &a,T b){ if(b<a)a=b; }
template<class T>inline T mgcd(T b,T d){ return b?mgcd(d%b,b):d; }//gcd模板，传入的参数必须是同一类型
//-------------------------------主代码--------------------------------------//

int mat[];
int mark[<<];//记录是否为独立集
int dp[][<<];
int A[<<],B[<<];
int n;
void gao(int _A[],int _B[],int C[])
{
    rep(i, , icc(n)){
        A[i] = _A[i];
        B[i] = _B[i];
    }
    //传说中的or卷积
    rep(i, , n){
        rep(j, , icc(n)){
            if( (icc(i)&j)){
                A[j] += A[j^icc(i)];
                B[j] += B[j^icc(i)];
            }
        }
    }
    rep(i, , icc(n)) C[i] = A[i]*B[i];
    //然后是逆着卷机
    rep(i, , n){
        rep(j, , icc(n)){
            if( (icc(i)&j)){
                C[j] -= C[j^icc(i)];
            }
        }
    }

}
int main()
{
    int T;
    rd(T);
    while(T--)
    {
        rd(n);
        rep(i, , n){
            mat[i] = ;
            rep(j, , n){
                char tmp;
                cin>>tmp;
                if(tmp == '')
                    mat[i] |= icc(j);
            }
        }
        //mset(dp, 0);
        //然后求出独立集
        mark[] = ;
        rep(i, , n){
            rep(j, , icc(i)){
                dp[][icc(i)|j] = mark[icc(i)|j] = mark[j]==?(mat[i]&j)==:;
            }
        }

        //然后开始dp

        rep(i, , n){
            if(dp[i][icc(n)-]) break;//这步优化很重要啊
            gao(dp[i],mark,dp[i+]);
            rep(j, , icc(n)) {
                dp[i+][j] = (dp[i+][j]!=);
                //pt(dp[i+1][j]); putchar(' ');
            }
            //puts("");
        }

        ui ans = ;
        ui tmp = ;
        rep(i, , icc(n))
        {
            ui mj = ;
            tmp *= ;
            rep(j, , n+){
                if(dp[j][i]!=){
                    mj = j;
                    break;
                }
            }
            ans += mj*tmp;
        }
        ptn(ans);
    }
    return ;
}

/*
//-----------Test Case------------//

*/