费用流

这种棋盘模型大概都是网络流吧

首先我们知道棋子之间不会影响到达目标的步数,那么就好做了,枚举终点,然后就是最小权匹配了,因为就是寻找总和最小,然后费用流就行了。

#include<bits/stdc++.h>
using namespace std;
const int N = , inf = 0x3f3f3f3f;
struct data {
int x, y;
} a[N];
struct edge {
int nxt, to, f, c;
} e[N * N];
int n, S, T, tot, cnt = , ans, kase;
int head[N], dis[N], q[N], Map[N][N], pree[N], prevv[N], inq[N];
void link(int u, int v, int f, int c)
{
e[++cnt].nxt = head[u];
head[u] = cnt;
e[cnt].to = v;
e[cnt].f = f;
e[cnt].c = c;
}
void insert(int u, int v, int f, int c)
{
link(u, v, f, c);
link(v, u, , -c);
}
bool spfa()
{
int l = , r = ;
memset(dis, 0x3f3f, sizeof(dis));
dis[] = ;
q[++r] = ;
while(l <= r)
{
int u = q[l++];
inq[u] = ;
for(int i = head[u]; i; i = e[i].nxt) if(e[i].f && dis[e[i].to] > dis[u] + e[i].c)
{
pree[e[i].to] = i;
prevv[e[i].to] = u;
dis[e[i].to] = dis[u] + e[i].c;
if(!inq[e[i].to])
{
inq[e[i].to] = ;
q[++r] = e[i].to;
}
}
}
return dis[T] != 0x3f3f3f3f;
}
int getflow()
{
int now = T, delta = inf;
while(now)
{
delta = min(delta, e[pree[now]].f);
now = prevv[now];
}
now = T;
while(now)
{
e[pree[now]].f -= delta;
e[pree[now] ^ ].f += delta;
now = prevv[now];
}
return delta * dis[T];
}
int maxcostflow()
{
int ret = ;
while(spfa()) ret += getflow();
return ret;
}
void build()
{
memset(head, , sizeof(head));
cnt = ;
for(int i = ; i <= n; ++i) insert(S, i, , );
for(int i = ; i <= n; ++i)
for(int x = ; x <= n; ++x)
for(int y = ; y <= n; ++y) if(Map[x][y])
insert(i, Map[x][y], , abs(a[i].x - x) + abs(a[i].y - y));
for(int i = ; i <= n; ++i)
for(int j = ; j <= n; ++j) if(Map[i][j])
insert(Map[i][j], T, , );
}
int main()
{
while(scanf("%d", &n))
{
if(!n) break;
T = * n + ;
ans = inf;
for(int i = ; i <= n; ++i) scanf("%d%d", &a[i].x, &a[i].y);
for(int i = ; i <= n; ++i)
{
memset(Map, , sizeof(Map));
tot = n;
for(int j = ; j <= n; ++j) Map[i][j] = ++tot;
build();
ans = min(ans, maxcostflow());
memset(Map, , sizeof(Map));
tot = n;
for(int j = ; j <= n; ++j) Map[j][i] = ++tot;
build();
ans = min(ans, maxcostflow());
}
memset(Map, , sizeof(Map));
tot = n;
for(int i = ; i <= n; ++i) Map[i][i] = ++tot;
build();
ans = min(ans, maxcostflow());
memset(Map, , sizeof(Map));
tot = n;
for(int i = ; i <= n; ++i) Map[i][n - i + ] = ++tot;
build();
ans = min(ans, maxcostflow());
printf("Board %d: %d moves required.\n\n", ++kase, ans);
}
return ;
}

最新文章

  1. discuz ucenter无法连接数据库
  2. 第14章 使用DHCP动态管理主机地址
  3. hdu 1520
  4. iOS之走进精益编程01
  5. Heritrix源码分析(二) 配置文件order.xml介绍(转)
  6. SQL自定义函数split分隔字符串
  7. Hibernate学习笔记--------2.一多|多多的CRUD
  8. 怎样向IT行业的朋友说明《圣经》的重要性
  9. js+jquery+html实现在三种不通的情况下,点击图片放大的效果
  10. 添加EF上下文对象,添加接口、实现类以及无处不在的依赖注入(DI)
  11. 数据库 SQL Server2012安装步骤详解
  12. git应用套路
  13. BZOJ 4129: Haruna’s Breakfast [树上莫队 分块]
  14. 搭建 vue2 单元测试环境(karma+mocha+webpack3)
  15. PHP PSR代码规范
  16. LeetCode 389 Find the Difference 解题报告
  17. eclispe设置workspace text file encoding
  18. python学习 day14 (3月19日)----
  19. poj-2421-最小生成树刷题
  20. Mongodb 安装(Windows)

热门文章

  1. Hibernate自动事务揪出的编码不规范
  2. django 加日志
  3. linux安装mysql可视化工具MySQL-workbench 连接数据库 执行sql
  4. 发布自己的nuget包
  5. Python 不定长参数、全局变量、局部变量 day4
  6. c++/c DEBUG宏
  7. oralce 创建表空间 和 查询进程
  8. php第十六节课
  9. [如何在Mac下使用gulp] 1.创建项目及安装gulp
  10. poj3176-Cow Bowling【dp】