สมมติว่าเรามีตารางที่มีแถว H และคอลัมน์ W ซึ่งแต่ละตารางจะเป็นระเบียบเรียบร้อยหรือไม่เป็นระเบียบ เราสามารถวางโคมไฟบนช่องสี่เหลี่ยมที่เป็นระเบียบเป็นศูนย์หรือมากกว่าในตารางนี้ โคมไฟสามารถทำให้เซลล์สว่างขึ้นในแต่ละทิศทางทั้งสี่ - ขึ้น ลง ซ้าย และขวา - จนถึงจุดก่อนถึงขอบของตารางหรือสี่เหลี่ยมที่ไม่เป็นระเบียบเป็นครั้งแรก (เซลล์ที่ไม่เป็นระเบียบจะไม่ถูกจุด) ). หลอดไฟจะทำให้เซลล์สว่างขึ้น ในตารางถ้า G[i, j] เป็น '.' เซลล์นั้นเป็นระเบียบและเมื่อเป็น '#' จะไม่เป็นระเบียบ ให้ K เป็นจำนวนกำลังสองที่เป็นระเบียบเรียบร้อย มี 2^K วิธีในการวางโคมไฟทั้งหมด สมมติว่า สำหรับแต่ละวิธี 2^K เหล่านี้ จำนวนเซลล์ที่แบ่งเบาโดยหนึ่งหลอดขึ้นไปจะถูกคำนวณ เราต้องหาผลรวมของตัวเลขเหล่านั้นแบบโมดูโล 10^9 + 7
ดังนั้นหากอินพุตเป็นแบบ
. | . | # |
# | . | . |
แล้วผลลัพธ์จะเป็น 52
ขั้นตอน
เพื่อแก้ปัญหานี้ เราจะทำตามขั้นตอนเหล่านี้ -
m := 10^9 + 7 N = 2003 Define 2D arrays u, l, r, d of order N x N, and another list p with N^2 elements. h := row count of matrix w := column count of matrix tidy := 0 p[0] := 1 for initialize i := 1, when i <= h * w, update (increase i by 1), do: p[i] := p[i - 1] * 2 mod m for initialize i := 0, when i < h, update (increase i by 1), do: for initialize j := 0, when j < w, update (increase j by 1), do: u[i, j] := i l[i, j] := j if i is non-zero, then: u[i, j] := u[i - 1, j] if j is non-zero, then: l[i, j] := l[i, j - 1] if matrix[i, j] is same as '#', then: u[i, j] := i + 1 l[i, j] := j + 1 Otherwise (increase tidy by 1) for initialize i := h - 1, when i >= 0, update (decrease i by 1), do: for initialize j := w - 1, when j >= 0, update (decrease j by 1), do: d[i, j] := i r[i, j] := j if i < h - 1, then: d[i, j] := d[i + 1, j] if j < w - 1, then: r[i, j] := r[i, j + 1] if matrix[i, j] is same as '#', then: d[i, j] := i - 1 r[i, j] := j - 1 cnt := 0 for initialize i := 0, when i < h, update (increase i by 1), do: for initialize j := 0, when j < w, update (increase j by 1), do: if matrix[i, j] is same as '#', then: Ignore following part, skip to the next iteration src := d[i, j] + r[i, j] - u[i, j] - l[i, j] + 1 cnt := (cnt + (p[src] - 1) * p[tidy - src]) mod m return cnt
ตัวอย่าง
ให้เราดูการใช้งานต่อไปนี้เพื่อความเข้าใจที่ดีขึ้น -
#include <bits/stdc++.h> using namespace std; const int m = 1e9 + 7, N = 2003; int u[N][N], l[N][N], r[N][N], d[N][N], p[N * N]; int solve(vector<vector<char>> matrix){ int h = matrix.size(); int w = matrix[0].size(); int tidy = 0; p[0] = 1; for (int i = 1; i <= h * w; ++i) p[i] = p[i - 1] * 2 % m; for (int i = 0; i < h; ++i){ for (int j = 0; j < w; ++j){ u[i][j] = i; l[i][j] = j; if (i) u[i][j] = u[i - 1][j]; if (j) l[i][j] = l[i][j - 1]; if (matrix[i][j] == '#'){ u[i][j] = i + 1; l[i][j] = j + 1; } else ++tidy; } } for (int i = h - 1; i >= 0; --i){ for (int j = w - 1; j >= 0; --j){ d[i][j] = i; r[i][j] = j; if (i < h - 1) d[i][j] = d[i + 1][j]; if (j < w - 1) r[i][j] = r[i][j + 1]; if (matrix[i][j] == '#'){ d[i][j] = i - 1; r[i][j] = j - 1; } } } int cnt = 0; for (int i = 0; i < h; ++i){ for (int j = 0; j < w; ++j){ if (matrix[i][j] == '#') continue; int src = d[i][j] + r[i][j] - u[i][j] - l[i][j] + 1; cnt = (cnt + (p[src] - 1) * p[tidy - src]) % m; } } return cnt; } int main(){ vector<vector<char>> matrix = { { '.', '.', '#' }, { '#', '.', '.' } }; cout << solve(matrix) << endl; }
อินพุต
3, 2, { 1, 5, 9 }, { 2, 4, 2 }
ผลลัพธ์
52