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

相关文章

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

poj 1287 Networking(prim or kruscal最小生成树)

题意给你点与点间距离,求最小生成树。 注意点是,两点之间可能有不同的路,输入的时候选择最小的,和之前有道最短路WA的题目类似。 prim代码: #include<stdio.h>const int MaxN = 51;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int P;int prim(){bool vis[MaxN];

poj 2349 Arctic Network uva 10369(prim or kruscal最小生成树)

题目很麻烦,因为不熟悉最小生成树的算法调试了好久。 感觉网上的题目解释都没说得很清楚,不适合新手。自己写一个。 题意:给你点的坐标,然后两点间可以有两种方式来通信:第一种是卫星通信,第二种是无线电通信。 卫星通信:任何两个有卫星频道的点间都可以直接建立连接,与点间的距离无关; 无线电通信:两个点之间的距离不能超过D,无线电收发器的功率越大,D越大,越昂贵。 计算无线电收发器D

poj 1502 MPI Maelstrom(单源最短路dijkstra)

题目真是长得头疼,好多生词,给跪。 没啥好说的,英语大水逼。 借助字典尝试翻译了一下,水逼直译求不喷 Description: BIT他们的超级计算机最近交货了。(定语秀了一堆词汇那就省略吧再见) Valentine McKee的研究顾问Jack Swigert,要她来测试一下这个系统。 Valentine告诉Swigert:“因为阿波罗是一个分布式共享内存的机器,所以它的内存访问