1. 程式人生 > >[USACO 06NOV]Corn Fields

[USACO 06NOV]Corn Fields

logs swa des getchar max void 其中 i++ utc

Description

題庫鏈接

給你一個 \(0,1\) 矩陣,只準你在 \(1\) 上放物品;並且要滿足物品不能相鄰。允許空放,問方案數,取模。

\(1\leq n,m\leq 12\)

Solution

狀壓 \(DP\)

\(f_{i,t}\) 為處理到第 \(i\) 行時放物品的狀態為 \(t\) 的方案數。

若第 \(i\) 行的 \(1\) 的狀態為 \(a\) 。顯然由於必須放在 \(1\) 上,所以 \(a~and~t = a\) ,其中 \(and\) 為按位與;並且由於同行之間不能相鄰,故 \((t>>1)~and~t=0\)

枚舉的上一行狀態為 \(k\)

。顯然 \(t~and~k=0\)

復雜度為 \(O(n\cdot (2^n)^2)\)

Code

//It is made by Awson on 2018.2.25
#include <bits/stdc++.h>
#define LL long long
#define dob complex<double>
#define Abs(a) ((a) < 0 ? (-(a)) : (a))
#define Max(a, b) ((a) > (b) ? (a) : (b))
#define Min(a, b) ((a) < (b) ? (a) : (b))
#define Swap(a, b) ((a) ^= (b), (b) ^= (a), (a) ^= (b)) #define writeln(x) (write(x), putchar('\n')) #define lowbit(x) ((x)&(-(x))) using namespace std; const int yzh = 100000000; const int size = 1<<12; void read(int &x) { char ch; bool flag = 0; for (ch = getchar(); !isdigit(ch) && ((flag |= (ch == '-'
)) || 1); ch = getchar()); for (x = 0; isdigit(ch); x = (x<<1)+(x<<3)+ch-48, ch = getchar()); x *= 1-2*flag; } void print(LL x) {if (x > 9) print(x/10); putchar(x%10+48); } void write(LL x) {if (x < 0) putchar('-'); print(Abs(x)); } int n, m, a[15], x, f[15][size+5], bin[15]; void work() { read(n), read(m); bin[0] = 1; for (int i = 1; i <= 12; i++) bin[i] = bin[i-1]<<1; for (int i = 1; i <= n; i++) for (int j = 1; j <= m; j++) read(x), a[i] = ((a[i]<<1)|x); f[0][0] = 1; for (int i = 1; i <= n; i++) { for (int j = 0; j < bin[m]; j++) if ((j == 0 || ((j|a[i]) == a[i])) && (j&(j>>1)) == 0) for (int k = 0; k < bin[m]; k++) if ((k == 0 || ((k|a[i-1]) == a[i-1])) && (k&(k>>1)) == 0) if ((j&k) == 0) f[i][j] = (f[i][j]+f[i-1][k])%yzh; } int ans = 0; for (int i = 0; i < bin[m]; i++) ans = (ans+f[n][i])%yzh; writeln(ans); } int main() { work(); return 0; }

[USACO 06NOV]Corn Fields