题意及思路：https://www.cnblogs.com/chaoswr/p/9460378.html

代码：

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const int maxn = 3010;
const LL mod = 1000000007;
LL dp[maxn][maxn];
char s[maxn][maxn];
int n, m;
LL ans[4];
void solve(int sx, int sy, LL ans[]) {
	if(s[sx][sy] == '#') return;
	memset(dp, 0, sizeof(dp));
	dp[sx][sy] = 1;
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= m; j++) {
			if(s[i][j] == '#') continue;
			if(s[i + 1][j] != '#' && i < n) dp[i + 1][j] = (dp[i + 1][j] + dp[i][j]) % mod;
			if(s[i][j + 1] != '#' && j < m) dp[i][j + 1] = (dp[i][j + 1] + dp[i][j]) % mod;
		}
	ans[0] = dp[n - 1][m], ans[1] = dp[n][m - 1];
}
int main() {
	scanf("%d%d", &n, &m);
	for (int i = 1; i <= n; i++)
		scanf("%s", s[i] + 1);
	solve(1, 2, ans);
	solve(2, 1, &ans[2]);
	printf("%lld\n", (ans[0] * ans[3] % mod - ans[1] * ans[2] % mod + mod) % mod);
}