POJ2195:Going Home(费用流入门)

http://poj.org/problem?id=2195

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <queue>
#include <math.h>
#define inf 0x3f3f3f3f
using namespace std;
struct node1//记录h或者m的坐标
{
    int x,y;
} a[402],b[402];
struct node
{
    int x,y,c,w;
    int next;
} eg[600001];
int n,m,tt,head[10001],v[10001],dis[10001],pre[10001],s,t;
char map[102][102];
void init()
{
    tt=0;
    memset(head,-1,sizeof(head));
    s=0;
    t=n*m+1;
}
void add(int xx,int yy,int cc,int ww)
{
    eg[tt].x=xx;
    eg[tt].y=yy;
    eg[tt].c=cc;
    eg[tt].w=ww;
    eg[tt].next=head[xx];
    head[xx]=tt++;
    eg[tt].x=yy;
    eg[tt].y=xx;
    eg[tt].c=0;
    eg[tt].w=-ww;
    eg[tt].next=head[yy];
    head[yy]=tt++;
}
int spfa(int s,int t)
{
    queue<int>q;
    for(int i=0; i<=t; i++)
    {
        v[i]=0;
        dis[i]=inf;
        pre[i]=-1;
    }
    v[s]=1;
    dis[s]=0;
    q.push(s);
    while(!q.empty())
    {
        int ff=q.front();
        q.pop();
        v[ff]=0;
        for(int i=head[ff]; i!=-1; i=eg[i].next)
        {
            int w1=eg[i].y;
            if(eg[i].c&&dis[w1]>dis[ff]+eg[i].w)
            {
                dis[w1]=dis[ff]+eg[i].w;
                pre[w1]=i;//用边表建图存储其信息
                if(!v[w1])
                {
                    v[w1]=1;
                    q.push(w1);
                }
            }
        }
    }
    if(dis[t]==inf)
    {
        return 0;
    }
    return 1;
}
void KM(int s,int t)
{
    int min1,ans=0;
    while(spfa(s,t)==1)
    {
        min1=inf;
        for(int i=pre[t]; i!=-1; i=pre[eg[i].x])
        {
            if(min1>=eg[i].c)
            {
                min1=eg[i].c;
            }
        }
        for(int i=pre[t]; i!=-1; i=pre[eg[i].x])
        {
            eg[i].c-=min1;
            eg[i+1].c+=min1;
        }
        ans+=min1*dis[t];
    }
    printf("%d\n",ans);
    return ;
}
int main()
{
    int t1,t2;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(n==0&&m==0) break;
        init();
        t1=0;
        t2=0;
        for(int i=1; i<=n; i++)
        {
            scanf("%s",map[i]+1);//因为要创建超级源点,所以从一比较好计算
            for(int j=1; j<=m; j++)
            {
                int l=(i-1)*m+j;//以点建图
                if(map[i][j]==‘m‘)
                {
                    add(s,l,1,0);
                    a[t1].x=i;
                    a[t1++].y=j;
                }
                else if(map[i][j]==‘H‘)
                {
                    add(l,t,1,0);
                    b[t2].x=i;
                    b[t2++].y=j;
                }
            }
        }
        for(int i=0; i<t1; i++)
        {
            int l1=(a[i].x-1)*m+a[i].y;//以点建图
            for(int j=0; j<t2; j++)
            {
                int l2=abs(a[i].x-b[j].x)+abs(a[i].y-b[j].y);
                int l3=(b[j].x-1)*m+b[j].y;//以点建图
                add(l1,l3,1,l2);
            }
        }
        KM(s,t);
    }
    return 0;
}
时间: 2024-11-05 04:53:17

POJ2195:Going Home(费用流入门)的相关文章

HDU3549(最大流入门模板题)

public static void main(String[] args) { String a=null; if("aa".equals(a))//这种情形,不出现空指针异常 //if(a.equals("aa"))//出现空指针异常 { System.out.println(true); } else { System.out.println(false); } } 上面的两句不同的比较语句测试,第一句不出现空指针异常,第二句出现. 所以在变量和常量比较的时候

POJ2195 Going Home[费用流|二分图最大权匹配]

Going Home Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 22088   Accepted: 11155 Description On a grid map there are n little men and n houses. In each unit time, every little man can move one unit step, either horizontally, or vertica

POJ2195(KB11-D 最小费用最大流)

Going Home Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 23515   Accepted: 11853 Description On a grid map there are n little men and n houses. In each unit time, every little man can move one unit step, either horizontally, or vertica

IO流入门-第十一章-PrintStream_PrintWriter

DataInputStream和DataOutputStream基本用法和方法示例 /* java.io.PrintStream:标准的输出流,默认打印到控制台,以字节方式 java.io.PrintWriter:以字符方式 */ import java.io.*; import java.util.Date; import java.text.*; public class PrintStreamTest01 { public static void main(String[] args) t

io流入门级别demo

import java.io.*; //import java.io.BufferedReader; //import java.io.BufferedWriter; //import java.io.File; //import java.io.FileInputStream; //import java.io.FileNotFoundException; //import java.io.FileOutputStream; //import java.io.FileReader; //imp

IO流入门-第十二章-ObjectInputStream_ObjectOutputStream

DataInputStream和DataOutputStream基本用法和方法示例,序列化和反序列化 import java.io.Serializable; //该接口是一个"可序列化"的 ,没有任何方法,是一个标识接口,还有Cloneable /* 标识接口的作用:标识作用,JVM如果看到对象实现了某个标识接口,会对它特殊待遇,会给该类添加一个属性,static final long serialVersionUID=xxx, 最好是自己指定一个唯一的,这样不会产生类的兼容问题.

IO流入门-第八章-BufferedWriter

BufferedWriter基本用法和方法示例 import java.io.*; public class BufferedWriterTest01 { public static void main(String[] args) throws Exception { //创建带有缓冲区的字符输出流 //BufferedWriter bw = new BufferedWriter(new FileWriter("temp04")); BufferedWriter bw = new B

最大流入门题目

HDU 3549 Flow Problem Time Limit: 5000/5000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)Total Submission(s): 8728    Accepted Submission(s): 4083 Problem Description Network flow is a well-known difficult problem for ACMers. Given a

IO流入门-第四部分-FileReader

FileReader基本用法和方法示例 /* java.io.Reader java.io.InputStreamReader 转换流(字节输入流---->字符输入流) java.io.FileReader 文件字符输入流 */ import java.io.*; public class FileReaderTest01 { public static void main(String[] args) throws Exception { //创建文件字符输入流 FileReader fr =

Drainage Ditches 最大流入门练习题,各种算法

Drainage Ditches 题目抽象:n个顶点,m条容量为ci的边组成的图,求源点为1,汇点为n的最大流. 分析:各种最大流算法. 1.ford() 1 #include <cstdio> 2 #include <algorithm> 3 #include <cmath> 4 #include <cstring> 5 using namespace std; 6 const int INF = 0x5fffffff; 7 const int MS =