TOJ 5020: Palindromic Paths

5020: Palindromic Paths 

Time Limit(Common/Java):10000MS/30000MS     Memory Limit:65536KByte
Total Submit: 8            Accepted:4

Description

Given an N×N grid of fields (1≤N≤500), each labeled with a letter in the alphabet. For example:

ABCD

BXZX
CDXB
WCBA

Each day, Susa walks from the upper-left field to the lower-right field, each step taking her either one field to the right or one field downward. Susa keeps track of the string that she generates during this process, built from the letters she walks across. She gets very disoriented, however, if this string is a palindrome (reading the same forward as backward), since she gets confused about which direction she had walked.

Please help Susa determine the number of distinct routes she can take that correspond to palindromes. Different ways of obtaining the same palindrome count multiple times. Please print your answer modulo 1,000,000,007.

Input

The first line of input contains N, and the remaining N lines contain the N rows of the grid of fields. Each row contains N characters that are in the range A...Z.

Output

Please output the number of distinct palindromic routes Susa can take, modulo 1,000,000,007.

Sample Input

4
ABCD
BXZX
CDXB
WCBA

Sample Output

12

Hint

Susa can make the following palindromes

1 x "ABCDCBA"

1 x "ABCWCBA"

6 x "ABXZXBA"

4 x "ABXDXBA"

Source

USACO 2015 US Open

一道不错的枚举+滚动数组，美滋滋，f[i][j][k]表示第一个点在第i行，第2个点在第j行都走了k步的方案数

#include <stdio.h>
#include <algorithm>
using namespace std;
typedef __int64 ll;
const int mod=1e9+;
char s[][];
ll dp[][][];
int main() {
    int n;
    scanf("%d",&n);
    for(int i=; i<=n; i++)
        scanf("%s",s[i]+);
    int now=,pre=;
    if(s[][]!=s[n][n]) {
        return ,printf("0\n");
    }
    dp[][n][pre]=;
    for(int k=; k<=n; k++) {
        for(int i=; i<=k; i++)
            for(int j=n; j>=i&&j>=n-k+; j--) {
                if(s[i][k-i+]==s[j][n-k+n-j+])
                    dp[i][j][now]=(dp[i-][j][pre]+dp[i][j][pre]+dp[i][j+][pre]+dp[i-][j+][pre])%mod;
                else dp[i][j][now]=;
            }
        swap(now,pre);
    }
    ll ans=;
    for(int i=; i<=n; i++) {
        ans=(ans+dp[i][i][pre])%mod;
    }
    printf("%lld",ans);
    return ;
}