7.23模拟赛总结 [数据结构优化dp] + [神奇建图]

复盘

浅复盘下吧…

7:40 开题

看 T1 ,起初以为和以前某道题有点像,子序列划分,注意到状态数很少,搜出来所有状态然后 dp,然后发现这个 T1 和那个毛关系没有

浏览了一下,感觉 T2 题面又臭又长,意思也很恶心;T3 树上 ds 看上去很亲切;T4 一眼就不太可做

回来看 T1,发现很 sb,每次肯定贪心选最优的后缀,至于选 m 的限制,填满一段后再贪心往前选就行

8:04 码完了,小样例顺利过掉,一测大样例发现小数点后面不对,想了想应该不会是精度问题,可能做法假了

然后 不断证明结论+修改一些毫无影响的细节,仍然和大样例不上,最后决定把题面重读一遍,知道看到了 “相对误差”,6,突然想到我学过高中物理。。。到 8:50 才交上去

快开 T2,朴素的想法当然是 二维状态记 Max,但这个显然不太有优化空间(这种状态设计只有 i − 1 i-1 i1 i i i 只有 O ( 1 ) O(1) O(1) 个决策需要单点修改时才能优化)

那么很快想到改成 f i f_i fi 一维状态,钦定 i i i 必选,往前找一个合法的,发现可以接上去的区间可以预处理出来,那么只需找 [ l , r ] [l,r] [l,r] 内小于 a i a_i ai 的最优决策… 三维偏序?树套树警告!

算了,先打暴力,竟然有 O ( 1 ) p t s O(1) pts O(1)pts 的高分,然鹅被一个 sb 细节卡住一直到 9:30 仍没调过,决定先看 T3

一眼发现是 sb 题,直接树上主席树就做完了,直接决定开码,10:30 过了小样例

然鹅大样例又挂了,反复调根本不明白怎么回事,一度以为自己主席树板子假了

最后突然想着把数组开大一倍,卧槽对了,不爆 RE 直接 WA 是吧,浪费了高达 O ( n ) O(n) O(n) 的时间

11:00 了,先把 T1 暴力调了吧,发现弱智错误… 然后看看能不能写点性质,得到了高达 5 p t s 5pts 5pts 的性质分

最后想 T2 感觉写不完,决定开 T4 暴力

最后获得了高达 7 p t s 7pts 7pts 的暴力分,赢!

所以 100+53+100+7=260 , rk O ( 2 n ) O(2^n) O(2n)

题解

T2

在这里插入图片描述
在这里插入图片描述
其实离正解不远了,回顾一下转移:

f i = m a x ( f j + 1 ) f_i=max(f_j+1) fi=max(fj+1),其中 j < i , a j < a i j<i,a_j<a_i j<i,aj<ai i , j i,j i,j 可以通过小于 a i a_i ai 的数 “连接” 起来

直接三维偏序做也没错,但是本题可以有更多性质

注意到:能转移的决策满足 a j < a i a_j<a_i aj<ai,能连接其它数的数也是 a j < a i a_j<a_i aj<ai ,二者可以同时维护

那么我们转到值域上做,从小到大插入新决策

每次处理询问时相当于要找连续段,这个可以 s e t set set 做(更好做法是并查集,但是我不会);然后询问区间 Max,线段树维护

遇到偏序问题可以尝试 在最外层换一维关键字排序,会有意想不到的效果

#include<bits/stdc++.h>
using namespace std ;

typedef long long LL ;
const int N = 3e5+100 ;

int n , D , a[N] ;
int f[N] , rhk[N] ;
struct Segtree
{
	int l , r , Max ;
}t[4*N] ;
bool cmp( int x , int y )
{
	return a[x] < a[y] ;
}
void build( int p , int l , int r )
{
	t[p].l = l , t[p].r = r ;
	if( l == r ) {
		return ;
	}
	int mid = ( l + r ) >> 1 ;
	build( p<<1 , l , mid ) ; build( p<<1|1 , mid+1 , r ) ;  
}
int ask( int p , int l , int r )
{
	if( l <= t[p].l && t[p].r <= r ) {
		return t[p].Max ;
	}
	int mid = ( t[p].l + t[p].r ) >> 1 , Max = 0 ;
	if( l <= mid ) Max = max( Max , ask(p<<1,l,r) ) ;
	if( r > mid ) Max = max( Max , ask(p<<1|1,l,r) ) ;
	return Max ;
}
void modify( int p , int x , int d )
{
	if( t[p].l == t[p].r ) {
		t[p].Max = max( t[p].Max , d ) ;
		return ;
	}
	int mid = ( t[p].l + t[p].r ) >> 1 ;
	if( x <= mid ) modify( p<<1 , x , d ) ;
	else modify( p<<1|1 , x , d ) ;
	t[p].Max = max( t[p<<1].Max , t[p<<1|1].Max ) ;
}
struct node
{
	int l , r ;
	friend bool operator < ( node x , node y ) {
		return x.r < y.r ;
	}
};
set<node> s ; // 维护连续段 
int query( int x ) // 查 x 可用决策 
{
	int P = max(1,x-D) ;
	auto it = s.lower_bound((node){0,P}) ;
	if( it == s.end() ) {
		return 0 ;
	}
	return ask( 1 , min(it->l,P) , x-1 ) ;
}
void Insert( int x )// 插入决策 x 
{
	modify( 1 , x , f[x] ) ;
	if( s.empty() ) {
		s.insert((node){x,x}) ;
		return ;
	}
	auto it = s.lower_bound((node){0,x}) ;
	int L = x , R = x , fg = 0 ;
	if( it != s.end() ) {
		if( x < it->l ) {
			if( x+D >= it->l ) {// 接上了
				R = it->r ;
				fg = 1 ;
			} 
			if( it != s.begin() ) {
				auto iit = it ;
				iit -- ;
				if( iit->r + D >= x ) {
					L = iit->l ;
					s.erase( iit ) ;
				}
			}
			if( fg ) s.erase( it ) ;
			s.insert((node){L,R}) ;
		}
		return ;
	}
	it -- ;
	if( it->r+D >= x ) {
		L = it->l ;
		s.erase( it ) ;
	}
	s.insert((node){L,x}) ;
}

int main()
{
	scanf("%d%d" , &n , &D ) ;
	for(int i = 1 ; i <= n ; i ++ ) {
		scanf("%d" , &a[i] ) ; 
		rhk[i] = i ;
	}
	// 转移有两维限制,按值域做,去掉一维
	sort( rhk+1 , rhk+n+1 , cmp ) ;
	build( 1 , 1 , n ) ;
	int ans = 0 ;
	for(int i = 1 ; i <= n ; i ++ ) {
		int j = i ;
		f[rhk[i]] = query( rhk[i] ) + 1 ;
		while( j+1 <= n && a[rhk[j]] == a[rhk[j+1]] ) {
			j ++ ;
			f[rhk[j]] = query( rhk[j] ) + 1 ;
		}
		for(int k = i ; k <= j ; k ++ ) {
			Insert( rhk[k] ) ;
		}
		i = j ;
	}
	for(int i = 1 ; i <= n ; i ++ ) {
		ans = max( ans , f[i] ) ;
	}
	printf("%d" , ans ) ;
	return 0 ;
}

其实很好写

T4

神奇题
在这里插入图片描述
尝试 d p dp dp ,发现根本做不了…

根据我们的经验,遇到十分之抽象且题面又臭又长又不知所云的题目大多是图论,所以这题考虑建图

老师讲的一个建图出发点:注意到最终序列中每个值 要么由 a a a 提供,要么由 b b b,这种 “由 A 或由 B” 的问题是有二分图的样子的,考虑分左右两部,对于每个 i i i,将左部 a i a_i ai 与 右部 b i b_i bi 连边

我理解的:题目中给的 划分 实际上约束是很多的,我们考虑最初选所有 a a a 中的集合,然后用 b b b 中的若干个集合替换,这样 替换者的并 和 被替换者的并 必须是一样的,想换掉某个子集 a [ ] a[] a[] 的话,必须选择 所有 包含 b b b 中对应位置 的区间,然后对应回来 a a a 中区间也必须选… 如此往复(口胡)

我们尝试形式化的去描述这个东西,发现正可以用上面那种建边方式来描述 ,一个连通块就代表了 可以相互替换的一组 a [ ] , b [ ] a[],b[] a[],b[]

对于 K = 0 K=0 K=0,我们直接对于每个连通块取较小的那一组替换

对于 K = 1 K=1 K=1,考虑枚举边,分讨一下:

若断的是非割边,然后连到另一个连通块,这样贡献的变化是 m i n ( a , b ) + m i n ( c , d ) → m i n ( a + c , b + d ) min(a,b)+min(c,d)\to min(a+c,b+d) min(a,b)+min(c,d)min(a+c,b+d)

显然有不等式 m i n ( a , b ) + m i n ( c , d ) ≤ m i n ( a + c , b + d ) min(a,b)+min(c,d)\leq min(a+c,b+d) min(a,b)+min(c,d)min(a+c,b+d),所以贡献不会变小,断开一定不优

这个不等式也告诉我们,要尽可能拆连通块,尽可能不合并连通块

若断割边,我们最好的当然是连回本身,只拆,不合并

注意!这里有一个 conner,如果一个连通块只有两个节点 u , v u,v u,v 怎么办?我们显然没办法把 u u u 向自己连边

然而有这么一个事情:不妨设 u u u 为左部点,那么一定可以把 u u u 连到一个 左部点数不小于右部点的连通块上(找不到的话考虑右部点 v v v ),当 n ≥ 1 n\geq 1 n1 时,这样的连通块一定能找到

n = 1 n=1 n=1,直接特判 答案为 1 1 1 即可

对于 K = 2 K=2 K=2,与上面相反,需要 尽可能合并连通块,尽可能不拆连通块

对于一个连通块,如果有 非割边 ,那么用它一定比用割边更优,因为不会拆,且能够合并

下面考虑应该和哪个连通块合并,对于两个连通块 ( a 1 , b 1 ) (a_1,b_1) (a1,b1) , ( a 2 , b 2 ) (a_2,b_2) (a2,b2)

a 1 < b 1 , a 2 < b 2 a_1<b_1,a_2<b_2 a1<b1,a2<b2 (不等号方向相同):

贡献变化 m i n ( a 1 + a 2 , b 1 + b 2 ) − m i n ( a 1 , b 1 ) − m i n ( a 2 , b 2 ) = 0 min(a_1+a_2,b_1+b_2)-min(a_1,b_1)-min(a_2,b_2)=0 min(a1+a2,b1+b2)min(a1,b1)min(a2,b2)=0,没影响

a 1 < b 1 , a 2 > b 2 a_1<b_1,a_2>b_2 a1<b1,a2>b2 (相反)

贡献变化 m i n ( a 1 + a 2 , b 1 + b 2 ) − m i n ( a 1 , b 1 ) − m i n ( a 2 , b 2 ) min(a_1+a_2,b_1+b_2)-min(a_1,b_1)-min(a_2,b_2) min(a1+a2,b1+b2)min(a1,b1)min(a2,b2)

= m i n ( a 1 + a 2 , b 1 + b 2 ) − a 1 − b 2 =min(a_1+a_2,b_1+b_2)-a_1-b_2 =min(a1+a2,b1+b2)a1b2
= m i n ( a 2 − b 2 , b 1 − a 1 ) =min(a_2-b_2,b_1-a_1) =min(a2b2,b1a1)
= m i n ( Δ 1 , Δ 2 ) =min(\Delta_1,\Delta_2) =min(Δ1,Δ2)

a 1 , b 1 a_1,b_1 a1,b1 确定时,要找一个不等号相反,且差最大的,直接扫一遍记 Max 就行

如果全是割边,我们直接遍历这棵缩点后的树,依次断开每条树边即可

不太好写

#include<bits/stdc++.h>
using namespace std ;

typedef long long LL ;
const int N = 2e5+100 ; 

int T , K , M1 , M2 ;
int n , a[N] , b[N] ;
struct nn
{
	int lst , to ;
}E[2*N] ;
int head[2*N] , tot = 1 ;
inline void add( int x , int y )
{
	E[++tot] = (nn){ head[x] , y } ;
	head[x] = tot ;
}
int dfn[2*N] , low[2*N] , tim ;
bool bri[2*N] ;
vector<int> e[2*N] ;
void tarjan( int x , int inE )
{
	dfn[x] = low[x] = ++tim ;
	for(int i = head[x] ; i ; i = E[i].lst ) {
		int t = E[i].to ;
		if( !dfn[t] ) {
			tarjan( t , i ) ;
			low[x] = min( low[x] , low[t] ) ;
			if( low[t] > dfn[x] ) {
				bri[i] = bri[i^1] = 1 ;
			}
		}
		else if( i != (inE^1) ) {
			low[x] = min( low[x] , dfn[t] ) ;
		}
	}
}
int cnt , bol , bel[2*N] , A[2*N] , B[2*N] , SA[2*N] , SB[2*N] ;
void dfs( int x )
{
	bel[x] = cnt ;
	if( x <= M1 ) A[cnt] ++ ;
	else B[cnt] ++ ;
	for(int i = head[x] ; i ; i = E[i].lst ) {
		int t = E[i].to ;
		if( bri[i] ) {
			if( bel[t] ) {
				e[bel[t]].push_back( bel[x] ) ;
				e[bel[x]].push_back( bel[t] ) ;
			}
		}
		else if( !bel[t] ) dfs( t ) ;
	}
}
int PA[2*N] , PB[2*N] , ans1 , ans2 , ans3 , rt[2*N] ;
bool can[2*N] ; // 是否有非割边 
void get_size( int x , int fa )
{
	SA[x] = A[x] , SB[x] = B[x] ;
	if( A[x]+B[x] > 1 ) can[bol] = 1 ;
	for(int t : e[x] ) {
		if( t == fa ) continue ;
		get_size( t , x ) ;
		SA[x] += SA[t] ; SB[x] += SB[t] ;
	}
}
void ddfs( int x , int fa )
{
	for(int t : e[x] ) {
		if( t == fa ) continue ;
		ans2 = max( ans2 , min(PA[bol],PB[bol])-min(SA[t],SB[t])-min(PA[bol]-SA[t],PB[bol]-SB[t]) ) ; // 最多减少多少 
		ddfs( t , x ) ;
	}
}
int I , Mx1 , Mx2 ;
inline int V ( int a , int b , int c , int d ) // 合并后增加多少 / 拆开后损失多少 
{
	return min(a+c,b+d)-min(a,b)-min(c,d) ;
}
inline int V1( int a , int b )
{
	if( a > b ) return min(a-b,Mx2) ;
	return min(b-a,Mx1) ;
}
void Get( int x , int fa ) // 处理割边 
{
	for(int t : e[x] ) {
		if( t == fa ) continue ;
		ans3 = max( ans3 , -V(SA[t],SB[t],PA[I]-SA[t],PB[I]-SB[t])+V1(SA[t],SB[t]) ) ;
		ans3 = max( ans3 , -V(SA[t],SB[t],PA[I]-SA[t],PB[I]-SB[t])+V1(PA[I]-SA[t],PB[I]-SB[t]) ) ;
		Get( t , x ) ;
	}
}
void solve()
{
	cnt = 0 ; bol = 0 ;
	for(int i = 1 ; i <= M1+M2 ; i ++ ) {
		if( !dfn[i] ) {
			tarjan( i , 0 ) ;
		}
	}
	for(int i = 1 ; i <= M1+M2 ; i ++ ) {
		if( !bel[i] ) {
			cnt ++ ;
			dfs( i ) ;
		}
	}
	ans1 = ans2 = ans3 = 0 ;
	for(int i = 1 ; i <= cnt ; i ++ ) {
		if( !SA[i]&&!SB[i] ) {
			bol ++ ;
			rt[bol] = i ;
			get_size( i , 0 ) ;
			PA[bol] = SA[i] , PB[bol] = SB[i] ;
			ddfs( i , 0 ) ;
			ans1 += min(PA[bol],PB[bol]) ;
		}
	}
	ans2 = ans1-ans2 ;
	if( K == 0 ) {
		printf("%d\n" , ans1 ) ;
		return ;
	}
	if( K == 1 ) {
		printf("%d\n" , ans2 ) ;
		return ;
	}
	if( K == 2 ) {
		Mx1 = 0 , Mx2 = 0 ;
		if( M1 < n ) Mx1 = 1 ;
		if( M2 < n ) Mx2 = 1 ;
		for(int i = 1 ; i <= bol ; i ++ ) {
			Mx1 = max( Mx1 , PA[i]-PB[i] ) ;
			Mx2 = max( Mx2 , PB[i]-PA[i] ) ;
		}
		for(int i = 1 ; i <= bol ; i ++ ) {
			if( can[i] ) {
				if( PA[i] > PB[i] ) ans3 = max( ans3 , min(PA[i]-PB[i],Mx2) ) ;
				else ans3 = max( ans3 , min(PB[i]-PA[i],Mx1) ) ;
			}
			else {
				I = i ;
				Get( rt[i] , 0 ) ;
			}
		}
		printf("%d\n" , ans3+ans1 ) ;
		return ;
	}
}
int nam[N][2] ;	
void Clear()
{
	for(int i = 2 ; i <= tot ; i ++ ) bri[i] = 0 ;
	tot = 1 ;
	for(int i = 1 ; i <= M1+M2 ; i ++ ) head[i] = dfn[i] = bel[i] = 0 ;
	for(int i = 1 ; i <= n ; i ++ ) nam[i][0] = nam[i][1] = 0 ;
	tim = 0 ; 
	for(int i = 1 ; i <= cnt ; i ++ ) e[i].clear() , A[i] = B[i] = SA[i] = SB[i] = 0 ;
	for(int i = 1 ; i <= bol ; i ++ ) can[i] = 0 ; 
}

int main()
{
	scanf("%d%d" , &T , &K ) ;
	while( T -- ) {
		scanf("%d" , &n ) ;
		if( n == 1 ) {
			printf("1\n") ;
			for(int i = 1 ; i <= n ; i ++ ) scanf("%d" , &a[i] ) ;
			for(int i = 1 ; i <= n ; i ++ ) scanf("%d" , &b[i] ) ;
			continue ;
		}
		M1 = 0 , M2 = 0 ;
		for(int i = 1 ; i <= n ; i ++ ) {
			scanf("%d" , &a[i] ) ;
			if( !nam[a[i]][0] ) {
				nam[a[i]][0] = ++M1 ;
			}
			a[i] = nam[a[i]][0] ;
		}
		for(int i = 1 ; i <= n ; i ++ ) {
			scanf("%d" , &b[i] ) ;
			if( !nam[b[i]][1] ) {
				nam[b[i]][1] = ++M2 ;
			}
			b[i] = nam[b[i]][1] ;
			add( a[i] , b[i]+M1 ) , add( b[i]+M1 , a[i] ) ;
		}
		solve() ;
		Clear() ;
	}
	return 0 ;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值