poj 3281 Dining 最大流dinic 模板题

2024-01-13 02:48
文章标签 模板 最大 poj 3281 dining dinic

本文主要是介绍poj 3281 Dining 最大流dinic 模板题,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

http://poj.org/problem?id=3281
题意:有n头牛,每头牛喜欢吃几种对应的饮料和食物(不同牛喜欢的种类可能不同),求怎样分配
使得能同时吃到喜欢的饮料和食物的牛的个数最多;

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <queue>
#include <map>
#include <algorithm>
#include <set>
using namespace std;
#define MM(a) memset(a,0,sizeof(a))
typedef long long ll;
typedef unsigned long long ULL;
const int mod = 1000000007;
const double eps = 1e-10;
const int inf = 0x3f3f3f3f;
const int big=50000;
int max(int a,int b) {return a>b?a:b;};
int min(int a,int b) {return a<b?a:b;};
struct Edge{
int to,cap,rev;
};
int level[500];
vector<Edge> G[500];
void add_edge(int f,int t)
{
G[f].push_back((Edge){t,1,G[t].size()});
G[t].push_back((Edge){f,0,G[f].size()-1});
}
void bfs(int s)
{
queue<int> q;
q.push(s);
level[s]=1;

//cout<<s<<" "<<level[s]<<endl;

while(q.size())
{
int now=q.front();q.pop();
for(int i=0;i<G[now].size();i++)
if(G[now][i].cap>0)
{
Edge e=G[now][i];
//printf("%d: %d\n",now,e.to);
if(level[e.to]<0)
{
level[e.to]=level[now]+1;
//printf("%d %d: %d %d\n",now,level[now],e.to,level[e.to]);
q.push(e.to);
}
}
}
}
int dfs(int s,int t)
{
if(s==t)
return 1;
for(int i=0;i<G[s].size();i++)
{
Edge &e=G[s][i];
if(level[e.to]>level[s]&&e.cap>0&&dfs(e.to,t))
{
e.cap-=1;
G[e.to][e.rev].cap+=1;
return 1;
}
}
return 0;
}
int max_flow(int s,int t)
{
int ans=0;
for(;;)
{
memset(level,-1,sizeof(level));
bfs(s);
//printf("%d\n",level[t]);
if(level[t]<0)
return ans;
while(dfs(s,t))
ans++;
}
}
void init()
{
for(int i=0;i<=500;i++)
G[i].clear();
}
int main()
{
int n,f,d,fx,dx,temp;
while(~scanf("%d %d %d",&n,&f,&d))
{
init();
for(int i=1;i<=n;i++)
{
scanf("%d %d",&fx,&dx);
for(int j=1;j<=fx;j++)
{
scanf("%d",&temp);
add_edge(temp,f+i);
}
for(int j=1;j<=dx;j++)
{
scanf("%d",&temp);
add_edge(i+n+f,2*n+f+temp);
}
}
for(int i=1;i<=f;i++)
add_edge(0,i);
for(int i=1;i<=d;i++)
add_edge(f+2*n+i,f+d+2*n+1);
for(int i=1;i<=n;i++)
add_edge(i+f,i+n+f);

/*
for(int i=0;i<=f+d+2*n+1;i++)
for(int j=0;j<G[i].size();j++)
if(G[i][j].cap>0)
printf("%d: %d %d\n",i,G[i][j].to,G[i][j].cap);
*/

printf("%d\n",max_flow(0,f+d+2*n+1));
}
return 0;
}

分析:因为同时有饮料喝食物所以不能二分匹配,只能跑最大流;
wa代码:
   
  • #include <iostream>
    #include <cstdio>
    #include <cstring>
    #include <cstdlib>
    #include <cmath>
    #include <vector>
    #include <queue>
    #include <map>
    #include <algorithm>
    #include <set>
    using namespace std;
    #define MM(a) memset(a,0,sizeof(a))
    typedef long long ll;
    typedef unsigned long long ULL;
    const int mod = 1000000007;
    const double eps = 1e-10;
    const int inf = 0x3f3f3f3f;
    const int big=50000;
    int max(int a,int b) {return a>b?a:b;};
    int min(int a,int b) {return a<b?a:b;};
    struct Edge{
    int to,cap,rev;
    };
    int level[500];
    vector<Edge> G[500];
    void add_edge(int f,int t)
    {
    G[f].push_back((Edge){t,1,G[t].size()});
    G[t].push_back((Edge){f,0,G[f].size()-1});
    }
    void bfs(int s)
    {
    queue<int> q;
    q.push(s);
    level[s]=1;
    while(q.size())
    {
    int now=q.front();q.pop();
    for(int i=0;i<G[now].size();i++)
    {
    Edge e=G[now][i];
    if(level[e.to]<0&&e.cap>0)
    {
    level[e.to]=level[s]+1;
    q.push(e.to);
    }
    }
    }
    }
    int dfs(int s,int t)
    {
    if(s==t)
    return 1;
    for(int i=0;i<G[s].size();i++)
    {
    Edge e=G[s][i];
    if(level[e.to]>level[s]&&dfs(e.to,t))
    {
    G[s][i].cap-=1;
    G[e.to][G[s][i].rev].cap+=1;
    return 1;
    }
    }
    return 0;
    }
    void max_flow(int s,int t)
    {
    memset(level,-1,sizeof(level));
    int ans=0;
    for(;;)
    {
    bfs(s);
    if(level[t]<0)
    break;
    while(dfs(s,t))
    ans++;
    }
    printf("%d\n",ans);
    }
    int main()
    {
    int n,f,d,fx,dx,temp;
    while(~scanf("%d %d %d",&n,&f,&d))
    {
    for(int i=1;i<=n;i++)
    {
    scanf("%d %d",&fx,&dx);
    for(int j=1;j<=fx;j++)
    {
    scanf("%d",temp);
    add_edge(temp,f+i);
    }
    for(int j=1;j<=dx;j++)
    {
    scanf("%d",&temp);
    add_edge(i+n+f,2*n+f+temp);
    }
    }
    for(int i=1;i<=f;i++)
    add_edge(0,i);
    for(int i=1;i<=d;i++)
    add_edge(f+2*n+i,f+d+2*n+1);
    max_flow(0,f+d+2*n+1);
    }
    return 0;
    }


这篇关于poj 3281 Dining 最大流dinic 模板题的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/599995

相关文章

C++中函数模板与类模板的简单使用及区别介绍

《C++中函数模板与类模板的简单使用及区别介绍》这篇文章介绍了C++中的模板机制,包括函数模板和类模板的概念、语法和实际应用,函数模板通过类型参数实现泛型操作,而类模板允许创建可处理多种数据类型的类,... 目录一、函数模板定义语法真实示例二、类模板三、关键区别四、注意事项 ‌在C++中,模板是实现泛型编程

基于Java实现模板填充Word

《基于Java实现模板填充Word》这篇文章主要为大家详细介绍了如何用Java实现按产品经理提供的Word模板填充数据,并以word或pdf形式导出,有需要的小伙伴可以参考一下... Java实现按模板填充wor编程d本文讲解的需求是:我们需要把数据库中的某些数据按照 产品经理提供的 word模板,把数据

如何提高Redis服务器的最大打开文件数限制

《如何提高Redis服务器的最大打开文件数限制》文章讨论了如何提高Redis服务器的最大打开文件数限制,以支持高并发服务,本文给大家介绍的非常详细,感兴趣的朋友跟随小编一起看看吧... 目录如何提高Redis服务器的最大打开文件数限制问题诊断解决步骤1. 修改系统级别的限制2. 为Redis进程特别设置限制

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

C++11第三弹:lambda表达式 | 新的类功能 | 模板的可变参数

🌈个人主页: 南桥几晴秋 🌈C++专栏: 南桥谈C++ 🌈C语言专栏: C语言学习系列 🌈Linux学习专栏: 南桥谈Linux 🌈数据结构学习专栏: 数据结构杂谈 🌈数据库学习专栏: 南桥谈MySQL 🌈Qt学习专栏: 南桥谈Qt 🌈菜鸡代码练习: 练习随想记录 🌈git学习: 南桥谈Git 🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈�

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

poj 1511 Invitation Cards(spfa最短路)

题意是给你点与点之间的距离,求来回到点1的最短路中的边权和。 因为边很大,不能用原来的dijkstra什么的,所以用spfa来做。并且注意要用long long int 来存储。 稍微改了一下学长的模板。 stack stl 实现代码: #include<stdio.h>#include<stack>using namespace std;const int M

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n