くれなゐの雑記

例を上げて 自分で手を動かして学習できる入門記事を多めに書いています

Codeforces Div.2 682D Alyona and Strings

問題要約

LCSを求めるんだけど, 連続している文字列に限りがある(っていえばいいんだろうか)
図の[]の数分しかとれないLCSみたいな


方針

LCSっぽいしDPっぽい
n*m*kのDPで間に合いそうな感じ

2段階でdpを行う
一回目は連続して連なってる部分一致文字列を計算し、
二回目で答えを算出する クッソ説明が難しい

ソースコード

#include <iostream>
#include <queue>ソースコード
#include <map>
#include <list>
#include <vector>
#include <string>
#include <limits>
#include <cassert>
#include <fstream>
#include <cstring>
#include <bitset>
#include <iomanip>
#include <algorithm>
#include <functional>
#include <cstdio>
#include <ciso646>
#include <array>
#include <cmath>

using namespace std;

#define FOR(i,a,b) for (int i=(a);i<(b);i++)
#define RFOR(i,a,b) for (int i=(b)-1;i>=(a);i--)
#define REP(i,n) for (int i=0;i<(n);i++)
#define RREP(i,n) for (int i=(n)-1;i>=0;i--)
#define FOREACH(i, a) for(int i=0;i<a.size();++i)

#define inf 0x3f3f3f3f
#define CLEAR(a) a = decltype(a)()
#define MP make_pair
#define ALL(a) (a).begin(),(a).end()
#define pii pair<int ,int>
#define pcc pair<char,char>
#define pic pair<int,char>
#define pci pair<char,int>
#define VS vector<string>
#define VI vector<int>
#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define pi 2*acos(0.0)ソースコード
#define INFILE() freopen("in.txt","r",stdin)
#define OUTFILE() freopen("out.txt","w",stdout)
#define ll long long
#define ull unsigned long long
#define eps 1e-14

int dpa[1001][1001] = {};
int dp[1001][1001][11] = {};

signed main(void) {
	int N, M, K; cin >> N >> M >> K;
	string s1, s2; cin >> s1 >> s2;

	s1 = " " + s1;
	s2 = " " + s2;

	memset(dpa, 0, sizeof(dpa));
	memset(dp, 0, sizeof(dp));

	//連続して連なるLCSみたいなやつ
	FOR(i, 1, N + 1) {
		FOR(j, 1, M + 1) {
			if (s1[i] == s2[j]) dpa[i][j] = dpa[i - 1][j - 1] + 1;
			else dpa[i][j] = 0;
		}
	}

	for (int n = 1; n <= N; ++n) {
		for (int m = 1; m <= M; ++m) {
			for (int k = K; k >= 1; --k) {
				dp[n][m][k] = max(dp[n - 1][m][k], dp[n][m - 1][k]);
				if (dpa[n][m] > 0) {
					int len = dpa[n][m];
					dp[n][m][k] = max(dp[n][m][k], dp[n - len][m - len][k - 1] + len);
				}
			}
		}
	}

	cout << dp[N][M][K] << endl;
	return 0;
}