HDU 4862 JUMP 最小费用最大流
2014 多校的B题,由于我不怎么搞图论,当时碰到这个题目,我怎么想都没往网络流方面弄,不过网络流真的是个好东西,对于状态多变,无法用动规或者数据结构来很好表示的时候,非常有用
这个题目要求每个点一定要访问到,并且每次访问的是没访问过的点,跳跃的方向为向右或者向下。
建图的时候,分成二分图,从一个超级源点向x部分建cap为1 cost为0的点,对所以可到达的点从x到y建cap为1,cost根据题目算出来,不过要算负值,因为我们要求得实际是最大费用,最后对结果求相反数即可。所有y部分的点对超级汇点T建cap为1,cost为0的点。
但是这样还是没满足只能走k次的条件,所以得额外建一个点,从S到该点建cap为k,cost为0的点,再用这个点对y部分所有点建cap为1,cost为0的点,因为每次从x部流到y部最终流到T的点,都是不包括起始点的,用这个额外点向y部分供应最多k次(不一定要k个)的起始点量即可
最后求得一定要流量为 N*M才算走通,此外,cost的相反数即为所求最大费用
用的是大白书上的最小费用最大模板,复杂度应该是n*m*某个系数
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
const int maxn = 310;
struct Edge
{
int from,to,cap,flow,cost;
};
struct MCMF
{
int n,m,s,t;
vector<Edge> edges;
vector<int> G[maxn];
int inq[maxn];
int d[maxn];
int p[maxn];
int a[maxn];
void init(int n){
this->n=n;
for (int i=0;i<=n;i++) G[i].clear();
edges.clear();
}
void add(int from,int to,int cap,int cost)
{
edges.push_back((Edge){from,to,cap,0,cost});
edges.push_back((Edge){to,from,0,0,-cost});
m=edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}
} TF;
char mat[15][15];
int N,M,K;
int id[15][15];
int S,T;
bool Bellman(int s,int t,int &flow,int &cost)
{
for (int i=0;i<=2*N*M+2;i++) TF.d[i]=(1<<30);
memset(TF.inq,0,sizeof TF.inq);
TF.d[s]=0;TF.inq[s]=1;TF.a[s]=1<<30;
queue<int> Q;
Q.push(s);
while (!Q.empty())
{
int u=Q.front();Q.pop();
TF.inq[u]=0;
for (int i=0;i<TF.G[u].size();i++){
Edge & e=TF.edges[TF.G[u][i]];
if (e.cap>e.flow && TF.d[e.to]>TF.d[u]+e.cost){
TF.d[e.to]=TF.d[u]+e.cost;
TF.p[e.to]=TF.G[u][i];
TF.a[e.to]=min(TF.a[u],e.cap-e.flow);
if (!TF.inq[e.to]) {Q.push(e.to);TF.inq[e.to]=1;}
}
}
}
if (TF.d[t]>=(1<<30)) return false;
flow+=TF.a[t];
cost+=TF.d[t];
int u=t;
while (u!=s){
TF.edges[TF.p[u]].flow+=TF.a[t];
TF.edges[TF.p[u]^1].flow-=TF.a[t];
u=TF.edges[TF.p[u]].from;
}
return true;
}
int mincost()
{
int flow=0,cost=0;
while (Bellman(S,T,flow,cost));
if (flow==N*M)
return -cost;
else return -1;
}
int main()
{
int t;
scanf("%d",&t);
int kase=0;
while (t--)
{
scanf("%d%d%d",&N,&M,&K);
S=0,T=2*N*M+2;
TF.init(T);
for (int i=1;i<=N;i++) scanf("%s",mat[i]+1);
for (int i=1;i<=N;i++){
for (int c=1;c<=M;c++){
for (int j=i+1;j<=N;j++){
int cost=0;
cost+=j-i-1;
if (mat[i][c]==mat[j][c]) cost-=mat[i][c]-'0';
TF.add((i-1)*M+c,N*M+(j-1)*M+c,1,cost);
}
for (int j=c+1;j<=M;j++){
int cost=0;
cost+=j-c-1;
if (mat[i][c]==mat[i][j]) cost-=mat[i][c]-'0';
TF.add((i-1)*M+c,N*M+(i-1)*M+j,1,cost);
}
}
}
for (int i=1;i<=N;i++){
for (int j=1;j<=M;j++){
TF.add(S,(i-1)*M+j,1,0);
TF.add(N*M+(i-1)*M+j,T,1,0);
}
}
TF.add(S,2*N*M+1,K,0);
for (int i=N*M+1;i<=2*N*M;i++){
TF.add(2*N*M+1,i,1,0);
}
printf("Case %d : %d\n",++kase,mincost());
}
return 0;
}
HDU 4862 JUMP 最小费用最大流的更多相关文章
- hdu 2686 Matrix 最小费用最大流
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2686 Yifenfei very like play a number game in the n*n ...
- hdu 4494 Teamwork 最小费用最大流
Teamwork Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=4494 ...
- POJ 2195 Going Home / HDU 1533(最小费用最大流模板)
题目大意: 有一个最大是100 * 100 的网格图,上面有 s 个 房子和人,人每移动一个格子花费1的代价,求最小代价让所有的人都进入一个房子.每个房子只能进入一个人. 算法讨论: 注意是KM 和 ...
- HDU 4862 Jump(更多的联合培训学校1)(最小费用最大流)
职务地址:pid=4862">HDU4862 最小费用流做的还是太少. 建图想不出来. . . 直接引用官方题解的话吧... 最小K路径覆盖的模型.用费用流或者KM算法解决,构造二部图 ...
- HDU 4862 Jump 费用流
又是一个看了题解以后还坑了一天的题…… 结果最后发现是抄代码的时候少写了一个负号. 题意: 有一个n*m的网格,其中每个格子上都有0~9的数字.现在你可以玩K次游戏. 一次游戏是这样定义的: 你可以选 ...
- hdu 1533 Going Home 最小费用最大流
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1533 On a grid map there are n little men and n house ...
- HDU 5988.Coding Contest 最小费用最大流
Coding Contest Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)To ...
- hdu 3667(拆边+最小费用最大流)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3667 思路:由于花费的计算方法是a*x*x,因此必须拆边,使得最小费用流模板可用,即变成a*x的形式. ...
- hdu 2686&&hdu 3376(拆点+构图+最小费用最大流)
Matrix Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Subm ...
随机推荐
- 二 基于java动态数组手写栈
package dataStucture2.stack; import dataStucture2.array.MyDynamicArray; /** * 基于动态数组手写栈 * 设计时,栈中仅栈顶对 ...
- C#学习之time控件和timer_tick事件 -----转载
Timer控件:Timer控件只有绑定了Tick事件,和设置Enabled=True后才会自动计时,停止计时可以用Stop()控制,通过Stop()停止之后,如果想重新计时,可以用Start()方法来 ...
- Spark教程——(5)PySpark入门
启动PySpark: [root@node1 ~]# pyspark Python 2.7.5 (default, Nov 6 2016, 00:28:07) [GCC 4.8.5 20150623 ...
- 十七、java内存模型_JVM_JDK_类加载
1.Java内存模型 共享内存模型指的就是Java内存模型(简称JMM),JMM决定一个线程对共享变量的写入时,能对另一个线程可见.从抽象的角度来看,JMM定义了线程和主内存之间的抽象关系:线程之间的 ...
- 「Luogu P2253 好一个一中腰鼓!」
就这道题的理论难度来说绿题是有点低了,但是这道题的实际难度来看,顶多黄题,所以建议加强数据或出数据升级版. 前置芝士 线段树:具体可以看我的另一篇文章. 具体做法 暴力的方法想必都会,所以来讲一下正解 ...
- English-Number
English-Number 1. 基数与序数 2. 数量级前缀 3. 与数量级前缀连用的常用单位 4. 数量的完整写法式例 5. 数学符号 6. 其他常用单位 7. 数字的用法 7.1 日期时间 7 ...
- 开源Web测试工具介绍
HtmlUnitHtmlUnit 是 JUnit 的扩展测试框架之一.HtmlUnit 将返回文档模拟成 HTML,这样您便可以直接处理这些文档了.HtmlUnit 使用例如 table.form 等 ...
- python2.7 操作ceph-cluster S3对象接口 实现: 上传 下载 查询 删除 顺便使用Docker装个owncloud 实现UI管理
python version: python2.7 需要安装得轮子: botofilechunkio command: yum install python-pip&& pip ...
- C++Review15_内存管理
一.野指针 定义指针变量时最好初始化为NULL: 内存回收后,指针也用完了,这时候也需要及时将指针置为NULL: 指针就像野狗一样,为了防止它乱指,除了在使用期间,别的时候都需要置为NULL.这样它就 ...
- 51nod 1380:夹克老爷的逢三抽一
1380 夹克老爷的逢三抽一 基准时间限制:1 秒 空间限制:131072 KB 分值: 320 难度:7级算法题 收藏 取消关注 又到了诺德县的百姓孝敬夹克大老爷的日子,带着数量不等的铜板的村民 ...