正题

题目链接:https://www.luogu.com.cn/problem/P3308


题目大意

三个\(n\)个数字的序列\(A,B,C\)。要求删除其中某些位置\(i\)使得\(A\)的最长上升子序列至少减少\(1\)且删去位置\(B\)的权值和最小的情况下满足删去位置的\(C\)值升序排序后字典序最小。


解题思路

首先\(B\)值最小很好求,跑一遍\(LIS\)的\(dp\),然后每个点拆成两个点,然后如果\(f[x]\)转移到\(f[y]\)是最优的就建边然后跑最小割就好了。

大体和P2766 最长不下降子序列问题差不多

也就是现在我们要求字典序最小的最小割,需要利用到最小割的性质。

如果一条边\(x,y\)是可行割,那么它满足在残量网络上\(x\)到达不了\(y\)。

首先如果\(x->y\)没有满流那么肯定不是最小割,其次如果满流了但是还有一条\(x\)到\(y\)的路径,那么证明如果走这条增广路一定可以使最大流更大,所以也不是最小割。

那么这样我们就可以判断一条边是否可行了,然后需要消去其他等价边的影响,大体方法是从\(T\)到\(y\)跑一次\(dinic\),再从\(x\)到\(S\)跑一次\(dinic\)。这个操作叫退流,这样残量网络就变成了满流边\(x->y\)的残量网络了。

先跑一次\(dinic\),然后按照\(C\)值排序,从小到大判断加入边即可。


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#define ll long long
using namespace std;
const ll N=710*2,inf=1e18;
struct node{
ll to,next,w;
}a[N*N];
ll T,n,tot,ls[N],dep[N],A[N],B[N],C[N],f[N],p[N];
vector<ll> prt;queue<ll> q;
void addl(ll x,ll y,ll w){
a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;
a[++tot].to=x;a[tot].next=ls[y];ls[y]=tot;a[tot].w=0;
}
bool bfs(ll s,ll t){
while(!q.empty())q.pop();q.push(s);
memset(dep,0,sizeof(dep));dep[s]=1;
while(!q.empty()){
ll x=q.front();q.pop();
for(ll i=ls[x];i;i=a[i].next){
ll y=a[i].to;
if(dep[y]||!a[i].w)continue;
dep[y]=dep[x]+1;
if(y==t)return 1;
q.push(y);
}
}
return 0;
}
ll dinic(ll x,ll t,ll flow){
if(x==t)return flow;
ll rest=0,k;
for(ll i=ls[x];i;i=a[i].next){
ll y=a[i].to;
if(dep[x]+1!=dep[y]||!a[i].w)continue;
rest+=(k=dinic(y,t,min(flow-rest,a[i].w)));
a[i].w-=k;a[i^1].w+=k;
if(rest==flow)return flow;
}
if(!rest)dep[x]=0;
return rest;
}
bool cmp(ll x,ll y)
{return C[x]<C[y];}
signed main()
{
scanf("%lld",&T);
while(T--){
memset(ls,0,sizeof(ls));tot=1;
scanf("%lld",&n);
for(ll i=1;i<=n;i++)scanf("%lld",&A[i]);
for(ll i=1;i<=n;i++)scanf("%lld",&B[i]);
for(ll i=1;i<=n;i++)scanf("%lld",&C[i]);
ll maxs=0,s=2*n+1,t=s+1;
for(ll i=1;i<=n;i++){
f[i]=1;p[i]=i;
for(ll j=1;j<i;j++)
if(A[i]>A[j])f[i]=max(f[i],f[j]+1);
maxs=max(maxs,f[i]);
}
for(ll i=1;i<=n;i++){
if(f[i]==1)addl(s,i,inf);
if(f[i]==maxs)addl(i+n,t,inf);
addl(i,i+n,B[i]);
for(ll j=i+1;j<=n;j++)
if(A[i]<A[j]&&f[i]+1==f[j])
addl(i+n,j,inf);
}
ll ans=0;prt.clear();
while(bfs(s,t))
ans+=dinic(s,t,inf);
printf("%lld ",ans);
sort(p+1,p+1+n,cmp);
for(ll i=1;i<=n;i++){
ll x=p[i];
if(bfs(x,x+n))continue;
while(bfs(t,x+n))dinic(t,x+n,inf);
while(bfs(x,s))dinic(x,s,inf);
prt.push_back(x);
}
printf("%lld\n",prt.size());
sort(prt.begin(),prt.end());
for(ll i=0;i<prt.size();i++)
printf("%lld ",prt[i]);
putchar('\n');
}
return 0;
}

P3308-[SDOI2014]LIS【最小割】的更多相关文章

  1. BZOJ.3532.[SDOI2014]LIS(最小割ISAP 退流)

    BZOJ 洛谷 \(LIS\)..经典模型? 令\(f_i\)表示以\(i\)结尾的\(LIS\)长度. 如果\(f_i=1\),连边\((S,i,INF)\):如果\(f_i=\max\limits ...

  2. P3308 [SDOI2014]LIS(最小割+退流)

    传送门 设\(f[i]\)为以\(i\)结尾的最长上升子序列.可以考虑建这样一张图,对于所有的\(i<j,f[j]=f[i+1]\)连边\((i,j)\),\(f[i]=1\)的话连边\((S, ...

  3. 3532: [Sdoi2014]Lis 最小字典序最小割

    3532: [Sdoi2014]Lis Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 865  Solved: 311[Submit][Status] ...

  4. 【BZOJ-3532】Lis 最小割 + 退流

    3532: [Sdoi2014]Lis Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 704  Solved: 264[Submit][Status] ...

  5. hdu3739 Anti LIS[最小割]

    长度为 n≤1000 的数列 ai,其中最长上升子序列的长度为 s.至少删去多少数使得最长上升子序列的长度小于 s. 其实这题和那个求有多少不重叠LIS是一样答案的. 先放个图. 图丑别说我. 原网络 ...

  6. 洛谷$P3308\ [SDOI2014]LIS$ 网络流

    正解:网络流 解题报告: 传送门$QwQ$ 恩先不考虑关于那个附加属性的限制,考虑这题怎么做? 首先这题从名字开始就让人忍不住联想起网络流24题里的那个最长不下降子序列?于是同样考虑预处理一个$f$呗 ...

  7. [bzoj3532][Sdoi2014]Lis——拆点最小割+字典序+退流

    题目大意 给定序列A,序列中的每一项Ai有删除代价Bi和附加属性Ci.请删除若 干项,使得4的最长上升子序列长度减少至少1,且付出的代价之和最小,并输出方案. 如果有多种方案,请输出将删去项的附加属性 ...

  8. bzoj千题计划141:bzoj3532: [Sdoi2014]Lis

    http://www.lydsy.com/JudgeOnline/problem.php?id=3532 如果没有字典序的限制,那么DP拆点最小割即可 加上字典序的限制: 按c从小到大枚举最小割边集中 ...

  9. BZOJ3532 : [Sdoi2014]Lis

    f[i]表示以i为结尾的LIS长度 对于所有f[i]=1的,由S向i连边 对于所有f[i]=maxf的,由i向T连边 对于j<i,a[j]<a[i],且f[j]+1=f[i]的,j向i连边 ...

随机推荐

  1. Clusternet - 新一代开源多集群管理与应用治理项目

    作者 徐迪,腾讯云容器技术专家. 汝英哲,腾讯云高级产品经理. 摘要 在过去的数年里,云计算领域经历了多次巨大的变革,当前越来越多的组织将应用部署在本地和云上的多个基础设施平台上,这些平台可能是两个公 ...

  2. 转: SIFT原理解释

    1.SIFT概述 SIFT的全称是Scale Invariant Feature Transform,尺度不变特征变换,由加拿大教授David G.Lowe提出的.SIFT特征对旋转.尺度缩放.亮度变 ...

  3. 【spring 注解驱动开发】扩展原理

    尚学堂spring 注解驱动开发学习笔记之 - 扩展原理 扩展原理 1.扩展原理-BeanFactoryPostProcessor BeanFactoryPostProcessor * 扩展原理: * ...

  4. 等待唤醒机制----线程池----lambda表达式

    1.等待唤醒机制 1.1线程间通信 概念:多个线程在处理同一个资源,但是处理的动作(线程的任务)却不相同. 比如:线程A用来生成包子的,线程B用来吃包子的,包子可以理解为同一资源,线程A与线程B处理的 ...

  5. ProjectEuler 006题

    题目: The sum of the squares of the first ten natural numbers is, 12 + 22 + ... + 102 = 385 The square ...

  6. mybatis根据表逆向自动化生成代码(自动生成实体类、mapper文件、mapper.xml文件)

    .personSunflowerP { background: rgba(51, 153, 0, 0.66); border-bottom: 1px solid rgba(0, 102, 0, 1); ...

  7. CNN的Pytorch实现(LeNet)

    CNN的Pytorch实现(LeNet)   上次写了一篇CNN的详解,可是累坏了老僧我.写完后拿给朋友看,朋友说你这Pytorch的实现方式对于新人来讲会很不友好,然后反问我说里面所有的细节你都明白 ...

  8. java和js中for循环的区别

    java中for循环,先执行语句后循环 for (int i=1;i<10;i++){ for(int b=1;b<=i;b++){ System.out.print(b+"*& ...

  9. ks.cfg文件相关

    原文转自:https://www.cnblogs.com/itzgr/p/10029631.html作者:木二 目录 一 图形化生成ks.cfg文件 二 ks.cfg文件相关项解析 一 图形化生成ks ...

  10. Linux centos7 复制,移动,删除文件或文件夹

    2021-08-121. 文件(文件夹)复制命令 # 命令格式 cp [-adfilprsu] 源文件(source) 目标文件(destination) cp [option] source1 so ...