CF739E Gosha is hunting(费用流/凸优化dp)
纪念合格考爆炸。
其实这个题之前就写过博客了,qwq但是不小心弄丢了,所以今天来补一下。
首先,一看到球的个数的限制,不难相当用网络流的流量来限制每个球使用的数量。
由于涉及到最大化期望,所以要使用最大费用最大流。
我们新建两个点\(ss,tt\),分别表示两种球。
那么我们现在考虑应该怎么计算期望呢。
首先,如果假设如果对于一个怪物用一个球,那么连边也就比较容易了
对于一个怪物\(x\)
我们\(ss -> x\),费用为\(p[i]\),流量为1。表示一个球在一个怪物上只能用一次。
\(tt\)也是同理。
然后对于每一个\(x->t\),费用是\(0\),流量是\(1\),表示一个怪物只能用一个球。
但是,要是每次不要求只能用一个球应该怎么做呢。
我们考虑,这条边的费用应该是多少。
两个球都用的期望应该是\(1-(1-p_i)(1-q_i)\)
经过展开,我们发现应该是\(p_i+q_i-p_i\times q_i\)
那么由于我们发现,由于用了两个球,所以已经获得了二者之和的收益,那么在这一侧,只需要在上述建图的基础上\(x->t\),费用是\(-p_i\times q_i\)即可。
最后跑一发最大费用最大流就能通过这个题qwq时间复杂度玄学。
#include<bits/stdc++.h>
#define pb push_back
#define mk make_pair
#define ll long long
#define db double
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
const int maxn = 4010;
const int maxm = 3e6+1e2;
const double eps = 1e-10;
int point[maxn],nxt[maxm],to[maxm],pre[maxm],from[maxn];
double dis[maxn];
int vis[maxn];
double cost[maxm];
int flow[maxm];
double ans;
int n,m,cnt=1;
int s,t;
void addedge(int x,int y,db w,int f)
{
nxt[++cnt]=point[x];
pre[cnt]=x;
to[cnt]=y;
cost[cnt]=w;
flow[cnt]=f;
point[x]=cnt;
}
void insert(int x,int y,db w,int f)
{
addedge(x,y,w,f);
addedge(y,x,-w,0);
}
queue<int> q;
bool spfa(int s)
{
for (int i=1;i<=maxn-3;i++) dis[i]=-1e9;
memset(vis,0,sizeof(vis));
q.push(s);
dis[s]=0;
while (!q.empty())
{
int x = q.front();
q.pop();
vis[x]=0;
for (int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (dis[p]-(dis[x]+cost[i])<-eps && flow[i]>0)
{
from[p]=i;
dis[p]=dis[x]+cost[i];
if (!vis[p])
{
q.push(p);
vis[p]=1;
}
}
}
}
if (dis[t]==-1e9) return false;
return true;
}
void mcf()
{
double x = 1e9;
for (int i=from[t];i;i=from[pre[i]]) x=min(x,1.0*flow[i]);
for (int i=from[t];i;i=from[pre[i]])
{
flow[i]-=x;
flow[i^1]+=x;
ans+=x*cost[i];
}
}
void solve()
{
while (spfa(s)) mcf();
}
db a[maxn],b[maxn];
int ss,sss;
int aa,bb;
int main()
{
n=read(),aa=read(),bb=read();
s=maxn-10;
ss=s+1;
t=s+3;
sss=ss+1;
insert(s,ss,0,aa);
insert(s,sss,0,bb);
for (int i=1;i<=n;i++) scanf("%lf",&a[i]);
for (int i=1;i<=n;i++) scanf("%lf",&b[i]);
for (int i=1;i<=n;i++)
{
insert(ss,i,a[i],1);
insert(sss,i,b[i],1);
insert(i,t,0,1);
insert(i,t,-a[i]*b[i],1);
}
solve();
printf("%.4lf\n",ans);
return 0;
}
但是其实这个题的正解是凸优化\(dp\)
首先,先做一个最\(naive\)的想法。
我们令\(dp[i][j][k]\)表示前\(i\)个怪物,用了\(j\)一号球,用了\(k\)个二号球
那么转移也是显然的。
每次只需要讨论一下对于当前的怪物是用几个球,用哪个即可。
但是这样的复杂度是\(O(n^3)\)的。
显然没有办法通过。
考虑怎么优化。
由于题目中涉及到的正好用几个球,并且通过打表发现函数是凸的,那么我们就可以直接用凸优化来优化掉一维。
(其实是可以直接优化两个的,但是我太懒,所以没写。)
我们对于当前二分的值,表示每选一个二号球,就可以多得到\(mid\)的期望。不限制选的个数。
那么不难得到下面的这个转移式子。
dp[i][j]=dp[i-1][j];
dp[i][j]=max(dp[i][j],dp[i-1][j]+(ymh){bb[i],1});
if (j)
{
dp[i][j]=max(dp[i][j],dp[i-1][j-1]+(ymh){a[i],0});
dp[i][j]=max(dp[i][j],dp[i-1][j-1]+(ymh){both[i],1});
}
然后通过调整\(mid\),通过正好选到\(k\)个二号球。
最后求一个\(dp\)数组,然后记得把贡献减去就行。
时间复杂度\(n^2log\),非常优秀。
(其实是如果精度太小会\(WA\),精度太大会\(TLE\))
但是完全可以做到\(nlog^2\)的。
给代码。
#include<bits/stdc++.h>
#define pb push_back
#define mk make_pair
#define ll long long
#define db double
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
const int maxn = 2010;
const db eps = 1e-6;
struct ymh{
db val;
int num;
ymh operator + (const ymh &b) const
{
return (ymh){val+b.val,num+b.num};
}
};
ymh dp[maxn][maxn];
db a[maxn],b[maxn];
int n;
db l=-4,r=4;
inline int dcmp(double x,double y)
{
return x-y<-eps ? -1 : (x-y>eps ? 1 : 0);
}
inline ymh max(ymh a,ymh b)
{
int now = dcmp(a.val,b.val);
if (now==0)
{
if (a.num<b.num) return a;
else return b;
}
else
{
if(now==-1) return b;
else return a;
}
}
int numa,numb;
db aa[maxn];
db bb[maxn];
db both[maxn];
bool check(db lim)
{
for (int i=1;i<=n;i++) aa[i]=a[i];
for (int i=1;i<=n;i++) bb[i]=b[i]+lim;
for (int i=1;i<=n;i++) both[i]=1.0-(1.0-a[i])*(1.0-b[i])+lim;
for (register int i=1;i<=n;++i)
{
for (register int j=0;j<=numa;++j)
{
dp[i][j]=dp[i-1][j];
dp[i][j]=max(dp[i][j],dp[i-1][j]+(ymh){bb[i],1});
if (j)
{
dp[i][j]=max(dp[i][j],dp[i-1][j-1]+(ymh){a[i],0});
dp[i][j]=max(dp[i][j],dp[i-1][j-1]+(ymh){both[i],1});
}
// if (dp[i][j].num>numb) return false;
}
}
return dp[n][numa].num<=numb;
}
int main()
{
n=read(),numa=read(),numb=read();
for (int i=1;i<=n;i++) scanf("%lf",&a[i]);
for (int i=1;i<=n;i++) scanf("%lf",&b[i]);
double ans=0;
while (r-l>=eps)
{
db mid = (l+r)/2;
// memset(dp,0,sizeof(dp));
if (check(mid)) l=mid,ans=mid;
else r=mid;
//printf("%.4lf %d\n",mid,dp[n][numa].num);
}
//cout<<1<<endl;
//printf("%.4lf\n",ans);
//memset(dp,0,sizeof(dp));
check(ans);
//printf("")
//printf("%.4lf %d\n",dp[n][numa].val,dp[n][numa].num);
printf("%.4lf",dp[n][numa].val-1.0*numb*ans);
return 0;
}
CF739E Gosha is hunting(费用流/凸优化dp)的更多相关文章
- CF739E Gosha is hunting(费用流,期望)
根据期望的线性性答案就是捕捉每一只精灵的概率之和. 捕捉一只精灵的方案如下: 1.使用一个\(A\)精灵球,贡献为\(A[i]\) 2.使用一个\(B\)精灵球,贡献为\(B[i]\) 3.使用一个\ ...
- HZOJ 赤(CF739E Gosha is hunting)
本来没有打算写题解的,时间有点紧.但是这个wqs二分看了好久才明白还是写点东西吧. 题解就直接粘dg的了: 赤(red) 本题来自codeforces 739E,加大了数据范围. 首先对一只猫不会扔两 ...
- hdu 2686 费用流 / 双线程DP
题意:给一个方阵,求从左上角出到右下角(并返回到起点),经过每个点一次不重复,求最大获益(走到某处获得改点数值),下来时每次只能向右或向下,反之向上或向左. 俩种解法: 1 费用流法:思路转化:从左 ...
- CF739E Gosha is hunting
法一: 匹配问题,网络流! 最大费用最大流,S到A,B流a/b费0,A,B到i流1费p[i]/u[i],同时选择再减p[i]*u[i]? 连二次!所以i到T流1费0流1费-p[i]*u[i] 最大流由 ...
- CF739E Gosha is hunting 【WQS二分 + 期望】
题目链接 CF739E 题解 抓住个数的期望即为概率之和 使用\(A\)的期望为\(p[i]\) 使用\(B\)的期望为\(u[i]\) 都使用的期望为\(p[i] + u[i] - u[i]p[i] ...
- CF739E Gosha is hunting DP+wqs二分
我是从其他博客里看到这题的,上面说做法是wqs二分套wqs二分?但是我好懒呀,只用了一个wqs二分,于是\(O(nlog^2n)\)→\(O(n^2logn)\) 首先我们有一个\(O(n^3)\)的 ...
- BZOJ2040[2009国家集训队]拯救Protoss的故乡——模拟费用流+线段树+树链剖分
题目描述 在星历2012年,星灵英雄Zeratul预测到他所在的Aiur行星在M天后会发生持续性暴雨灾害,尤其是他们的首都.而Zeratul作为星灵族的英雄,当然是要尽自己最大的努力帮助星灵族渡过这场 ...
- 【BZOJ5252】林克卡特树(动态规划,凸优化)
[BZOJ5252]林克卡特树(动态规划,凸优化) 题面 BZOJ(交不了) 洛谷 题解 这个东西显然是随着断开的越来越多,收益增长速度渐渐放慢. 所以可以凸优化. 考虑一个和\(k\)相关的\(dp ...
- 【CF739E】Gosha is hunting(动态规划,凸优化)
[CF739E]Gosha is hunting(动态规划,凸优化) 题面 洛谷 CF 题解 一个\(O(n^3)\)的\(dp\)很容易写出来. 我们设\(f[i][a][b]\)表示前\(i\)个 ...
随机推荐
- validity属性返回对象中的属性值
- 微信小程序--聊天室小程序(云开发)
微信小程序 -- 聊天室小程序(云开发) 从微信小程序开发社区更新watch接口之后,一直在构思这个项目.项目已经完成很久,但是一直都没有空写一篇博客记录展示一下. 开源地址 wx-cloud-im: ...
- openresty lua-resty-string md5 sha aes random string
安装 https://github.com/openresty/lua-resty-string $ sudo opm get openresty/lua-resty-string $ ls -al ...
- MyBatis学习总结(五)——关联表查询的实现
一.一对一关联 1.1.提出需求 根据班级id查询班级信息(带老师的信息) 1.2.创建表和数据 创建一张教师表和班级表,这里我们假设一个老师只负责教一个班,那么老师和班级之间的关系就是一种一对一的关 ...
- Struts2 的 OGNL
ONGL简介: OGNL 的全称是对象图导航语言( Object-Graph Navigation Language),它是一种功能强大的开源表达式语言,使用这种表达式语言,可以通过某种表达式语法,存 ...
- 浅谈 JVM GC 的安全点与安全区域
OopMap 前文我们说到,JVM 采用的可达性分析法有个缺点,就是从 GC Roots 找引用链耗时. 都说他耗时,他究竟耗时在哪里? GC 进行扫描时,需要查看每个位置存储的是不是引用类型,如果是 ...
- web安全性测试——XSS跨站攻击
1.跨站攻击含义 XSS:(Cross-site scripting)全称"跨站脚本",是注入攻击的一种.其特点是不对服务器端造成任何伤害,而是通过一些正常的站内交互途径,例如发布 ...
- 基于Ubuntu18.04一站式部署(python-mysql-redis-nginx)
基于Ubuntu18.04一站式部署 Python3.6.8的安装 1. 安装依赖 ~$ sudo apt install openssl* zlib* 2. 安装python3.6.8(个人建议从官 ...
- JDK 1.7 正式发布,Oracle 官宣免费提供!“新版任你发,我用JDK 8”或成历史?
Oracle公司JDK 17正式发布,JDK 17属于长期支持(LTS)版本,也就是获得8年的技术支持,自2021年9月至2029年9月截止. JDK 17版本更新了很多比较实用的新特性,关于此版本的 ...
- python+echarts+flask实现对全国疫情数据的爬取并可视化展示
用Python进行数据爬取并存储到数据库,3.15学习总结(Python爬取网站数据并存入数据库) - 天岁 - 博客园 (cnblogs.com) 通过echarts+flask实现数据的可视化展示 ...