poj3308最小割

这题做的稀里糊涂 首先 建图不会,然后 dinic姿势不对 ,还不知道 为啥。。。。还有这尼玛怎么看出来是 最小割的

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <string>
#include <iostream>
#include <map>
#include <cstdlib>
#include <list>
#include <set>
#include <queue>
#include <stack>
#include<math.h>
using namespace std;
const double INF=1e10;
const double eps=1e-9;
double Map[222][222];
int level[222];
int s,e;
double Min(double a,double b)
{
    return a>b?b:a;
}
bool bfs()
{
    memset(level,0,sizeof(level));
    queue<int> q;q.push(s);
    level[s]=1;
    while(!q.empty()){
        int cur=q.front();q.pop();
        for(int i=0;i<=e;i++){
            if(Map[cur][i]>eps&&!level[i]){
                level[i]=level[cur]+1;
                q.push(i);
            }
        }
    }
    return level[e];
}

double dfs(int x,double val)
{
    double sum=0;
    if(x==e) return val;
    for(int i=0;i<=e;i++){
        if(val>=eps&&Map[x][i]>=eps&&level[i]==level[x]+1){
            double t=dfs(i,Min(val,Map[x][i]));
            sum+=t;  val-=t;
            Map[x][i]-=t;Map[i][x]+=t;
        }
    }
    if(sum<eps) level[x]=-1;//优化  此时的层次网络 进过某点的流量为0说明 不用再更新这点了即使还有边连
    return sum;
}

double dinic()
{
    double ans=0;
    while(bfs()) ans+=dfs(s,INF);
    return ans;
}
int main()
{
    int Icase;
    int n,m,l;
    scanf("%d",&Icase);
    while(Icase--){
        scanf("%d%d%d",&m,&n,&l);
        memset(Map,0,sizeof(Map));
        s=0;e=n+m+1;
        for(int i=1;i<=m;i++){
            double a;
            scanf("%lf",&a);
            a=log(a);
            Map[s][i]=a;
        }
        for(int i=m+1;i<e;i++){
           double a;
           scanf("%lf",&a);
           a=log(a);
           Map[i][e]=a;
        }
        for(int i=0;i<l;i++){
            int a;int b;
            scanf("%d%d",&a,&b);
            b+=m;
            Map[a][b]=INF;
        }
        printf("%.4f\n",exp(dinic()));
    }
    return 0;
}

poj3308最小割

时间: 2024-10-10 02:44:53

poj3308最小割的相关文章

poj3308 最小割

因为行可以了,那列就不行,所以根据行列建立最小割模型. 然后这题精妙之处在于把乘法取对数后转化为加法,瞬间就简单了. 保证精度,C++AC ,16MS G++WA. #include<stdio.h> #include<string.h> #include<queue> #include<math.h> #define maxn 120 #define INF 10000000 using namespace std; struct node { int t

poj3308 Paratroopers --- 最小点权覆盖-&gt;最小割

题目是一个很明显的二分图带权匹配模型, 添加源点到nx建边,ny到汇点建边,(nx,ny)=inf建边,求最小割既得最小点权覆盖. 在本题中由于求的是乘积,所以先全部取log转换为加法,最后再乘方回来. #include <iostream> #include <cstring> #include <string> #include <cstdio> #include <cmath> #include <algorithm> #inc

poj3308 Paratroopers --- 最小点权覆盖-&amp;gt;最小割

题目是一个非常明显的二分图带权匹配模型, 加入源点到nx建边,ny到汇点建边,(nx.ny)=inf建边.求最小割既得最小点权覆盖. 在本题中因为求的是乘积,所以先所有取log转换为加法,最后再乘方回来. #include <iostream> #include <cstring> #include <string> #include <cstdio> #include <cmath> #include <algorithm> #in

POJ3308 Paratroopers(最小割/最小点权覆盖)

把入侵者看作边,每一行每一列都是点,选取某一行某一列都有费用,这样问题就是选总权最小的点集覆盖所有边,就是最小点权覆盖. 此外,题目的总花费是所有费用的乘积,这时有个技巧,就是取对数,把乘法变为加法运算,最后再还原. 另外还可以从最小割的思路去这么理解: 每一行与源点相连,容量为该行的花费:每一列与汇点相连,容量为该列的花费:对于每个入侵者的坐标,该行该列连接一条容量INF的边. 要让源点汇点不连通,割边集必然与所有入侵者的行或列相关,而这样建模后的最小割就是最小的花费(容量INF的边必然不是最

poj3308--Paratroopers(最小割)

poj3308:题目链接 题目大意:给出一个n*m的矩阵,矩阵的有l个格子会出现外星人,每行的开头和每列的开头都可以装备武器,可以消灭该行或该列的所有外星人,但是每装备一种武器需要有花费,如果装备多种武器,需要的花费是各种花费的乘积.问消灭所有外星人的最小的花费. 输入:给出n m l 然后一行n个数,表示每行武器的花费,之后m个数,是每列的花费,最后是外星人的坐标. 思路: 首先从题意中就可以知道这是一个二分图的覆盖问题,求最小覆盖,但是让求的不是最少的武器数,而是花费数,所以求二分匹配的哪一

【BZOJ2039】【2009国家集训队】人员雇佣 [最小割]

人员雇佣 Time Limit: 20 Sec  Memory Limit: 259 MB[Submit][Status][Discuss] Description 作为一个富有经营头脑的富翁,小L决定从本国最优秀的经理中雇佣一些来经营自己的公司.这些经理相互之间合作有一个贡献指数,(我们用Ei,j表示i经理对j经理的了解程度),即当经理i和经理j同时被雇佣时,经理i会对经理j做出贡献,使得所赚得的利润增加Ei,j.当然,雇佣每一个经理都需要花费一定的金钱Ai,对于一些经理可能他做出的贡献不值得

hdoj 4289 Control 【拆点 求最小割】

Control Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2295    Accepted Submission(s): 961 Problem Description You, the head of Department of Security, recently received a top-secret informati

There is a war (hdu 2435 最小割+枚举)

There is a war Time Limit: 5000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 970    Accepted Submission(s): 277 Problem Description There is a sea. There are N islands in the sea. There are some directional

HDU 3657 Game(取数 最小割)经典

Game Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1065    Accepted Submission(s): 449 Problem Description onmylove has invented a game on n × m grids. There is one positive integer on each g