SDOI 2014 LIS 题解

题目传送门

题目大意： 在一个序列中删去若干个元素，在总代价最小的情况下让删去元素的附加属性排列后字典序最小。

题解

这题很像网络流 $24$ 题中的那个最长不下降子序列，我们模仿它的建图，唯一不同的是，我们将每个元素拆点后，中间连的边的流量为 $b_i$，而其他的边流量为 $inf$，这样可以保证只有点会被割掉，而割掉第 $i$ 个点的代价也就是 $b_i$。

由于需要字典序最小，所以建完图，跑一次最大流后，我们将所有点按照 $c_i$ 排序后一个一个看，如果这个点是割点，那么就将它记录到答案中。然后我们需要抹去这个点在网络中的贡献，重新建图的话会 $TLE$，所以要选择其他的方式。

抹去这个点的贡献，其实就是相当于删掉这个点与它的分身中间的边，那么这个时候用一个小技巧退流即可，设这个点和分身为 $u$ 和 $v$，那么只需要从 $u$ 往 $S$ 跑一次最大流，再从 $T$ 往 $v$ 跑一次最大流即可。

因为对于一条经过 $(u,v)$ 的增广路，上面能割掉的边有很多条，然而我们只需要割掉 $1$ 条，当我们选择割掉 $(u,v)$ 时，就不能割其它的了，这时候我们跑一次退流，把流到 $u,v$ 的流量流回去，这样他们就不会被割掉了。

判断一条边是否是割点很简单，看看跑完最大流之后，还能不能找到 $u$ 到 $v$ 的增广路即可。

题目稍微有一点卡常，稍稍优化一下就能过了：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 2010
#define inf 2147483640
#define ll long long
#define reg register

int t,n,a[maxn],b[maxn];
int f[maxn],max_f;
struct node{int x,y;}c[maxn];
struct edge{int y,z,next;};
edge e[200010];
int first[maxn],len;
void buildroad(reg int x,reg int y,reg int z)
{
	e[++len]=(edge){y,z,first[x]};
	first[x]=len;
}
bool cmp(reg node x,reg node y){return x.x<y.x;}
int q[maxn],st,ed,h[maxn];
bool bfs(reg int S,reg int T,reg int type)
{
	for(int i=1;i<=2*n+2;i++)h[i]=0;
	st=ed=1;q[st]=S;h[S]=1;
	while(st<=ed)
	{
		int x=q[st++];
		for(reg int i=first[x];i;i=e[i].next)
		{
			int y=e[i].y;
			if(!h[y]&&e[i].z)
			{
				if(type&&y==T)return true;
				q[++ed]=y;h[y]=h[x]+1;
			}
		}
	}
	return h[T];
}
int dfs(reg int x,reg int T,reg int flow)
{
	if(x==T)return flow;
	int tt=0;
	for(reg int i=first[x];i;i=e[i].next)
	{
		int y=e[i].y;
		if(h[y]==h[x]+1&&e[i].z)
		{
			int p=dfs(y,T,min(e[i].z,flow-tt));tt+=p;
			e[i].z-=p;e[i^1].z+=p;
			if(tt==flow)break;
		}
	}
	if(!tt)h[x]=0;
	return tt;
}
ll dinic(int S,int T)
{
	ll re=0;
	while(bfs(S,T,0))re+=dfs(S,T,inf);
	return re;
}
int S,T,ans[maxn],tot;

int main()
{
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d",&n);max_f=0;
		for(reg int i=1;i<=n;i++)
		{
			scanf("%d",&a[i]);f[i]=1;
			for(reg int j=1;j<i;j++)if(a[j]<a[i])
			f[i]=max(f[i],f[j]+1),max_f=max(max_f,f[i]);
		}
		for(int i=1;i<=2*n+2;i++)
		first[i]=0;len=1;
		S=2*n+1;T=S+1;
		for(reg int i=1;i<=n;i++)
		{
			if(f[i]==1)buildroad(S,i,inf),buildroad(i,S,0);
			if(f[i]==max_f)buildroad(i+n,T,inf),buildroad(T,i+n,0);
			for(reg int j=i+1;j<=n;j++)
			if(a[i]<a[j]&&f[j]==f[i]+1)
			buildroad(i+n,j,inf),buildroad(j,i+n,0);
		}
		for(reg int i=1;i<=n;i++)scanf("%d",&b[i]),
		buildroad(i,i+n,b[i]),buildroad(i+n,i,0);
		for(reg int i=1;i<=n;i++)scanf("%d",&c[i].x),c[i].y=i;
		sort(c+1,c+n+1,cmp); tot=0;
		printf("%lld ",dinic(S,T));
		for(reg int i=1;i<=n;i++)
		{
			int x=c[i].y;
			if(bfs(x,x+n,1))continue;
			dinic(x,S);dinic(T,x+n);
			ans[++tot]=x;
		}
		sort(ans+1,ans+tot+1);
		printf("%d\n",tot);
		for(int i=1;i<=tot;i++)
		printf("%d ",ans[i]);
		printf("\n");
	}
}