洛谷P6054:开门大吉

Description\text{Description}Description

P6054 开门大吉

nnn 位选手去参加节目“开门大吉”。共有 mmm 套题,每套题包含 ppp 个题目,第 iii 位选手答对第 jjj 套题中第 kkk 道的概率为 fi,j,kf_{i,j,k}fi,j,k
若一位选手答对第 iii 题,会在已得到奖励的基础上,再得到 cic_ici 元奖励。选手总是从第一道开始,按顺序答题的。
同时,为了防止过多的选手做同一套题,还有 yyy 条形如“第 iii 位选手的套题编号必须至少比第 jjj 位的大 kkk”的限制。
你需要给每一位选手分配一套题(不同选手可以相同),使得所有人的期望奖励和最小。

Solution\text{Solution}Solution

首先,我们可以简单的算出 iii 选手选套题 jjj 的期望价值 vi,jv_{i,j}vi,j
关键就在于如何表示限制。
建出 nnn 条长度为 m+1m+1m+1 的链,设链上的点为 p1...n,1..m+1p_{1...n,1..m+1}p1...n,1..m+1,连边 (S,pi,1,INF),(pi,m+1,T,INF),(pi,j,pi,j+1,vi,j)(S,p_{i,1},INF),(p_{i,m+1},T,INF),(p_{i,j},p_{i,j+1},v_{i,j})(S,pi,1,INF),(pi,m+1,T,INF),(pi,j,pi,j+1,vi,j)。断掉 (pi,j,pi,j+1,vi,j)(p_{i,j},p_{i,j+1},v_{i,j})(pi,j,pi,j+1,vi,j) 的边即表示令 iii 选手做 jjj 号题。
这样建图之后限制就容易表示了,设 iii 选手比 jjj 选手题目编号至少大 kkk,就对于所有合法的下标 xxx,均连边 (pj,x,pi,x+k,INF)(p_{j,x},p_{i,x+k},INF)(pj,x,pi,x+k,INF) 即可。

但是这样是无法通过本题的。主流题解都说还需要加上 (pi,j+1,pi,j,INF)(p_{i,j+1},p_{i,j},INF)(pi,j+1,pi,j,INF) 的边,以实现限制的传递性(或者说保证每条链只断一条边),但我认为这并不是问题的本质,这个连边方法并不是必须的,实践结果也证明了这一点。
我们来看看原来的做法到底为什么出错。
原来的连边方式其实是很健全的,它本身已经有了传递性,如:
请添加图片描述红线和绿线限制已经可以进行传递,形成黄色的限制。
那么为什么我们难以通过一些测试点呢?
因为我们的限制没有加全!
题解和我一开始的写法对于限制的加边一般都是写成类似下边的形式:

for(int j=max(1,1-k);j<=m+1&&j+k<=m+1;j++) add(id[y][j],id[x][j+k],inf);

但正确的写法应该是:

for(int j=max(1,1-k);j<=m+1;j++) add(id[y][j],id[x][min(j+k,m+1)],inf);

这两个具象一下就是下面两张图的差别:
请添加图片描述请添加图片描述
注意到,在第一种加法中,后面的两个点没有被限制,成了“法外之地”,这也是为什么加完反向边可以把这个问题修正的原因。
即使可以这么修正,但从问题的本质来看还是不够优美的。因此,建议直接把加限制边的方式修改,而不是进行侧面的修正。

Code\text{Code}Code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define ok debug("OK\n")
inline ll read(){ll x(0),f(1);char c=getchar();while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}while(isdigit(c)){x=(x<<1)+(x<<3)+c-'0';c=getchar();}return x*f;
}
const int N=2e5+100;
const int M=2e6+100;
const double inf=1e11;
const double eps=1e-6;int n,m,num,y;
int flag;int tot,s,t;
struct node{int to,nxt;double cap;
}p[M<<1];
int fi[N],cur[N],cnt;
inline void addline(int x,int y,double c){p[++cnt]=(node){y,fi[x],c};fi[x]=cnt;
}
inline void add(int x,int y,double c){//printf("%d->%d cap=%.2lf\n\n",x,y,c);addline(x,y,c);addline(y,x,0);
}
int bel[N];
queue<int>que;
int bfs(){fill(bel,bel+1+tot,0);bel[s]=1;que.push(s);while(!que.empty()){int now=que.front();que.pop();for(int i=cur[now]=fi[now];~i;i=p[i].nxt){int to=p[i].to;if(p[i].cap<eps||bel[to]) continue;bel[to]=bel[now]+1;que.push(to);}}return bel[t];
}
double dfs(int x,double lim){if(x==t||lim<eps) return lim;double res(0);for(int &i=cur[x];~i;i=p[i].nxt){int to=p[i].to;if(bel[to]!=bel[x]+1) continue;double add=dfs(to,min(lim,p[i].cap));res+=add;lim-=add;p[i].cap-=add;p[i^1].cap+=add;//printf("x=%d to=%d add=%lf\n",x,to,add);if(lim<eps) break;}if(res<eps) bel[x]=-1;return res;
}
void dinic(){double flow(0),tmp(0);while(bfs()){while((tmp=dfs(s,inf))>eps) flow+=tmp;if(flow>=inf){puts("-1");return;}}printf("%lf\n",flow);return;
}
int id[90][90];
int c[90];
double ans;
double f[90][90][90],pp[90][90][90],val[90][90];
void init(){ans=0;tot=0;for(int i=1;i<=n;i++){for(int j=1;j<=m+1;j++) id[i][j]=++tot;}s=++tot;t=++tot;memset(fi,-1,sizeof(fi));cnt=-1;
}
void calc(){for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){pp[i][j][0]=1;val[i][j]=0;for(int k=1;k<=num;k++) pp[i][j][k]=pp[i][j][k-1]*f[i][j][k];for(int k=1;k<=num;k++){pp[i][j][k]=pp[i][j][k]*(1-f[i][j][k+1]);val[i][j]+=pp[i][j][k]*c[k];}//printf("i=%d j=%d val=%lf\n",i,j,val[i][j]);}}
}void work(){n=read();m=read();num=read();y=read();init();for(int i=1;i<=num;i++) c[i]=read()+c[i-1];for(int j=1;j<=m;j++){for(int i=1;i<=n;i++){for(int k=1;k<=num;k++) scanf("%lf",&f[i][j][k]);}}calc();for(int i=1;i<=n;i++){add(s,id[i][1],inf);add(id[i][m+1],t,inf);for(int j=1;j<=m;j++){add(id[i][j],id[i][j+1],val[i][j]);//add(id[i][j+1],id[i][j],inf);}}for(int i=1;i<=y;i++){int x=read(),y=read(),k=read();for(int j=max(1,1-k);j<=m+1;j++) add(id[y][j],id[x][min(j+k,m+1)],inf);}dinic();return;
}
signed main(){#ifndef ONLINE_JUDGE//freopen("a.in","r",stdin);//freopen("a.out","w",stdout);#endifint T=read();while(T--) work();return 0;
}
/*
1
3 2 4 1
10 10 10 20
0.3 0.5 0.7 0.4
0.2 0.6 0.2 0.2
0.7 0.1 0.8 0.2
0.5 0.5 0.5 0.5
0.2 0.5 0.3 0.6
0.3 0.5 0.4 0.1
2 3 1
*/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/317055.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Docker最全教程之Ubuntu下安装Docker(十五)

前言Ubuntu是一个以桌面应用为主的开源GNU/Linux操作系统&#xff0c;应用很广。本篇主要讲述Ubuntu下使用SSH远程登录并安装Docker&#xff0c;并且提供了Docker安装的两种方式&#xff0c;希望对大家有所帮助。拥抱Linux&#xff0c;大家可以从Ubuntu开始&#xff01;Ubuntu下…

群论学习笔记

文章目录前言群基本定义&#xff1a;子群陪集拉格朗日定理正规子群交换群商群阶置换定义置换的乘法循环置换群群作用等价类不动点Burnside引理内容证明法1 轨道-稳定子定理法2Polya 定理所谓群论&#xff0c;就是对群体行为问题的讨论。 &#xff08;逃&#xff09; 前言 个人…

Super Jumping! Jumping! Jumping! HDU - 1087

Super Jumping! Jumping! Jumping! HDU - 1087 题意&#xff1a; 给定一条长度为n的序列&#xff0c;其中一定存在一条元素和最大的严格上升子序列&#xff0c;求这条序列的元素和。 题解&#xff1a; 最长上升序列模板题 代码&#xff1a; #include<bits/stdc.h> …

Codeforces Round #727 (Div. 2) 题解

文章目录A. Contest StartB. Love SongC. Stable GroupsD. PriceFixedE. Game with CardsF. Strange Array#727-Div.2A. Contest Start 数学题&#xff0c;分类讨论 一般的&#xff0c;一段区间[l,r][l,r][l,r]会对后面固定人数造成影响&#xff0c;假设是kkk最后kkk个人&…

潘淳(寒树Office):不务正业的公众号满月了,都写了些啥?

三喜临门一 喜今天真是个好日子&#xff0c;“流浪太阳”又回来了&#xff0c;阴雨绵绵长恨无期&#xff0c;今天苏州终于天晴&#xff0c;于是心情大好&#xff01;都说好心情会带来好运气&#xff0c;冥冥感觉要写点啥了&#xff0c;果不其然今天还有另外两喜。大早起来得…

G List it all

传送 题意&#xff1a; 题解&#xff1a; 我们来考虑以下样例&#xff1a;1&#xff0c;1&#xff0c;2 我们先考虑1的贡献&#xff1a;如图(图中只花了) 2&#xff01;表示还剩两个空位&#xff0c;还有两个数未填入&#xff0c;所以是2&#xff01;个 对于n个数重复&#x…

洛谷P4727:图的同构计数(Polya引理)(dfs)

解析 《关于我想了半天 dp 结果看题解 dfs 就行这回事》 我就说 gcd⁡\gcdgcd 这玩意 dp 个锤子啊… 拆分数的增长速度远没有想像中那么大&#xff0c;事实上&#xff0c;n60n60n60 也就 1e6 左右。 据题解说&#xff0c;这玩意的增长速度仅有 O(enn)O(\frac{e^{\sqrt n}}{n})…

ASP.NET Core 实战:基于 Dapper 扩展你的数据访问方法

一、前言在非静态页面的项目开发中&#xff0c;必定会涉及到对于数据库的访问&#xff0c;最开始呢&#xff0c;我们使用 Ado.Net&#xff0c;通过编写 SQL 帮助类帮我们实现对于数据库的快速访问&#xff0c;后来&#xff0c;ORM&#xff08;Object Relational Mapping&#x…

Codeforces Round #699 (Div. 2) 题解

文章目录A. Space NavigationB. New ColonyC. Fence PaintingD. AB GraphE. Sorting BooksF. AB Tree#699-Div.2A. Space Navigation 对于最终位置(x,y)&#xff0c;我们只关心那两个方向的字符是否足够即可 #include <cstdio> #include <cstring> #define maxn …

洛谷P4271:New Barns P(倍增)(LCT)(直径)

解析 倍增真香 关键性质&#xff1a;树上距离一个点最远的点必定是直径两端点其一。 本题限制好&#xff0c;要求少动态维护倍增数组暴力维护直径即可。 如果每次合并的是两棵树&#xff0c;而不是一棵树加一个点&#xff0c;可以先离线下来&#xff0c;照样能做。 如果每次强…

HDU 5510 Bazinga

HDU 5510 Bazinga 题意&#xff1a; 依次给你n个字符串&#xff0c;让你找到编号最大的字符串&#xff0c;存在一个比他编号小的字符串且不是其子串 题解&#xff1a; string中有find查找功能&#xff0c; 思路是用一个vector来存之前所有字符串&#xff0c;数组book用来表…

微软发布 Visual Studio 2019年第二季度路线图

微软近日发布了 Visual Studio 2019 年第二季度的路线图&#xff0c;路线图介绍了目前 VS 致力于在 VS 2019 发布的一些重要功能。官方表示&#xff0c;Visual Studio 2019 将继续按照 Visual Studio 发行周期流程提供更新&#xff0c;也就是约每 6 周推出一次次要更新&#xf…

Educational Codeforces Round 107 (Rated for Div. 2) 题解

文章目录A. Review SiteB. GCD LengthC. Yet Another Card DeckD. Min Cost StringE. Colorings and DominoesF. ChainwordG. Chips on a BoardEducational-Round-107A. Review Site 都给了两台机子&#xff0c;直接把所有只会投②的扔到一台&#xff0c;其余的全是另一台 就…

李争——一个骨子里是极客的程序员

我的业余作品《IT 英雄传》&#xff0c;聚焦身边的英雄&#xff0c;以文字采访的形式记录奇人趣事&#xff0c;笑看风云变幻。所写的人都是我见过面且比较熟悉的&#xff0c;绝大部分都是交往很久的&#xff0c;其中为了避嫌&#xff0c;我很少写微软同事&#xff0c;但今天这一…

Pagodas HDU - 5512

Pagodas HDU - 5512 题意&#xff1a; 一开始给你两个数a和b&#xff0c;你可以得到c通过&#xff0c;cab&#xff0c;或者ca-b&#xff0c;你所能得到的数的范围是1~n&#xff0c;两个人轮流操作&#xff0c;当有一方无法操作时&#xff0c;另一方获胜 题解&#xff1a; c…

期望学习笔记

前言 突然发现自己没有系统学过期望。 做一本通的时候是从二分图开始听的课&#xff0c;dp这一章只是四处搜题解而已。 做期望题基本都是靠玄学和《感性理解》 都是很简单的东西&#xff0c;但系统很重要&#xff0c;该补的还是要补的。 期望的基本性质 E(c)cE(c)cE(c)cE(cx)…

仅此一文让你明白事务隔离级别、脏读、不可重复读、幻读

网络上关于这方面的博文有些偏理论&#xff0c;有些通篇代码&#xff0c;都不能深入浅出。本文用图文并茂的方式&#xff0c;配上行云流水般的代码&#xff0c;非要摆清楚这个问题。相关代码已提交至码云&#xff08;点击这里下载&#xff09;。事务是现代关系型数据库的核心之…

[AtCoder Regular Contest 123] 题解

文章目录A - Arithmetic SequenceB - Increasing TriplesC - 1, 2, 3 - DecompositionD - Inc, Dec - DecompositionE - TrainingF - Insert AdditionARC123A - Arithmetic Sequence 大讨论 只能111&#xff0c;先固定中间的数&#xff0c;看两边加谁&#xff0c;如果都是加负…

2.5:模拟总结

文章目录前言考场题目解析T1T2T3总结代码T1T2T3前言 50pts 30020 rnk19 … 把1000ms看成10s我也真是个人才。 T3自然溢出50带模数T成20有点离谱。 但倒没有因为WA失分。 就是菜罢了 考场 这次时间管理还是比较合理的。 乍看三题觉得T1似乎是个伞兵题 这离线下来可持久化数组…

Meeting HDU - 5521

Meeting HDU - 5521 题意&#xff1a; 一共有n个点&#xff0c;有m个块&#xff0c;每个块内有Si个点&#xff0c;块内点彼此到达费用为wi&#xff0c;两个人分别位于1和n号块&#xff0c;两者同时出发问最短时间遇到是多少&#xff1f;在哪些地方可以遇到&#xff1f; ΣSi&…