矩阵【NOIP2016提高A组模拟9.21】

题目

在麦克雷的面前出现了一个有n* m个格子的矩阵，每个格子用“.”或“#”表示，“.”表示这个格子可以放东西，“#”则表示这个格子不能放东西。现在他拿着一条1* 2大小的木棒，好奇的他想知道对于一些子矩阵，有多少种放木棒的方案

样例输入：
第一行包含 2 个正整数 n，m。
接下来 n 行每行包含 m 个字符“.”或“#”。
第n+1行包含1个正整数q，表示询问次数。
接下来q行每行包含4个正整数r1，c1，r2，c2，分别表示询问的子矩阵的左上格子和右下格子的位置。
5 8
….#..#
.#……
##.#….
##..#.##
……..
4
1 1 2 3
4 1 4 1
1 2 4 5
2 5 5 8

样例输出：
4
0
10
15

数据范围：
30%：q<=100
100%：q<=10^5,1<=n,m<=500

---

剖解题目

。。。。。。

---

思路

方案数一般也可以是dp吧。。。。。

---

解法

dp
设$f_{i,j}$表示右下角格子左边为（i，j），左上角位置是（1，1）的矩阵的方案数。
转移就是

$$f_{i,j}=f_{i-1,j}+f_{i,j-1}-f_{i-1,j-1}+k$$
k表示他与左边与上面的格子是否能够成新的方案的数量，只有0，1，2三个值。
至于求和就稍微复杂点，先如同一般的矩阵那样，右下角减左边和上面在加左上角，但这里要把这个矩阵边上的一个点构成的方案数要加回来，这个自己想一想吧。

---

代码

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)

using namespace std;

const int maxn=505;
int f[maxn][maxn],n,m,a[maxn][maxn],sumh[maxn][maxn],suml[maxn][maxn];
char s[maxn];

int main()
{
    //freopen("T.in","r",stdin);
    scanf("%d%d",&n,&m); 
    fo(i,1,n){
        scanf("%s",&s);
        fo(j,1,m){
            a[i][j]=s[j-1];
            f[i][j]=f[i-1][j]+f[i][j-1]-f[i-1][j-1];
            sumh[i][j]=sumh[i][j-1]; suml[j][i]=suml[j][i-1];
            if (a[i][j]==46){
                if (a[i-1][j]==46) {
                    ++f[i][j]; ++suml[j][i];
                }
                if (a[i][j-1]==46) {
                    ++f[i][j]; ++sumh[i][j];
                }
            }
        }
    }
    int T;
    scanf("%d",&T);
    fo(i,1,T){
        int x1,y1,x2,y2;
        scanf("%d%d%d%d",&x2,&y2,&x1,&y1);
        int ans1=f[x1][y1]-f[x2][y1]-f[x1][y2]+f[x2][y2];
        int ans2=sumh[x2][y1]-sumh[x2][y2]+suml[y2][x1]-suml[y2][x2];
        int ans=ans1+ans2;
        printf("%d\n",ans);
    }
}