[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的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了, ...
随机推荐
- Hive的数据库和表
本文介绍一下Hive中的数据库(Database/Schema)和表(Table)的基础知识,由于篇幅原因,这里只是一些常用的.基础的. Hive的数据库和表 先看一张草图: Hive结构 从图上可以 ...
- ArrayList & Vector的源码实现
#ArrayList & Vector #####前言: 本来按照计划,ArrayList和Vector是分开讲的,但是当我阅读了ArrayList和Vector的源码以后,我就改变了注意,把 ...
- windows禁用/启用hyper-V,解决hyper-V与模拟器同时启用时造成冲突
- Android 自定义光标样式
今天自定义光标,自己切图,不过怎么切都是很宽.不是一个很细的条.我用ps花了一个像素的直线,放上去还是不行.后来在网上找到方法,那就是用shape.不得不说,shape真的是太吊了. 给EditTex ...
- 把实体bean对象转换成DBObject工具类
import java.lang.reflect.Field; import java.lang.reflect.InvocationTargetException; import java.util ...
- TerminateProcess
Remarks The TerminateProcess function is used to unconditionally cause a process to exit. The state ...
- centos使用--centos7.3配置LNMP
目录 1 源的配置 2 安装软件 2.1 安装php7 2.2 安装nginx 2.3 安装mysql 2.4 安装vsftp (ftp登录配置) 3 开机启动设置 4 其它一些配置 4.1 git的 ...
- DataGridView重查后,返回原来所在行
首先记录选中行 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 //查询前记录选中行 int _currentRow = 0; //int _cu ...
- Google浏览器历史版本下载地址和驱动器对应关系地址分享
Google浏览器历史版本下载地址https://www.slimjet.com/chrome/google-chrome-old-version.php google webdriver下载地址分享 ...
- sentry
https://docs.sentry.io/quickstart/?platform=javascript