[bzoj2127]happiness——最小割
这个题太恶心了。。。并不想继续做了。。。
本代码在bzoj上TLE!
大致说一下思路: 建立ST,首先由S连边(S,u,a)a代表学文的分数,连向T(u,T,b)b表示学理的分数,这样构造出了两个人独立的分数。
然后考虑联合分数,对于相邻的两个点xy,看下图(盗个图:
设xy都学文的分数为w1,都学理的分数为w2,则a=w1/2,b=w1/2,c=w2/2,d=w2/2,e=(w1+w2)/2,每一种割与其对应的亏损分数如下:
a+b -w1 都学理->w2
c+d -w2 都学文->w1
a+d+e -w1-w2 不同-> 0
c+d+e -w1-w2 ...
注意双向边e,我们是变成两条有向边加入网络,而又因为我们求最小割用的是最大流的算法,所以这条边可以看作是一条双向且权值为e的边。
然后把权值*2,解决精度问题。
#include <bits/stdc++.h>
using namespace std;
const int maxn = 105; /////////////////////////////////////////////
const int maxv = maxn * maxn * 2;
const int inf = INT_MAX;
int n, m, s, t, v;
struct edge {
int from;
int to;
int cap;
};
vector<edge> edges;
vector<int> G[maxv];
int dist[maxv], iter[maxv];
int z[maxv][maxv];
bool zz[maxv][maxv];
inline int read() {
char c = getchar();
int f = 1, x = 0;
while (!isdigit(c)) {
if (c == '-')
f = -1;
c = getchar();
}
while (isdigit(c))
x = x * 10 + c - '0', c = getchar();
return x * f;
}
inline void add_edge(int from, int to, int cap) {
edges.push_back((edge){from, to, cap});
edges.push_back((edge){to, from, 0});
int m = edges.size();
G[from].push_back(m - 2);
G[to].push_back(m - 1);
}
inline void bfs(int s) {
memset(dist, -1, sizeof(dist));
dist[s] = 0;
queue<int> q;
q.push(s);
while (!q.empty()) {
int u = q.front();
q.pop();
for (int i = 0; i < G[u].size(); i++) {
edge &e = edges[G[u][i]];
if (e.cap > 0 && dist[e.to] == -1) {
dist[e.to] = dist[u] + 1;
q.push(e.to);
}
}
}
}
inline int dfs(int s, int t, int flow) {
if (s == t)
return flow;
for (int &i = iter[s]; i < G[s].size(); i++) {
edge &e = edges[G[s][i]];
if (e.cap > 0 && dist[e.to] > dist[s]) {
int d = dfs(e.to, t, min(flow, e.cap));
if (d > 0) {
e.cap -= d;
edges[G[s][i] ^ 1].cap += d;
return d;
}
}
}
return 0;
}
inline int dinic(int s, int t) {
int flow = 0;
while (1) {
bfs(s);
if (dist[t] == -1)
return flow;
memset(iter, 0, sizeof(iter));
int F;
while (F = dfs(s, t, inf))
flow += F;
}
}
int main() {
// freopen("input", "r", stdin); //////////////////////////
memset(z, 0, sizeof(z));
memset(zz, 0, sizeof(zz));
scanf("%d %d", &n, &m);
s = 0; //文科
t = n * m + 1; //理科
v = t + 1;
int ans = 0;
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++) {
int x;
x = read();
x <<= 1;
z[s][(i - 1) * m + j] += x;
zz[s][(i - 1) * m + j] = 1;
}
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++) {
int x;
x = read();
x <<= 1;
z[(i - 1) * m + j][t] += x;
zz[(i - 1) * m + j][t] = 1;
}
for (int i = 1; i < n; i++) {
for (int j = 1; j <= m; j++) {
int x;
x = read();
z[(i - 1) * m + j][i * m + j] += x;
z[s][(i - 1) * m + j] += x;
z[s][i * m + j] += x;
zz[(i - 1) * m + j][i * m + j] = 1;
zz[s][(i - 1) * m + j] = 1;
zz[s][i * m + j] = 1;
add_edge((i - 1) * m + j, i * m + j, x);
add_edge(i * m + j, (i - 1) * m + j, x);
}
}
for (int i = 1; i < n; i++) {
for (int j = 1; j <= m; j++) {
int x;
x = read();
z[(i - 1) * m + j][i * m + j] += x;
z[(i - 1) * m + j][t] += x;
z[i * m + j][t] += x;
zz[(i - 1) * m + j][i * m + j] = 1;
zz[(i - 1) * m + j][t] = 1;
zz[i * m + j][t] = 1;
add_edge((i - 1) * m + j, i * m + j, x);
add_edge(i * m + j, (i - 1) * m + j, x);
}
}
for (int i = 1; i <= n; i++) {
for (int j = 1; j < m; j++) {
int x;
x = read();
z[(i - 1) * m + j][(i - 1) * m + j + 1] += x;
z[s][(i - 1) * m + j] += x;
z[s][(i - 1) * m + j + 1] += x;
zz[(i - 1) * m + j][(i - 1) * m + j + 1] = 1;
zz[s][(i - 1) * m + j] = 1;
zz[s][(i - 1) * m + j + 1] = 1;
add_edge((i - 1) * m + j, (i - 1) * m + j + 1, x);
add_edge((i - 1) * m + j + 1, (i - 1) * m + j, x);
}
}
for (int i = 1; i <= n; i++) {
for (int j = 1; j < m; j++) {
int x;
x = read();
z[(i - 1) * m + j][(i - 1) * m + j + 1] += x;
z[(i - 1) * m + j][t] += x;
z[(i - 1) * m + j + 1][t] += x;
zz[(i - 1) * m + j][(i - 1) * m + j + 1] = 1;
zz[(i - 1) * m + j][t] = 1;
zz[(i - 1) * m + j + 1][t] = 1;
add_edge((i - 1) * m + j, (i - 1) * m + j + 1, x);
add_edge((i - 1) * m + j + 1, (i - 1) * m + j, x);
}
}
for (int i = 1; i < t; i++) {
ans += z[s][i];
ans += z[i][t];
if (zz[s][i])
add_edge(s, i, z[s][i]);
if (zz[i][t])
add_edge(i, t, z[i][t]);
}
for (int i = 1; i < t; i++) {
for (int j = 1; j < t; j++) {
if (zz[i][j]) {
add_edge(i, j, z[i][j]);
add_edge(j, i, z[i][j]);
}
}
}
/* for (int i = 0; i < v; i++) {
cout << "For " << i << ':' << endl;
for (int j = 0; j < G[i].size(); j++) {
edge &e = edges[G[i][j]];
if (e.cap > 0)
cout << "to " << e.to << " cap " << e.cap << endl;
}
}*/
ans -= dinic(s, t);
printf("%d\n", ans >> 1);
}
[bzoj2127]happiness——最小割的更多相关文章
- 【BZOJ2127】happiness 最小割
题目大意:有一个$n\times m$的矩阵,矩阵的每个位置上有一个同学,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了,每个同学对于选择文科与理科有着自己的喜悦 ...
- [置顶] [BZOJ]2127: happiness 最小割
happiness: Description 高一一班的座位表是个n*m的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了,每个同学对于选择文科与理科有着自己 ...
- [BZOJ2127]happiness-[网络流-最小割]
Description 传送门 Solution 按照最小割的思路考虑. 根据题意,当两个人都选文(理),需要砍掉两个人都选理(文)的加成:如果两个人选的不一样,就要都砍掉. 这是一个网络流建模的套路 ...
- [国家集训队]happiness 最小割 BZOJ 2127
题目描述 高一一班的座位表是个n*m的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了,每个同学对于选择文科与理科有着自己的喜悦值,而一对好朋友如果能同时选文 ...
- BZOJ 2127 / Luogu P1646 [国家集训队]happiness (最小割)
题面 BZOJ传送门 Luogu传送门 分析 这道题又出现了二元关系,于是我们只需要解方程确定怎么连边就行了 假设跟SSS分在一块是选文科,跟TTT分在一块是选理科,先加上所有的收益,再来考虑如何让需 ...
- luogu P1646 [国家集训队]happiness (最小割)
高一一班的座位表是个n*m的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了,每个同学对于选择文科与理科有着自己的喜悦值,而一对好朋友如果能同时选文科或者理科 ...
- BZOJ 2127: happiness [最小割]
2127: happiness Time Limit: 51 Sec Memory Limit: 259 MBSubmit: 1815 Solved: 878[Submit][Status][Di ...
- BZOJ 2127: happiness(最小割解决集合划分)
Time Limit: 51 Sec Memory Limit: 259 MBSubmit: 2350 Solved: 1138[Submit][Status][Discuss] Descript ...
- 【BZOJ2127】happiness(最小割)
[BZOJ2127]happiness(最小割) 题面 Description 高一一班的座位表是个n*m的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了, ...
随机推荐
- PHP时间日期操作增减(date strtotime) 加一天 加一月
date_default_timezone_set('PRC'); //默认时区 //当前的时间增加5天 $date1 = "2014-11-11"; echo date('Y-m ...
- Android面试收集录16 Android动画总结
一.Android 动画分类 总的来说,Android动画可以分为两类,最初的传统动画和Android3.0 之后出现的属性动画: 传统动画又包括 帧动画(Frame Animation)和补间动画( ...
- Android面试收集录7 AsyncTask详解
1.Android中的线程 在操作系统中,线程是操作系统调度的最小单元,同时线程又是一种受限的系统资源,即线程不可能无限制地产生, 并且 **线程的创建和销毁都会有相应的开销.**当系统中存在大量的线 ...
- TouTiao开源项目 分析笔记16 新闻评论
1.要达到的效果 1.1.主要效果图 点击了标题栏的消息图标后,然后会跳转到评论详情的页面. 1.2.触发的点击事件 在新闻详情的片段中的菜单点击事件中 设置上方标题栏的消息标的监听事件 case R ...
- 3 破解密码,xshell连接
1.破解root密码 (1)启动电脑,按上下键进入启动菜单界面,选择第二个Red Hat Enterprise Linux Server, with Linux 0-rescue-* (2)按 'e' ...
- Android 布局错乱 Android花屏
最近做项目,妈的,有个一个很难受的bug. 这个bug ,自己这里没有手机,没有办法复现,找到了手机之后.解决了. 我先给大家看下什么叫布局错乱,花屏: 来张正常的图片: 正常情况下是这样的.然后, ...
- async/await 实现协程
2. 基本了解 在了解异步协程之前,我们首先得了解一些基础概念,如阻塞和非阻塞.同步和异步.多进程和协程. 2.1 阻塞 阻塞状态指程序未得到所需计算资源时被挂起的状态.程序在等待某个操作完成期间,自 ...
- Android学习记录(7)—Intent中显示意图和隐式意图的用法
Intent(意图)主要是解决Android应用的各项组件之间的通讯. Intent负责对应用中一次操作的动作.动作涉及数据.附加数据进行描述,Android则根据此Intent的描述,负责找到对应的 ...
- IIS 部署网站--浏览--“该页无法显示”
解决办法: 打开IIS管理器--web站点(网站)--右键点击对应的站点--属性--主目录--执行权限改为(脚本和和执行文件) 点击“应用”--确定. 重启一下站点,OK.
- 程序员必备PC维修法(软件篇)
学会使用专业软件检测与修复电脑硬件故障问题也是程序员的一种软技能. windows篇 情景:如何获取电脑硬件的真实信息.(如何检验选购回来的硬件是否正品) 自检:使用AIDA64软件检查电脑硬件,能详 ...