[分块 + DP] Candies and Stones

【题目描述】
小A和小B玩游戏，有n堆石子和m堆糖果，小B先动。
小A的操作，取走一堆石子或一堆糖果。
小B的操作，若截止到此轮取走了i堆石子，j堆糖果，则分值加上(xi+yj)%p
要求：直至全部拿完，分值最大，并输出方案
【输入文件】
第一行三个数n,m,p,意义如题所述
第二行n个数x0~xn-1
第三行m个数y0~ym-1
【输出文件】
第一行一个数表示最大分值
第二行输出方案，S表示拿石头（stone）,C表示拿糖果（candy）
【样例输入1】
2 2 10
0 0
0 1
【样例输出1】
2
SC
【样例输入2】
3 3 10
0 2 0
0 0 2
【样例输出2】
10
CSSC
【数据规模】
1≤n,m≤20000, 1≤p≤109
0≤xi≤20000 0≤yi≤20000

这是考试时的一道题，原题见 CF 101E 。鉴于 15s 的时限，平方的方程可以秒出了；但是只有 45M 的空间限制，可以计算一下，如果全部强压到 long long 里面去，空间刚刚好 48M。于是就被卡住了。。。。。怎么办呢？不知道是谁想出来的，竟然分块记状态：DP 中对于每一块记初始状态，DP 完之后从后往前对于每一个块分别倒着在 DP 一次，由于可以知道最终状态从何转入，每次经过一个块都可以将路径求出。于是，这道看似很难很难的题就被解决了。

然后分块的时候有很多细节。。。。非常坑爹。。。。以后写这样的题目一定要细心。。

Code :

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>

#define maxn 20005
#define maxl 150

int n, m, p, q;
int x[maxn], y[maxn], s[maxn], t[maxn];
int a[maxl + 1][maxn], f[maxl + 1][maxn], g[maxl + 1][maxn];
char ans[maxn * 2], * tot = ans;

inline void getit(int s, int r)
{
	for (int i = 1; i <= r; ++ i)
		for (int j = 1; j <= m; ++ j)
			if (j == 1 || f[i - 1][j] > f[i][j - 1])
				f[i][j] = f[i - 1][j] + (x[s + i] + y[j]) % p, g[i][j] = 0;
			else
				f[i][j] = f[i][j - 1] + (x[s + i] + y[j]) % p, g[i][j] = 1;
}

inline void getans(int r, int & j)
{
	for (int i = r; i >= 1; )
		if (g[i][j] == 0) -- i, * (++ tot) = 'C';
		else -- j, * (++ tot) = 'S';
}

int main()
{
	freopen("E.in", "r", stdin);
	freopen("E.out", "w", stdout);

	scanf("%d%d%d", & n, & m, & p);
	for (int i = 1; i <= n; ++ i) scanf("%d", x + i);
	for (int i = 1; i <= m; ++ i) scanf("%d", y + i);
	memset(a[0], 0xEE, m + 1 << 2), a[0][1] = 0;
	for (int i = 1; i <= n; i += maxl)
	{
		int j = std :: min(i + maxl - 1, n);
		memcpy(f[0], a[q], m + 1 << 2);
		++ q, s[q] = i, t[q] = j, getit(i - 1, j - i + 1);
		memcpy(a[q], f[j - i + 1], m + 1 << 2);
	}
	printf("%d\n", f[t[q] - s[q] + 1][m]);
	int kk = m;
	getans(t[q] - s[q] + 1, kk);
	for (int i = q - 1; i; -- i)
	{
		memcpy(f[0], a[i - 1], m + 1 << 2);
		getit(s[i] - 1, t[i] - s[i] + 1);
		getans(t[i] - s[i] + 1, kk);
	}
	for (-- tot; tot > ans; -- tot) putchar(* tot);

	return 0;
}