伞兵(Paratroopers)

时间限制: 1 Sec  内存限制: 128 MB

题目描述

公元 2500 年,地球和火星之间爆发了一场战争。最近,地球军队指挥官获悉火星入侵者将派一些伞兵来摧毁地球的兵工厂,兵工厂是一个 m×n 大小的网格。他还获悉每个伞兵将着陆的具体位置(行和列)。由于火星的伞兵个个都很强壮、而且组织性强,只要有一个伞兵存活了,就能摧毁整个兵工厂。因此,地球军队必须在伞兵着陆后瞬间全部杀死他们。

为了完成这个任务,地球军队需要利用高科技激光枪。他们能在某行(或某列)安装一架激光枪,一架激光枪能杀死该行(或该列)所有的伞兵。在第 i 行安装一架激光枪的费用是 Ri,在第 i 列安装的费用是 Ci。要安装整个激光枪系统,以便能同时开火,总的费用为这些激光枪费用的乘积。现在,你的任务是选择能杀死所有伞兵的激光枪,并使得整个系统的费用最小。

输入

输入文件的第 1 行为整数 T,表示测试数据的数目,接下来有 T 个测试数据。每个测试数据的第 1 行为 3 个整数 m、n 和 L,1≤m≤50,1≤n≤50,1≤L≤500,分别表示网格的行和列、以及伞兵的数目;接下来一行为 m 个大于或等于 1.0 的实数,第 i 个实数表示 Ri;再接下来一行为 n 个大于或等于 1.0 的实数,第 i 个实数表示 Ci;最后 L 行,每行为两个整数,描述了每个伞兵的着陆位置。

输出

对每个测试数据,输出搭建整个激光枪系统的最小费用,精确到小数点后面 4 位有效数字。

样例输入

1
4 4 5
2.0 7.0 5.0 2.0
1.5 2.0 2.0 8.0
1 1
2 2
3 3
4 4
1 4

样例输出

16.0000
题解:
首先可以看出这是一道网络流的题目,至于怎么输出所选边的乘积,稍后再说,这里先讲一讲建图。
1.虚构源点和汇点。
2.源点和行连边,汇点和列连边。
3.根据士兵的位置将相对应的行和列连边,权值为正无穷大。
然后讲一讲怎么输出乘积。
1.a*b*c=e^(log(a*b*c)) ,又log(a*b*c)=log(a)+log(b)+log(c)
2.a*b*c最小则log(a)+log(b)+log(c)最小,也就是以log(a),log(b),log(c)为边的图的网络流的最小割。
3.最大流=最小割,于是乎就成了求log(a),log(b),log(c)为边的图的最大流。
4.当然最后答案要用exp()还原一下。
5.完美!!!
代码仅供参考:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<queue>
#include<stack>
#include<ctime>
#include<vector>
using namespace std;
int n,m,l,t;
struct node
{
int next,to;
double cap;
}edge[];
int size=,head[];
void putin(int from,int to,double cap)
{
size++;
edge[size].to=to;
edge[size].cap=cap;
edge[size].next=head[from];
head[from]=size;
}
void in(int from,int to,double cap)
{
putin(from,to,cap);
putin(to,from,);
}
int dist[],numbs[];
void bfs(int src,int des)
{
int i;
for(i=;i<=n+m+;i++)dist[i]=n+m+;
memset(numbs,,sizeof(numbs));
queue<int>mem;
mem.push(des);
dist[des]=;numbs[]++;
while(!mem.empty())
{
int x=mem.front();mem.pop();
for(i=head[x];i!=-;i=edge[i].next)
{
int y=edge[i].to;
if(edge[i].cap==&&dist[y]==n+m+)
{
dist[y]=dist[x]+;
numbs[dist[y]]++;
mem.push(y);
}
}
}
return;
}
double dfs(int src,double flow,int des)
{
if(src==des)return flow;
int i,mindist=n+m+;
double low=;
for(i=head[src];i!=-;i=edge[i].next)
{
int y=edge[i].to;
if(edge[i].cap)
{
if(dist[y]==dist[src]-)
{
double t=dfs(y,min(flow-low,edge[i].cap),des);
edge[i].cap-=t;
edge[i^].cap+=t;
low+=t;
if(dist[]>=n+m+)return low;
if(flow==low)break;
}
mindist=min(mindist,dist[y]+);
}
}
if(!low)
{
if(!(--numbs[dist[src]]))dist[]=n+m+;
++numbs[dist[src]=mindist];
}
return low;
}
double ISAP(int src,int des)
{
double ans=;
bfs(src,des);
while(dist[]<n+m+)ans+=dfs(src,2e8,des);
return ans;
}
int main()
{
int i,j;
scanf("%d",&t);
while(t--)
{
size=;
memset(head,-,sizeof(head));
scanf("%d%d%d",&n,&m,&l);
for(i=;i<=n;i++)
{
double q;
scanf("%lf",&q);
in(,i,log(q));
}
for(i=;i<=m;i++)
{
double q;
scanf("%lf",&q);
in(n+i,n+m+,log(q));
}
for(i=;i<=l;i++)
{
int from,to;
scanf("%d%d",&from,&to);
in(from,n+to,);
}
double maxflow=ISAP(,n+m+);
printf("%.4lf\n",exp(maxflow));
}
return ;
}

伞兵(Paratroopers)的更多相关文章

  1. POJ 3308 Paratroopers(最小割EK(邻接表&矩阵))

    Description It is year 2500 A.D. and there is a terrible war between the forces of the Earth and the ...

  2. POJ - 3308 Paratroopers(最大流)

    1.这道题学了个单词,product 还有 乘积 的意思.. 题意就是在一个 m*n的矩阵中,放入L个敌军的伞兵,而我军要在伞兵落地的瞬间将其消灭.现在我军用一种激光枪组建一个防御系统,这种枪可以安装 ...

  3. ACM/ICPC 之 伞兵-最小割转最大流(POJ3308)

    //以行列建点,伞兵位置为单向边-利用对数将乘积转加法 //最小割转最大流 //Time:63Ms Memory:792K #include<iostream> #include<c ...

  4. Paratroopers

    Paratroopers Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 7881 Accepted: 2373 Descript ...

  5. POJ3308 Paratroopers(网络流)(最小割)

                                                     Paratroopers Time Limit: 1000MS   Memory Limit: 655 ...

  6. POJ 3308 Paratroopers(最小割EK)

    题目链接 题意 : 有一个n*m的矩阵,L个伞兵可能落在某些点上,这些点的坐标已知,需要在某些位置安上一些枪,然后每个枪可以将一行或者一列的伞兵击毙.把这种枪安装到不同行的行首.或者不同列的列首,费用 ...

  7. zoj 2874 &amp; poj 3308 Paratroopers (最小割)

    意甲冠军: 一m*n该网络的规模格.详细地点称为伞兵着陆(行和列). 现在,在一排(或列) 安装激光枪,激光枪可以杀死线(或塔)所有伞兵.在第一i安装一排 费用是Ri.在第i列安装的费用是Ci. 要安 ...

  8. POJ 3308 Paratroopers 最大流,乘积化和 难度:2

    Paratroopers Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 7267   Accepted: 2194 Desc ...

  9. POJ 3308 Paratroopers(最小点权覆盖)(对数乘转加)

    http://poj.org/problem?id=3308 r*c的地图 每一个大炮可以消灭一行一列的敌人 安装消灭第i行的大炮花费是ri 安装消灭第j行的大炮花费是ci 已知敌人坐标,同时消灭所有 ...

随机推荐

  1. iOS关于Cookie验证登录状态

    1.第一次进入应用,登录获取Cookie,此时如果用到的是AFN去获取接口数据,Cookie已经写入了,所以无需处理,每次请求的时候,会自动将该cookie传给后台去验证 2.将Cookie缓存到本地 ...

  2. Android VideoView使用小记

    在Android中播放视频一般采用VideoView,当然也可以自己使用MediaPlayer+SurfaceView,但是比较麻烦.这里记录一些我使用VideoView时的疑惑 1.如何监听播放完成 ...

  3. python3 time模块与datetime模块

    time模块 在Python中,通常有这几种方式来表示时间:1)时间戳 2)格式化的时间字符串 3)元组(struct_time)共九个元素.由于Python的time模块实现主要调用C库,所以各个平 ...

  4. 小谈-—ServletConfig对象和servletContext对象

    一.servletContext概述 servletContext对象是Servlet三大域对象之一,每个Web应用程序都拥有一个ServletContext对象,该对象是Web应用程序的全局对象或者 ...

  5. Java学习笔记——序列化和反序列化

    寒雨连江夜入吴,平明送客楚山孤. 洛阳亲友如相问,一片冰心在玉壶. --芙蓉楼送辛渐 持久化数据的第一种方式.在序列化之前也可以把数据打散逐行存储在文件中,然后在逐行读取. 比如定Student类 用 ...

  6. Maven学习-简介、安装

    Maven是一个项目管理工具,它包含了一个项目对象模型,一组标准集合,一个项目声明周期,一个依赖管理系统和用来运行定义在生命周期阶段中插件目标的逻辑.Maven采用了约定优于配置这一基本原则.在没有自 ...

  7. 为Play初学者准备的Scala基础知识

    1 前言 本文的主要目的是为了让Play Framework的初学者快速了解Scala语言,算是一篇Play Framework的入门前传吧.使用PlayFramework可以极大的提高开发效率,但是 ...

  8. 每天一道Java题[6]

    题目 String字符串怎么转换为Date,Date又怎么转换成String字符串 解答 String->Date 主要用到类SimpleDateFormat及其抽象父类DateFormat中的 ...

  9. OpenCV探索之路(六):边缘检测(canny、sobel、laplacian)

    边缘检测的一般步骤: 滤波--消除噪声 增强--使边界轮廓更加明显 检测--选出边缘点 Canny算法 Canny边缘检测算法被很多人推崇为当今最优秀的边缘检测算法,所以我们第一个就介绍他. open ...

  10. 软件工程期末考试 AHNU

    1.      数据流图:一种图形化技术,它描绘信息流和数据从输入移动到输出的过程中所经受的变换.在数据流图中没有任何具体的物理部件,它只是描绘数据在软件中流动和被处理的逻辑过程,是系统逻辑功能的图形 ...