[gmoj 5102]【GDOI2017 day2】小学生语文题

在这里插入图片描述
思路:
首先我们很容易想到,每个字符一定只会移一次,这样才是最优的。
然后最关键的是,假设现在把b串(目标串)移动后符合了a串从后面开始的几位,那我们犯不着把已经匹配好的后面几位打乱,于是我们可以想到倒着dp。

设f[i][j]表示a串 i~n都匹配完成,b串用到了 j~n(有些可能不动,有些可能抽了出来)

考虑怎么转移:

  1. f[j][k]=f[j+1][k+1],那么就是说当前的这个位置a[i]=b[j],可以不用动上面的是不用抽出来的情况(不用贡献+1),所以要先更新状态,否则会有问题

  2. f[i][j]=f[i+1][j],那么就是说抽出来的可以把i填满,那么判断一下现在b串(j+1 ~ n)这个字符的数量是否比a串(i+1 ~ n)串这个字符的数量大。

  3. f[i][j]=f[i][j+1]表示将b[j+1]放弃合法匹配,而是跟以前的某个a中没能进行合法匹配的位置去匹配然而这个转移就要有条件了。并不是想放弃就能放弃的,前提是A中有能和它进行不合法匹配的位置

我们还要求移动方案。

首先找到那些不用动的点:
记录每个状态是由哪个状态推来的,如果f[i][j] 是由f[i+1][j+1]推来的,显然Bj和Ai就不会移动。

然后我们直接模拟一遍,从A串(目标串)的左边到右边,对于每一个Ai,如果它是被移过来的,直接在b串里随便找一个没被移过且要移动的,移过来即可,这里可以直接暴力把B串右移,反正都是O( n 2 n^2 n2

#include <cstdio>
#include <iostream>
#include <cstring>

using namespace std;

const int N = 2e3 + 10;
char s[N], t[N];
int m, n, bz[N], bz2[N], sum1[30][N], sum2[30][N], f[N][N], g[N][N][2], len;

int main()
{
	freopen("chinese.in", "r", stdin);
	freopen("chinese.out", "w", stdout);
	scanf("%d", &m);
	while(m--)
	{
		memset(g, 0, sizeof(g));
		memset(bz, 0, sizeof(bz)),memset(bz2, 0, sizeof(bz2));
		memset(sum1, 0, sizeof(sum1)), memset(sum2, 0, sizeof(sum2));
		memset(f, 127 / 3, sizeof(f));
		scanf("%s%s", s + 1, t + 1);
		n = strlen(s + 1);
		for(int i = n; i >= 1; i--)
		{
			for(int j = 0; j <= 25; j++)
				sum1[j][i] = sum1[j][i + 1], sum2[j][i] = sum2[j][i + 1];
			sum1[s[i] - 'a'][i]++; sum2[t[i] - 'a'][i]++;	
		}
		for(int i = 1; i <= n + 1; i++)
		{
			f[n + 1][i] = n - i + 1, 
			g[n + 1][i][0] = n + 1, g[n + 1][i][1] = i + 1;
		}
		for(int i = n; i >= 1; i--)
		{
			for(int j = i; j >= 1; j--)
			{	
				if(sum2[s[i] - 'a'][j] >= sum1[s[i] - 'a'][i] && f[i][j] > f[i + 1][j])
				{
					f[i][j] = f[i + 1][j];
					g[i][j][0] = i + 1, g[i][j][1] = j;
				}
				if(s[i] == t[j] && f[i + 1][j + 1] < f[i][j])
                    f[i][j] = f[i + 1][j + 1], g[i][j][0] = i + 1, g[i][j][1] = j + 1;
                if(f[i][j + 1] + 1 < f[i][j])
                    f[i][j] = f[i][j + 1] + 1, g[i][j][0] = i, g[i][j][1] = j + 1;
			}
		} 
		printf("%d\n", f[1][1]);
		int x = 1, y = 1;
		while(x <= n && y <= n)
		{
			int xx = g[x][y][0], yy = g[x][y][1];
            if(x + 1 == xx && y + 1 == yy)
                bz[x] = 1, bz2[y] = 1;
            x = xx; y = yy;
		}
		for(int i = 1; i <= n; i++)
		 if(!bz[i]) 
		 {
            for(int j = i; j <= n; j++)
			 if(!bz2[j] && s[i] == t[j])
			 {
                for(int k = j; k >= i + 1; k--)
                    bz2[k] = bz2[k - 1], t[k] = t[k - 1];
                printf("%d %d\n", j, i);
                break;
             }
        }
	}
	return 0;
} 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值