Codeforces Round 646 (Div. 2) F. Rotating Substrings(基础dp 最长公共子序列lcs思想+hall定理思想)

题目

给出一个长度为n(n<=2e3)的小写字母字符串s，和等长小写字母字符串t，

你可以执行以下操作任意次，每次操作，可以选择s的一个子区间[l,r]，将s[r]拽到s[l]左边，

具体来说，操作前是[s[l],s[l+1],...,s[r-1],s[r]]，操作后是[s[r],s[l],s[l+1],...,s[r-1]]

求将s换成t的最小操作次数，如果无法交换，输出-1

实际为t(t<=2e3)组样例，保证sumn不超过2e3

思路来源

b站heyuhhh

Submission #101216700 - Codeforces

题解

第一个做法，第二个做法

第一个做法是，

先判一下不合法的情况，两个串排序之后如果不等说明无解，

否则肯定有解，把s中第一个t1字符拽到s1的位置，然后考虑[2,n]，以此类推，就是一组解

考虑s和t的最长公共子序列，类似最长公共子序列的转移，

最长公共子序列是不需要换的，一一对应的，其他位置都是需要换的，

s[i]和t[j]能匹配，当且仅当，匹配了这个字符之后，后面的每种字符都还够用，

由于s可以把字符换到左边，根据类似hall定理的结论，就是s后缀中每种字符都不少于t后缀

第二个做法是，

记忆化搜索，dfs(i,j)表示当前倒着考虑到s的第i个字符，t的第j个字符时的代价

然后考虑三种转移，

①s[i]=t[j]，直接匹配上，继续考虑(i-1,j-1)

②把s[i]换到左边，只是暂且拿出来了，换到哪里并没有确定，转移到(i-1,j)，代价为1

③把t[j]和s之前拿出来的相同字母匹配，

能这样匹配，当且仅当在[j+1,n]中t[j]的出现次数cnt2，小于在[i+1,n]中相同字母的出现次数cnt1

只有cnt1>cnt2，说明靠后的s字符才需要被换到前面来，而其代价，在前面的转移中已经加过了

这种同时利用下标和总量限制的转移，感觉还是很值得学习的

类似括号匹配，仿佛②是加上一个(，而③是加上一个)

代码1

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define per(i,a,b) for(int i=(a);i>=(b);--i)
typedef long long ll;
typedef double db;
typedef pair<int,int> P;
#define fi first
#define se second
#define dbg(x) cerr<<(#x)<<":"<<x<<" ";
#define dbg2(x) cerr<<(#x)<<":"<<x<<endl;
#define SZ(a) (int)(a.size())
#define sci(a) scanf("%d",&(a))
#define pb push_back
#define all(a) a.begin(),a.end()
#define pt(a) printf("%d",a);
#define pte(a) printf("%d\n",a)
#define ptlle(a) printf("%lld\n",a)
#define debug(...) fprintf(stderr, __VA_ARGS__)
std::mt19937_64 gen(std::chrono::system_clock::now().time_since_epoch().count());
ll get(ll l, ll r) { std::uniform_int_distribution<ll> dist(l, r); return dist(gen); }
const int N=2e3+10,M=26,INF=0x3f3f3f3f;
int T,n,dp[N][N],suf[M][N],suf2[M][N];
char s[N],t[N],ss[N],tt[N];
bool no(){
	strcpy(ss+1,s+1);
	strcpy(tt+1,t+1);
	sort(ss+1,ss+n+1);
	sort(tt+1,tt+n+1);
	return strcmp(ss+1,tt+1);
}
int main(){
	sci(T);
	while(T--){
		sci(n);
		scanf("%s%s",s+1,t+1);
		if(no()){
			puts("-1");
			continue;
		}
		rep(i,0,n+1){
			rep(j,0,n+1){
				dp[i][j]=0;
			}
			rep(j,0,3){
				suf[j][i]=suf2[j][i]=0;
			}
		}
		per(i,n,1){
			rep(j,0,25){
				suf[j][i]=suf[j][i+1];
				suf2[j][i]=suf2[j][i+1];
			}
			suf[s[i]-'a'][i]++;
			suf2[t[i]-'a'][i]++;
		}
		rep(i,0,n){
			rep(j,0,n){
				if(i && j && s[i]==t[j]){
					bool ok=1;
					rep(k,0,25){
						ok&=(suf[k][i+1]>=suf2[k][j+1]);
					}
					if(ok)dp[i][j]=max(dp[i][j],dp[i-1][j-1]+1);
				}
				if(j)dp[i][j]=max(dp[i][j],dp[i][j-1]);
				if(i)dp[i][j]=max(dp[i][j],dp[i-1][j]);
			}
		}
		printf("%d\n",n-dp[n][n]);
	}
	return 0;
}

代码2

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define per(i,a,b) for(int i=(a);i>=(b);--i)
typedef long long ll;
typedef double db;
typedef pair<int,int> P;
#define fi first
#define se second
#define dbg(x) cerr<<(#x)<<":"<<x<<" ";
#define dbg2(x) cerr<<(#x)<<":"<<x<<endl;
#define SZ(a) (int)(a.size())
#define sci(a) scanf("%d",&(a))
#define pb push_back
#define all(a) a.begin(),a.end()
#define pt(a) printf("%d",a);
#define pte(a) printf("%d\n",a)
#define ptlle(a) printf("%lld\n",a)
#define debug(...) fprintf(stderr, __VA_ARGS__)
std::mt19937_64 gen(std::chrono::system_clock::now().time_since_epoch().count());
ll get(ll l, ll r) { std::uniform_int_distribution<ll> dist(l, r); return dist(gen); }
const int N=2e3+10,M=26,INF=0x3f3f3f3f;
int T,n,dp[N][N],suf[M][N],suf2[M][N];
char s[N],t[N];
int dfs(int i,int j){
	//if(i<=0 && j)return -1;
	if(!i && !j)return 0;
	if(~dp[i][j])return dp[i][j];
	int &ans=dp[i][j];
	ans=INF;
	if(i){
		ans=min(ans,dfs(i-1,j)+1);//让i和前面的换 把i挖走
	}
	int v=t[j]-'a';
	if(j && suf[v][i+1]>suf2[v][j+1]){//让i后面的和j换 前提是后边的多
		ans=min(ans,dfs(i,j-1));
	}
	if(i && j){//匹配 不换 直接考虑下一个
		if(s[i]==t[j])ans=min(ans,dfs(i-1,j-1));
	}
	return ans;
}
int main(){
	sci(T);
	while(T--){
		sci(n);
		scanf("%s%s",s+1,t+1);
		rep(i,0,n+1){
			rep(j,0,n+1){
				dp[i][j]=-1;
			}
			rep(j,0,3){
				suf[j][i]=suf2[j][i]=0;
			}
		}
		per(i,n,1){
			rep(j,0,25){
				suf[j][i]=suf[j][i+1];
				suf2[j][i]=suf2[j][i+1];
			}
			suf[s[i]-'a'][i]++;
			suf2[t[i]-'a'][i]++;
		}
		int ans=dfs(n,n);
		if(ans==INF)ans=-1;
		printf("%d\n",ans);
	}
	return 0;
}