BZOJ1566 [NOI2009]管道取珠

这是一道思维复杂度很高的DP题

看题目，为什么是取两次序列一样呢？YY一下，其实等价于两个人一起取，最后序列一样。

然后就水了:

令f[i, j, k]表示取到第i个珠子，第一个人在1号管道取了j个珠子，第二个人在1号管道取了k个珠子时，他们取出的序列相等的方案数

于是真水了！！！（方程请自行脑补或看程序呗）

然后我那坑爹的编程能力，1个小时啊大爷的。。。

/**************************************************************
    Problem: 1566
    User: rausen
    Language: Pascal
    Result: Accepted
    Time:11920 ms
    Memory:3592 kb
****************************************************************/
 
const prime = 1024523;
 
var
  s1, s2 : array[0..600] of char;
  m, n, i, j, k, j1, k1 : longint;
  p, q : longint;
  f : array[0..1, 0..600, 0..600] of longint;
 
procedure reverse;
var
  s : ansistring;
  i : longint;
 
begin
  readln(s);
  for i := 1 to m do
    s1[m - i + 1] := s[i];
  readln(s);
  for i := 1 to n do
    s2[n - i + 1] := s[i];
end;
 
begin
  readln(m, n);
  reverse;
  f[0, 1, 1] := 1;
  for i := 1 to m + n + 1 do begin
    p := i mod 2;
    q := p xor 1;
    fillchar(f[p], sizeof(f[p]), 0);
    for j := 1 to m + 1 do
      for k := 1 to m + 1 do
        while f[q, j, k] >= prime do
          dec(f[q, j, k], prime);
    for j := 1 to m + 1 do
      for k := 1 to m + 1 do begin
        j1 := i - j + 1;
        k1 := i - k + 1;
        if (j1 > 0) and (k1 > 0) and (j1 <= n + 1) and (k1 <= n + 1) then begin
          if s1[j] = s1[k] then inc(f[p, j + 1, k + 1], f[q, j, k]);
          if s1[j] = s2[k1] then inc(f[p, j + 1, k], f[q, j ,k]);
          if s2[j1] = s1[k] then inc(f[p, j, k + 1], f[q, j ,k]);
          if s2[j1] = s2[k1] then inc(f[p, j, k], f[q, j, k]);
        end;
      end;
  end;
  writeln(f[q, m + 1, m + 1]);
end.

（p.s. 上面↑是以前打的，所以还是Pascal写的）