给出一种最小割的方法。

设\(num1[i]\),\(num2[i]\)为第i种形状的点心的两种口味的数量

设\(type[i]\),\(type[i]\)为第i种形状的点心的两种口味

假设\(num1[i]<num2[i]\)

考虑几种最优的购买方案:

1.买\(num1[i]+1\)个点心。这样一定可以得到\(type2[i]\)。

2.买\(num2[i]+1\)个点心。这样一定可以得到\(type1[i]\)和\(type2[i]\)。

3.不买。这样连毛都得不到。

然后根据这几个购买方案建图。

把每一个点i,拆成\(i\),\(i+n\)

设\(S\)为源点,\(T\)为汇点。

对于每一个\(i\),\(S\)到\(i\)连容量为\(num2[i]+1\)的边,\(i+n\)到\(T\)连容量为\(num2[i]+1\)的边。

对于每一个口味,设它对应的形状为\(i\)和\(j\),从\(i\)到\(j+n\)和\(j\)到\(i+n\)连容量为通过买形状为i或j的点心总而买到这个口味点心的最小代价的较小值。

然后跑最小割。

因为图是对称的,答案就是最小割/2

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<queue>
#include<vector>
using namespace std;
const int N=1010;
const int INF=1e9;
int n,cnt,head[N*2],ans,S,T,dis[N*2];
int read(){
int sum=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){sum=sum*10+ch-'0';ch=getchar();}
return sum*f;
}
struct node{
int type,num;
node(int t=0,int n=0){
type=t,num=n;
}
};
vector<node> vec[N];
struct edge{
int to,nxt,flow;
}e[N*8];
void add_edge(int u,int v,int flow){
cnt++;
e[cnt].nxt=head[u];
e[cnt].to=v;
e[cnt].flow=flow;
head[u]=cnt;
}
bool bfs(){
queue<int> q;
q.push(S);
memset(dis,-1,sizeof(dis));
dis[S]=0;
while(!q.empty()){
int u=q.front();
q.pop();
for(int i=head[u];i;i=e[i].nxt){
int v=e[i].to;
if(dis[v]==-1&&e[i].flow){
dis[v]=dis[u]+1;
q.push(v);
}
}
}
if(dis[T]==-1)return false;
return true;
}
int dfs(int u,int f){
if(u==T||f==0)return f;
int used=0;
for(int i=head[u];i;i=e[i].nxt){
int v=e[i].to;
if(dis[v]==dis[u]+1&&e[i].flow){
int w=dfs(v,min(f-used,e[i].flow));
if(w){
e[i].flow-=w;
e[i^1].flow+=w;
used+=w;
if(used==f)return f;
}
}
}
if(used==0)dis[u]=-1;
return used;
}
void Dinic(){
while(bfs())ans+=dfs(S,INF);
}
int main(){
n=read();
cnt=1;
S=0,T=n*2+1;
for(int i=1;i<=n;i++){
int type1=read()+1,num1=read(),type2=read()+1,num2=read();
if(num1>num2)swap(num1,num2),swap(type1,type2);
vec[type1].push_back(node(i,num2+1)) ;
vec[type2].push_back(node(i,num1+1));
add_edge(S,i,num2+1),add_edge(i,S,0);
add_edge(i+n,T,num2+1),add_edge(T,i+n,0);
}
for(int i=1;i<=n;i++){
int a=vec[i][0].type,b=vec[i][1].type,c=min(vec[i][0].num,vec[i][1].num);
add_edge(a,b+n,c);add_edge(b+n,a,0);
add_edge(b,a+n,c);add_edge(a+n,b,0);
}
Dinic();
printf("%d",ans/2);
return 0;
}

51nod 1321 收集点心(最小割)的更多相关文章

  1. 51nod 1499 (最小割)

    题意 分析 将一些点分成两个集合,很明显的最小割问题 设一个S.T,和S相连的点表示在B集合中,和T相连的点表示在A集合中 因为原题是完美值最大,我们转换一下,变成损失的价值最小,那么就是最小割问题了 ...

  2. BZOJ 1391: [Ceoi2008]order [最小割]

    1391: [Ceoi2008]order Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 1509  Solved: 460[Submit][Statu ...

  3. BZOJ-2127-happiness(最小割)

    2127: happiness(题解) Time Limit: 51 Sec  Memory Limit: 259 MBSubmit: 1806  Solved: 875 Description 高一 ...

  4. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  5. BZOJ3438 小M的作物(最小割)

    题目 Source http://www.lydsy.com/JudgeOnline/problem.php?id=3438 Description 小M在MC里开辟了两块巨大的耕地A和B(你可以认为 ...

  6. 最大流-最小割 MAXFLOW-MINCUT ISAP

    简单的叙述就不必了. 对于一个图,我们要找最大流,对于基于增广路径的算法,首先必须要建立反向边. 反向边的正确性: 我努力查找了许多资料,都没有找到理论上关于反向边正确性的证明. 但事实上,我们不难理 ...

  7. bzoj1412最小割

    太羞耻了,m n写反了(主要是样例n m相等) 建图方法比较高(ji)端(chu),对于可以加栅栏的地方连上1的边,然后求最小割即可 为了让代码优(suo)美(duan),我写了一个check,避免多 ...

  8. 【BZOJ1497】[NOI2006]最大获利 最小割

    裸的最小割,很经典的模型. 建图:要求总收益-总成本最大,那么将每条弧与源点相连,流量为成本,每个收益与汇点相连,流量为收益,然后每条弧与它所能到达的收益相连,流量为inf. 与源点相连的是未被选中的 ...

  9. 二分图&网络流&最小割等问题的总结

    二分图基础: 最大匹配:匈牙利算法 最小点覆盖=最大匹配 最小边覆盖=总节点数-最大匹配 最大独立集=点数-最大匹配 网络流: 技巧: 1.拆点为边,即一个点有限制,可将其转化为边 BZOJ1066, ...

随机推荐

  1. SQL中Group By的使用(转)

    1.概述 “Group By”从字面意义上理解就是根据“By”指定的规则对数据进行分组,所谓的分组就是将一个“数据集”划分成若干个“小区域”,然后针对若干个“小区域”进行数据处理. 2.原始表 3.简 ...

  2. org.xml.sax.SAXParseException: Failed to read schema document 的原因分析与解决方法

    现象: org.xml.sax.SAXParseException: schema_reference.4: Failed to read schema documen t 'http://www.s ...

  3. 安装idea

    1.下载idea https://www.jetbrains.com/idea/download/#section=linux 2.解压 sudo tar -zxvf  ideaIC-2018.3.2 ...

  4. qt quick中qml编程语言

    Qt QML 入门 — 使用C++定义QML类型 发表于 2013 年 3 月 11 日   注册C++类 注册可实例化的类型 注册不实例化的QML类型 附带属性 注册C++类 注册可实例化的类型 如 ...

  5. 使用展开操符作替代 .apply() (prefer-spread)

    在ES2015以前,你必须使用Function.prototype.apply()来调用可变函数. var args = [1, 2, 3, 4]; Math.max.apply(Math, args ...

  6. 【henuacm2016级暑期训练-动态规划专题 B】Coloring Trees

    [链接] 我是链接,点我呀:) [题意] 在这里输入题意 [题解] f[i][j][k]前i个位置,第i个位置放j这个颜色,然后形成了k个联通块的最小花费 分第i个位置有没有已经放颜色两种情况考虑. ...

  7. 【转】C# winform 加载网页 模拟键盘输入自动接入访问网络

    [转]C# winform 加载网页 模拟键盘输入自动接入访问网络 声明: 本文原创,首发于博客园 http://www.cnblogs.com/EasyInvoice/p/6070563.html  ...

  8. hbase的几种访问方式

    Hbase的访问方式 1.Native Java API:最常规和高效的访问方式: 2.HBase Shell:HBase的命令行工具,最简单的接口,适合HBase管理使用: 3.Thrift Gat ...

  9. ECMAScript 6新特性之Proxy

    ECMAScript 6中新增了一个全局构造函数:Proxy.该构造函数能够接收两个參数:一个目标对象.一个处理对象. 代码演示样例: var target = {}; var handler = { ...

  10. 对Java、C#转学swift的提醒:学习swift首先要突破心理障碍。

    网上非常多都说swift是一门新手友好的语言. 但以我当年从Java转学Ruby的经验,swift对于从Java.C#转来的程序猿实际并不友好.原因就在于原来总有一种错觉:一个语言最重要的就是严谨,而 ...