Solid Dominoes Tilings HDU - 5731 (輪廓線 容斥 )

比較好玩的容斥,,,

如果不加限制的話,是一個經典的輪廓線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;
}

 

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章