【省内训练2018-12-21】Cards

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/qq_39972971/article/details/85215789

【思路要点】

  • 考虑最长公共子串的形式,应当由若干条长度为 iii1i-1 的链穿插而成,举个例子:
    1,4,2,5,3,y,b1,4,2,5,3,y,b
    a,x,1,4,2,5,3a,x,1,4,2,5,3
    就由一条长度为 44 的链 b321ab-3-2-1-a 和一条长度为 33 的链 y54xy-5-4-x 穿插而成。
  • 长度为 ii 的链对答案的贡献为 i1i-1 ,即若最长公共子串由 aa 个长度为 ii 的链和 bb 个长度为 i1i-1 的链拼接而成,答案应当为 a(i1)+b(i2)a(i-1)+b(i-2)
  • 考虑枚举 ii ,将所有置换环拆成若干长度为 iii1i-1 的链,并最大化每一个置换环对答案的贡献,这一点可以通过一个简单的背包来实现。
  • 背包部分的时间复杂度为 O(N2)O(N^2) ,单次询问我们需要枚举 ii ,并花费 O(Ni)O(\frac{N}{i}) 的时间计算答案,因此总时间复杂度为 O(N2+TNLogN)O(N^2+TNLogN)

【代码】

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 1005;
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
template <typename T> void chkmax(T &x, T y) {x = max(x, y); }
template <typename T> void chkmin(T &x, T y) {x = min(x, y); } 
template <typename T> void read(T &x) {
	x = 0; int f = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';
	x *= f;
}
template <typename T> void write(T x) {
	if (x < 0) x = -x, putchar('-');
	if (x > 9) write(x / 10);
	putchar(x % 10 + '0');
}
template <typename T> void writeln(T x) {
	write(x);
	puts("");
}
int n, m, a[MAXN], b[MAXN], c[MAXN], nxt[MAXN];
int dp[MAXN][MAXN]; bool vis[MAXN];
void init(int n) {
	dp[1][1] = 1;
	for (int i = 2; i <= n; i++) {
		for (int j = i; j <= n; j++) {
			if (j >= i) chkmax(dp[i][j], dp[i][j - i] + (i - 1));
			if (j >= i + 1) chkmax(dp[i][j], dp[i][j - i - 1] + i);
		}
	}
}
int main() {
	freopen("cards.in", "r", stdin);
	freopen("cards.out", "w", stdout);
	init(1e3);
	int T; read(T);
	while (T--) {
		int n; read(n), m = 0;
		for (int i = 1; i <= n; i++)
			read(a[i]);
		for (int i = 1; i <= n; i++)
			read(b[i]);
		for (int i = 1; i <= n; i++)
			nxt[a[i]] = b[i];
		memset(vis, false, sizeof(vis));
		for (int i = 1; i <= n; i++)
			if (!vis[i]) {
				int pos = i, cnt = 0;
				while (!vis[pos]) {
					vis[pos] = true;
					pos = nxt[pos];
					cnt++;
				}
				c[++m] = cnt;
			}
		sort(c + 1, c + m + 1);
		reverse(c + 1, c + m + 1);
		int ans = 0;
		for (int i = 1; i <= n; i++) {
			while (m >= 1 && c[m] < i) m--;
			int now = 0;
			for (int j = 1; j <= m; j++)
				now += dp[i][c[j]];
			chkmax(ans, now);
		}
		printf("%d\n", ans);
	}
	return 0;
}

没有更多推荐了,返回首页

私密
私密原因:
请选择设置私密原因
  • 广告
  • 抄袭
  • 版权
  • 政治
  • 色情
  • 无意义
  • 其他
其他原因:
120
出错啦
系统繁忙,请稍后再试