SRM 573 div2

250...

500...

1000

　　比赛的时候没思路。倒是想到一种trick，想cha人来着，后来发现房间里没有人提交。。。

只想到了预处理出从一个点到另一个点走m步所要用到的情况数，发现dp空间复杂度太高。今天看了一下别人怎么预处理的，我还是太嫩啊。。。

dp[dx][dy][m]表示走m步坐标移动了(dx, dy)所有的情况数。

记忆化搜索一下。

const int MAXN = 55;
const int MOD = 1e9+7;

int dp[MAXN][MAXN][MAXN];

int dir[4][2] = {{0, 1}, {0, -1}, {1, 0}, {-1, 0}};

class WolfPackDivTwo {
public:
    int dfs(int dx, int dy, int m) {
        if(dp[dx][dy][m] != -1) return dp[dx][dy][m];
        int res = 0;
        if(m == 0) {
            res = (dx == 0 && dy == 0);
        } else {
            for(int i = 0; i < 4; ++i) {
                int xx = abs(dx + dir[i][0]);
                int yy = abs(dy + dir[i][1]);
                if(xx + yy > m - 1) continue;
                res = (res + dfs(xx, yy, m - 1))%MOD;
            }
        }
        return dp[dx][dy][m] = res;
    }

    int calc(vector <int> x, vector <int> y, int m) {
        int n = x.size(), tx, ty, i, ans = 0;
        CL(dp, -1);
        for(tx = x[0] - m; tx <= x[0] + m; ++tx) {
            for(ty = y[0] - m; ty <= y[0] + m; ++ty) {
                int tmp = 1;
                for(i = 0; i < n; ++i) {
                    int dx = abs(tx - x[i]);
                    int dy = abs(ty - y[i]);
                    if(dx + dy > m) {tmp = 0; continue;}
                    tmp = (tmp*dfs(dx, dy, m))%MOD;
                }
                ans = (ans + tmp)%MOD;
            }
        }
        return ans;
    }
};