[DP套DP] HDU5079: Square

题意

给出一个n*n的网格，每个格子一开始都是黑的，且有些格子是损坏的。
现在可以任意选择一些非损坏格子染成白色。定义一个染色方案的优美度为白色最大子正方形的边长。
对于每个i=0~n，求优美度恰为i的染色方案。
n<=10

题解

显然是DP套DP。
首先思考最大子正方形这个全局的东西比较难用一个较简单的状态直接体现出来。
所以就分开来做，考虑容斥，如果我们求出了ans[i]表示最大子正方形边长小于i的方案数，ans[i+1]-ans[i]即是优美度恰为i的方案数。
现在考虑如何求ans[sz]:
在外层我们一行一行做，f(i,s)表示做了前i行，状态为s的方案数。
里面这个状态如何设计呢？关键就是状态要能体现出最大的子正方形边长是否超出限制sz。
我们就从当前i行的1~sz，st+1~2sz，2sz+1~3sz……往上看看能推多远即可。
所以状态就是一个sz进制数，每位表示从这里往上推能推多远。转移时超过sz-1就不合法了。
然后转移就是暴力枚举这行放哪些就好了。

一个小地方打错查了半天…不爽…

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MOD=1000000007;
int _test,n,f[15][2000],ans[15],p_sz[15],all_n;
int map[15];
int main(){
    freopen("hdu5079.in","r",stdin);
    freopen("hdu5079.out","w",stdout);
    scanf("%d",&_test);
    while(_test--){
        scanf("%d",&n);
        memset(map,0,sizeof map); all_n=1;
        for(int i=1;i<=n;i++){
            while(getchar()!='\n');
            for(int j=1;j<=n;j++){
                char ch=getchar(); if(ch=='o') (all_n*=2)%=MOD;
                map[i]|=(ch=='*')<<(j-1);
            }
        }
        memset(ans,0,sizeof ans); ans[0]=0; ans[1]=1; ans[n+1]=all_n;
        for(int sz=2;sz<=n;sz++){
            p_sz[0]=1; for(int i=1;i<=n;i++) p_sz[i]=p_sz[i-1]*sz;
            int alls=0; for(int i=1;i+sz-1<=n;i++) alls=alls*sz+(sz-1);
            memset(f,0,sizeof f); f[0][0]=1;
            for(int i=0;i<=n-1;i++)
             for(int s=0;s<=alls;s++) if(f[i][s]){
                for(int k=0;k<=(1<<n)-1;k++) if((map[i+1]&k)==0){
                    int nxt=0;
                    for(int tmp=k,j=1;j+sz-1<=n;j++,tmp>>=1){
                        int now=(((tmp&((1<<sz)-1))==(1<<sz)-1)?(s/p_sz[j-1]%sz+1):0); 
                        if(now>=sz){ nxt=-1; break; }
                        nxt+=now*p_sz[j-1];
                    }
                    if(nxt!=-1) f[i+1][nxt]=(f[i+1][nxt]+f[i][s])%MOD;
                }
             }
            for(int s=0;s<=alls;s++) ans[sz]=(ans[sz]+f[n][s])%MOD;
        }
        for(int i=0;i<=n;i++) printf("%d\n",((ans[i+1]-ans[i])%MOD+MOD)%MOD);
    }
    return 0;
}