WC模拟(1.10) T3 Subsequence

标签: NOI DP 复杂度分析
111人阅读 评论(0) 收藏 举报
分类:

Subsequence

题目背景:

1.10 WC模拟T3

分析:DP + 复杂度分析

 

这个题感觉难点在复杂度分析,看似总状态数n * m转移复杂度k,但是实际上,可以证明得到答案数与k的乘积是小于(n + m)的,证明:假设当前匹配到A的第j位,B的第i位,那么第j + 1位到第j + k - 1位,那么一定有一个数没有出现过,那么选择这个数就会让i + j增加k,所以每一次操作至少能增加k,所以最后的答案不超过(n + m) / k,所以直接按照O(n3)的方式DP就可以了,实际复杂度为O(n2),定义dp[i][j]表示,当前枚举第i位,在A串中匹配到jdp[i][j]表示在B串中匹配到的最后的位置。预处理next_a[i][j]表示第i位之后第一个j出现的位置,最后结尾设为n + 1,那么当dp[i][n + 1] = m + 1时,说明i就是可行的长度了。

 

Source:

 

/*
	created by scarlyw
*/
#include <cstdio>
#include <string>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <cmath>
#include <cctype>
#include <vector>
#include <set>
#include <queue>
#include <ctime>
#include <bitset>

inline char read() {
	static const int IN_LEN = 1024 * 1024;
	static char buf[IN_LEN], *s, *t;
	if (s == t) {
		t = (s = buf) + fread(buf, 1, IN_LEN, stdin);
		if (s == t) return -1;
	}
	return *s++;
}

///*
template<class T>
inline void R(T &x) {
	static char c;
	static bool iosig;
	for (c = read(), iosig = false; !isdigit(c); c = read()) {
		if (c == -1) return ;
		if (c == '-') iosig = true;	
	}
	for (x = 0; isdigit(c); c = read()) 
		x = ((x << 2) + x << 1) + (c ^ '0');
	if (iosig) x = -x;
}
//*/

const int OUT_LEN = 1024 * 1024;
char obuf[OUT_LEN], *oh = obuf;
inline void write_char(char c) {
	if (oh == obuf + OUT_LEN) fwrite(obuf, 1, OUT_LEN, stdout), oh = obuf;
	*oh++ = c;
}

template<class T>
inline void W(T x) {
	static int buf[30], cnt;
	if (x == 0) write_char('0');
	else {
		if (x < 0) write_char('-'), x = -x;
		for (cnt = 0; x; x /= 10) buf[++cnt] = x % 10 + 48;
		while (cnt) write_char(buf[cnt--]);
	}
}

inline void flush() {
	fwrite(obuf, 1, oh - obuf, stdout);
}

/*
template<class T>
inline void R(T &x) {
	static char c;
	static bool iosig;
	for (c = getchar(), iosig = false; !isdigit(c); c = getchar())
		if (c == '-') iosig = true;	
	for (x = 0; isdigit(c); c = getchar()) 
		x = ((x << 2) + x << 1) + (c ^ '0');
	if (iosig) x = -x;
}
//*/

const int MAXN = 4000 + 10;

int n, m, k;
int dp[MAXN][MAXN], a[MAXN], b[MAXN];
int next_a[MAXN][MAXN], next_b[MAXN][MAXN];

inline void read_in() {
	R(n), R(m), R(k);
	for (int i = 1; i <= n; ++i) R(a[i]);
	for (int i = 1; i <= m; ++i) R(b[i]);
	for (int i = 1; i <= k; ++i) {
		next_a[n][i] = next_a[n + 1][i] = n + 1;
		next_b[m][i] = next_b[m + 1][i] = m + 1;
	}
	for (int i = n - 1; i >= 0; --i) {
		for (int j = 1; j <= k; ++j)
			next_a[i][j] = next_a[i + 1][j];
		next_a[i][a[i + 1]] = i + 1;
	}
	for (int i = m - 1; i >= 0; --i) {
		for (int j = 1; j <= k; ++j)
			next_b[i][j] = next_b[i + 1][j];
		next_b[i][b[i + 1]] = i + 1;
	}
}

inline void solve() {
	memset(dp, -1, sizeof(dp));
	dp[0][0] = 0;
	for (int i = 0; ; ++i) {
		for (int j = 0; j <= n + 1; ++j) {
			if (dp[i][j] == -1) continue ;
			if (j == n + 1 && dp[i][j] == m + 1) std::cout << i, exit(0);
			for (int c = 1; c <= k; ++c)
				dp[i + 1][next_a[j][c]] = std::max(dp[i + 1][next_a[j][c]], 
					next_b[dp[i][j]][c]);
		}
	}
}

int main() {
	freopen("subsequence.in", "r", stdin);
	freopen("subsequence.out", "w", stdout);
	read_in();
	solve();
	return 0;
}


查看评论

WC模拟(1.10) T2 Permutation

结论推理
  • scar_lyw
  • scar_lyw
  • 2018年01月17日 20:23
  • 105

WC模拟(1.10) T1 Function

莫比乌斯反演
  • scar_lyw
  • scar_lyw
  • 2018年01月11日 16:21
  • 132

WC模拟(1.4) T3 字符串

AC自动机 + DP
  • scar_lyw
  • scar_lyw
  • 2018年01月11日 15:35
  • 111

2017.1.17WC模拟赛总结

boom zero了,很不开心,总结一下问题。 T1boom zero了,原因是常数太大自己又没有去卡。 T2boom zero了,原因是没看好题,不符合它输出要求。 T3boom zero了,...
  • WerKeyTom_FTD
  • WerKeyTom_FTD
  • 2017年01月17日 14:47
  • 512

WC模拟(1.8) T3 送你一朵圣诞树

并查集 + 堆 + 贪心
  • scar_lyw
  • scar_lyw
  • 2018年01月11日 16:09
  • 171

WC模拟(1.6) T3 送你一棵圣诞树

状压DP + 子集枚举
  • scar_lyw
  • scar_lyw
  • 2018年01月11日 15:51
  • 110

WC模拟(1.2) T3 取石子

结论分析
  • scar_lyw
  • scar_lyw
  • 2018年01月03日 14:13
  • 151

Wannafly模拟赛3 贝伦卡斯泰露(搜索+剪枝)

题目链接:点击打开链接 思路:搜索+剪枝 // 贝伦卡斯泰露 运行/限制:3ms/1000ms #include #include #include using namespace std;...
  • ccDLlyy
  • ccDLlyy
  • 2017年10月14日 18:14
  • 201

bzoj 2509 送分题 && SNOI省选模拟赛Round3 T3 正方形square 树状数组

题目链接:戳这里2509: 送分题Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 62  Solved: 26[Submit][Status][Disc...
  • qq_39791208
  • qq_39791208
  • 2018年01月17日 21:32
  • 86

省赛模拟 水题 挑战密室

挑战密室 时间限制:1000 ms  |  内存限制:65535 KB 难度:4 描述 R组织的特工Dr. Kong 为了寻找丢失的超体元素,不幸陷入WTO密室。Dr. Kong必须...
  • FengTwoYear
  • FengTwoYear
  • 2015年11月06日 14:17
  • 297
    个人资料
    专栏达人 持之以恒
    等级:
    访问量: 2万+
    积分: 2169
    排名: 2万+
    博客专栏
    最新评论