Cogs 731. [网络流24题] 最长递增子序列(最大流)
- [网络流24题] 最长递增子序列
★★★☆ 输入文件:alis.in 输出文件:alis.out 简单对比
时间限制:1 s 内存限制:128 MB
«问题描述:
给定正整数序列x1,…, xn。
(1)计算其最长递增子序列的长度s。
(2)计算从给定的序列中最多可取出多少个长度为s的递增子序列。
(3)如果允许在取出的序列中多次使用x1和xn,则从给定序列中最多可取出多少个长
度为s的递增子序列。
注意:这里的最长递增子序列即最长不下降子序列!!!
«编程任务:
设计有效算法完成(1)(2)(3)提出的计算任务。
«数据输入:
由文件alis.in提供输入数据。文件第1 行有1个正整数n(n<=500),表示给定序列的长度。接
下来的1 行有n个正整数x1,…, xn。
«结果输出:
程序运行结束时,将任务(1)(2)(3)的解答输出到文件alis.out中。第1 行是最长
递增子序列的长度s。第2行是可取出的长度为s 的递增子序列个数。第3行是允许在取出
的序列中多次使用x1和xn时可取出的长度为s 的递增子序列个数。
输入文件示例 输出文件示例
alis.in
4
3 6 2 5
alis.out
2
2
3
/*
好题.
没想出来。。。
第一问求DPLIS,貌似费用流也行2333.
f[i]表示以i为结尾的LIS长度.
第二问和第三问网络流.
因为每个位置只允许用到一次,所以拆点连边.
然后有大小正序关系的连边
(这点想到了,但是没想到是当且仅当f值相差为1的时候.
因为第二问和第三问求的是最长的,
所以f值相差不为1的点显然没必要建边.
第三问的x1,xn疯狂用无非就是该相关边的流量为INF.
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define INF 1e9
#define MAXN 1010
using namespace std;
int n,m,S,T,cut=1,f[MAXN],a[MAXN],head[MAXN],dis[MAXN],len,ans;
struct edge{int v,next,c;}e[MAXN*MAXN*2];
queue<int>q;
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9') {if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
return x*f;
}
void add(int u,int v,int c)
{
e[++cut].v=v;e[cut].c=c;e[cut].next=head[u];head[u]=cut;
e[++cut].v=u;e[cut].c=0;e[cut].next=head[v];head[v]=cut;
}
void LIS()
{
for(int i=1;i<=n;i++)
{
f[i]=1;
for(int j=1;j<i;j++)
if(a[i]>=a[j]) f[i]=max(f[i],f[j]+1);
len=max(len,f[i]);
}
printf("%d\n",len);
}
bool bfs()
{
for(int i=S;i<=T;i++) dis[i]=-1;dis[S]=0;q.push(S);
while(!q.empty())
{
int u=q.front();q.pop();
for(int i=head[u];i;i=e[i].next)
{
int v=e[i].v;
if(dis[v]==-1&&e[i].c) dis[v]=dis[u]+1,q.push(v);
}
}
return dis[T]!=-1;
}
int dfs(int u,int y)
{
if(u==T) return y;
int rest=0;
for(int i=head[u];i&&rest<y;i=e[i].next)
{
int v=e[i].v;
if(dis[v]==dis[u]+1&&e[i].c)
{
int x=dfs(v,min(y-rest,e[i].c));
rest+=x;
e[i].c-=x;
e[i^1].c+=x;
}
}
if(!rest) dis[u]=-1;
return rest;
}
void dinic()
{
while(bfs()) ans+=dfs(S,INF);
printf("%d\n",ans);
}
void slove()
{
S=0,T=2*n+1;
for(int i=1;i<=n;i++) add(i,i+n,1);
for(int i=1;i<=n;i++) if(f[i]==1) add(S,i,1);
for(int i=1;i<=n;i++) if(f[i]==len) add(i+n,T,1);
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++)
{
if(a[j]>=a[i]&&f[j]==f[i]+1) add(i+n,j,1);
}
dinic();
}
void slove2()
{
memset(head,0,sizeof head);
S=0,T=2*n+1;cut=1;ans=0;add(1,n+1,INF),add(n,2*n,INF);;
if(f[1]==1) add(S,1,INF);
if(f[n]==len) add(2*n,T,INF);
for(int i=2;i<=n-1;i++) add(i,i+n,1);
for(int i=2;i<=n-1;i++) if(f[i]==1) add(S,i,1);
for(int i=2;i<=n-1;i++) if(f[i]==len) add(i+n,T,1);
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++)
{
if(a[j]>=a[i]&&f[j]==f[i]+1) add(i+n,j,1);
}
dinic();
}
int main()
{
freopen("alis.in","r",stdin);
freopen("alis.out","w",stdout);
n=read();
for(int i=1;i<=n;i++) a[i]=read();
LIS();
if(len==1){
printf("%d\n%d",n,n);return 0;
}
slove();
slove2();
return 0;
}
Cogs 731. [网络流24题] 最长递增子序列(最大流)的更多相关文章
- COGS731 [网络流24题] 最长递增子序列(最大流)
给定正整数序列x1,..., xn (n<=500).(1)计算其最长递增子序列的长度s.(2)计算从给定的序列中最多可取出多少个长度为s的递增子序列.(3)如果允许在取出的序列中多次使用x1和 ...
- loj6005 [网络流24题]最长递增子序列
题意:给你一个序列,求不严格上升lcs长度/最多有几个没有重复元素的lcs/如果x1和xn可以多次出现,求最多有几个lcs?n<=500. 标程: #include<cstdio> ...
- COGS743. [网络流24题] 最长k可重区间集
743. [网络流24题] 最长k可重区间集 ★★★ 输入文件:interv.in 输出文件:interv.out 简单对比时间限制:1 s 内存限制:128 MB «问题描述: «编 ...
- [网络流24题]最长k可重区间集[题解]
最长 \(k\) 可重区间集 题目大意 给定实心直线 \(L\) 上 \(n\) 个开区间组成的集合 \(I\) ,和一个正整数 \(k\) ,试设计一个算法,从开区间集合 \(I\) 中选取开区间集 ...
- Cogs 14. [网络流24题] 搭配飞行员
这道题其实蛮好想的,因为分为正,副飞行员.所以就把正飞行员当作Boy,副飞行员当作Girl.然后做Hungry即可. #include<bits/stdc++.h> using names ...
- [COGS 0014][网络流24题] 搭配飞行员
先贴题面 14. [网络流24题] 搭配飞行员 ★★☆ 输入文件:flyer.in 输出文件:flyer.out简单对比时间限制:1 s 内存限制:128 MB [问题描述] 飞行 ...
- Cogs 732. [网络流24题] 试题库(二分图)
[网络流24题] 试题库 ★★ 输入文件:testlib.in 输出文件:testlib.out 评测插件 时间限制:1 s 内存限制:128 MB «问题描述: 假设一个试题库中有n道试题.每道试题 ...
- Cogs 739. [网络流24题] 运输问题(费用流)
[网络流24题] 运输问题 ★★ 输入文件:tran.in 输出文件:tran.out 简单对比 时间限制:1 s 内存限制:128 MB «问题描述: «编程任务: 对于给定的m 个仓库和n 个零售 ...
- Cogs 727. [网络流24题] 太空飞行计划(最大权闭合子图)
[网络流24题] 太空飞行计划 ★★☆ 输入文件:shuttle.in 输出文件:shuttle.out 简单对比 时间限制:1 s 内存限制:128 MB [问题描述] W 教授正在为国家航天中心计 ...
随机推荐
- 音量调节条-封装通用的ProgressBar组件
import React, { Component } from 'react' import PropTypes from 'prop-types' import assign from 'obje ...
- poj 3252 数位dp
题意:一个二进制的数,如果0的个数大于1的个数,那么我们称这个数为Round Numbers,求给定区间(十进制表示)中Round Numbers的个数 题解:数位dp,不过这里枚举的时候lead标记 ...
- RuntimeError: Model class users.models.UserProfile doesn't declare an explicit app_label and isn't in an application in INSTALLED_APPS.
Django启动的时候报错 File "/home/hehecat/PycharmProjects/MxShop/MxShop/urls.py", line 23, in from ...
- module 'cv2' has no attribute 'KNearest_create'
python版本:3.6.5 opencv版本:3.2.0 使用的是jupyter notebook 源代码如下: import cv2 import numpy as np import matpl ...
- Python基础数据类型(数字、字符串、布尔、列表、元组、字典、集合等)
数据类型 计算机顾名思义就是可以做数学计算的机器,因此,计算机程序理所当然地可以处理各种数值.但是,计算机能处理的远不止数值,还可以处理文本.图形.音频.视频.网页等各种各样的数据,不同的数据,需要定 ...
- css 盒子 取值
盒子:当我们设置一个标签宽高时,默认设置的是盒子里面content大小. 内容盒:content 填充盒:content+padding(overflow截取的区域) 边框盒:content+padd ...
- HTML5页面介绍
1.<!DOCTYPE html> 文档声明:用于告诉浏览器使用html哪个版本的标准解析页面,此写法代表使用html5的标准去解析 2.<html> 根标签, ...
- $.ajaxSetup()与$.ajax()区别
ajaxSetup()其使用方法与ajax()并无区别,只是在其之后的ajax()不用再重复添加相同参数,节省了代码量. 附:API中关于jQuery.ajaxSetup([options])的描述和 ...
- 100行代码打造属于自己的代理ip池
经常使用爬虫的朋友对代理ip应该比较熟悉,代理ip就是可以模拟一个ip地址去访问某个网站.我们有时候需要爬取某个网站的大量信息时,可能由于我们爬的次数太多导致我们的ip被对方的服务器暂时屏蔽(也就是所 ...
- Ubuntu18.04开机启动sudo命令
首先接前文:ubuntu18.04 下启动Android Studio报错KVM is required to run this AVD. /dev/kvm device: permission de ...