关于模拟网络流
题意
其实是CF 724 E
这其实是师兄改编了题目QMQ,真实范围应该是\(n<=10000\)
思路
这道题目很明显可以用网络流来做:
但是范围直接T了。
然后我们可以用贪心
然后我们发现这个图貌似有神奇的性质,利用最大流=最小割,我们可以枚举最小割。
最小割的性质就是每个点要么位于\(st\)集合或者\(ed\)集合,所以我们就可以先枚举出一种情况:
我们发现这个中间的边也要处理,所以对于\(st,ed\)我们还要处理一下,但是我们要如何处理最小割呢?
我们设\(f[i][j]\)表示的是到了第\(i\)个点,有\(j\)个归到\(st\)集合。
那么不计\(C\)的影响的话,不难列出:
\(f[i][j]=min(f[i-1][j-1]+a[i],f[i-1][j]+b[i])\)。
但是如何处理\(C\)呢,不难发现对于\(i,j(i<j)\)如果\(i\)选了\(ed\),\(j\)选了\(st\),那么就会有个\(C\)。
这个该归到\(st\)还是\(ed\)呢,很明显为了DP没有后效性,我们就归到\(ed\)吧。
那么DP转移方程变成了:\(f[i][j]=min(f[i-1][j-1]+a[i],f[i-1][j]+b[i]+j*C)\)。
就可以转移了,时间复杂度\(O(n^2)\),当然,这个又叫模拟网络流。
我的思路
我是采用贪心,这道题目\(n\)才\(2000\),还开\(2s\),那么我们不难想到一种贪心思路,就是使得所有卖完粮食后剩余粮食大于\(C\)的位置的粮食尽量的平衡。
也就是尽可能的榨干这个\(C\),使得后面每个位置都可以堆满\(C\)。
至于平衡的值,我们采用二分查找,那么时间复杂度\(O(n^2log值域)\)。
而且还有优化,\(1s\)内过掉。
#include<cstdio>
#include<cstring>
#include<queue>
#define N 3100
using namespace std;
typedef long long LL;
template <class T>
inline T mymin(T x,T y){return x<y?x:y;}
int a[N],b[N];
priority_queue<int>q;//储存粮食信息,当然是大于C的才给放
int sta[N],top;//表示可以移动的粮食
int list[N],tail;
int n,m;
LL zans;
bool check(int x,int k/*原本的数字*/)
{
for(int i=1;i<=top;i++)
{
if(sta[i]<=x)return false;
k+=mymin(sta[i]-x,m);//可以给你多少的数字
if(k>=x)return true;
}
return false;
}
void work(int x,int y)
{
top=0;
while(!q.empty())sta[++top]=q.top(),q.pop();
int k=x-y;//可能是负数
if(top)//有大于C的粮食位置
{
int l=k+1,r=mymin((LL)sta[1],(LL)k+(LL)top*(LL)m),mid,ans=k/*就是不变*/;//表示范围,而且能防止m过小时时间过大
while(l<=r)
{
mid=(l+r)/2;
if(check(mid,k)==true)ans=mid,l=mid+1;
else r=mid-1;
}
for(int i=1;i<=top;i++)
{
if(k==ans)break;
int zjj=mymin(mymin(ans-k,m),sta[i]-ans);sta[i]-=zjj;k+=zjj;
}
}
if(k<0)
{
while(k<0 && tail)k+=list[tail--];
}
zans+=mymin(y,k+y);
if(k>0)
{
if(k<=m)list[++tail]=k;
else q.push(k);
}
for(int i=1;i<=top;i++)
{
if(sta[i]<=m)list[++tail]=sta[i];
else q.push(sta[i]);
}
}
int main()
{
// freopen("c.in","r",stdin);
// freopen("c.out","w",stdout);
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%d",&a[i]);
for(int i=1;i<=n;i++)scanf("%d",&b[i]);
for(int i=1;i<=n;i++)
{
work(a[i],b[i]);
}
printf("%lld\n",zans);
return 0;
}
原文地址:https://www.cnblogs.com/zhangjianjunab/p/11789135.html
- 洛谷P3195 [HNOI2008]玩具装箱TOY(单调队列优化DP)
- SQL Server 深入解析索引存储(下)
- 2751: [HAOI2012]容易题(easy)
- codevs3002 石子归并 3
- 算法模板——计算几何2(二维凸包——Andrew算法)
- 算法模板——splay区间反转 2
- 算法模板——Dinic网络最大流 2
- 1935: [Shoi2007]Tree 园丁的烦恼
- 1339 / 1163: [Baltic2008]Mafia
- 4010: [HNOI2015]菜肴制作
- 4052: [Cerc2013]Magical GCD
- ElasticSearch搜索引擎在SpringBoot中的实践
- 2292: 【POJ Challenge 】永远挑战
- 四边形不等式优化DP
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- [WPF] WPF做的漂亮的登陆界面[附源码]
- C#自定义控件的创建
- Oracle参数解析(processor_group_name)
- [Oracle故障处理]记一次PX msg pool 4031错误的处理
- WPF的布局-Grid(表格布局)
- C# 通过遍历设置控件属性
- C# 软件开机启动
- Java之映射
- C# 武汉肺炎全国疫情实时信息图
- brew报错:`initialize': Version value must be a string; got a NilClass () (TypeError)
- gitlab和gitlab项目迁移
- Qt5+VS2017点击按钮打开新窗口
- C#连接MySQL数据库
- 了不起的 IoC 与 DI
- QT之UDP通信