比較好玩的容斥,,,
如果不加限制的話,是一個經典的輪廓線dp。
加上限制之後,考慮利用容斥求出不合法的方案數。
最開始想的時候,想的是求出所有的方案數,減去所有列不合法的方案數,減去所有行不合法但列合法的方案數。。
但是好像過去麻煩,可以考慮再求得時候就保證行合法,這樣就只需要減去列不合法的方案數。
只需要容斥一下,枚舉所有可能的列狀態,減去一列不合法的方案,加上兩列不合法方案,,,
那麼怎麼保證行合法?
還是減去行不合法的方案數。如果再進行一遍容斥的複雜度難以接受。
但是只需要枚舉最上面不合法的一行,此時,這一行上面的行一定合法,下面沒有要求。
那麼爲什麼我們不先這樣枚舉列,在這樣枚舉行呢?
如果我們先求出所以列不合法的狀態,再枚舉行時,上面的行是列合法的,下面的行是列合法的,但實際上上面下面是可以不合法的,只要對接起來合法即可。
所以我們必須先枚舉所以的列狀態,對一個固定的列狀態,再減去不合法的行。
#include <bits/stdc++.h>
using namespace std;
#define ms(a,b) memset(a,b,sizeof a)
#define add(a,b) (a+=((b)%mod+mod)%mod)%=mod
#define b(i) (1<<(i))
#define go(i,a,b) for(int i=(a);i<=(b);i++)
#define mod 1000000007
int dp[2][1<<17],bas[20][20],ans[20][20],t1[20],t2[20],p[20],n,m;
void init(){
go(n,1,16){ int cur=0; ms(dp,0); dp[cur][b(n)-1]=1;
go(i,1,n){
go(j,0,n-1){ cur^=1; ms(dp[cur],0);
go(k,0,b(n)-1) {
add(dp[cur][k^b(j)],dp[cur^1][k]);
if(j&&(k&b(j))&&!(k&b(j-1)))add(dp[cur][k|b(j-1)],dp[cur^1][k]);
}
} bas[n][i]=bas[i][n]=dp[cur][b(n)-1];
} }
go(n,1,16)go(k,0,b(n-1)-1){ int r=0,las=0;
go(i,0,n-2)if(k&b(i))p[++r]=i+1-las,las=i+1; p[++r]=n-las;
go(i,1,16){
int w=1;go(j,1,r)w=1ll*w*bas[p[j]][i]%mod;t1[i]=t2[i]=w;
go(j,1,i-1)add(t1[i],-1ll*t1[j]*t2[i-j]);
add(ans[i][n],t1[i]*(r&1?1:-1));
}
}
}
int main()
{
ios::sync_with_stdio(0);
init(); while(cin>>n>>m)cout<<ans[n][m]<< endl;
}