这个题太恶心了。。。并不想继续做了。。。

本代码在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——最小割的更多相关文章

  1. 【BZOJ2127】happiness 最小割

    题目大意:有一个$n\times m$的矩阵,矩阵的每个位置上有一个同学,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了,每个同学对于选择文科与理科有着自己的喜悦 ...

  2. [置顶] [BZOJ]2127: happiness 最小割

    happiness: Description 高一一班的座位表是个n*m的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了,每个同学对于选择文科与理科有着自己 ...

  3. [BZOJ2127]happiness-[网络流-最小割]

    Description 传送门 Solution 按照最小割的思路考虑. 根据题意,当两个人都选文(理),需要砍掉两个人都选理(文)的加成:如果两个人选的不一样,就要都砍掉. 这是一个网络流建模的套路 ...

  4. [国家集训队]happiness 最小割 BZOJ 2127

    题目描述 高一一班的座位表是个n*m的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了,每个同学对于选择文科与理科有着自己的喜悦值,而一对好朋友如果能同时选文 ...

  5. BZOJ 2127 / Luogu P1646 [国家集训队]happiness (最小割)

    题面 BZOJ传送门 Luogu传送门 分析 这道题又出现了二元关系,于是我们只需要解方程确定怎么连边就行了 假设跟SSS分在一块是选文科,跟TTT分在一块是选理科,先加上所有的收益,再来考虑如何让需 ...

  6. luogu P1646 [国家集训队]happiness (最小割)

    高一一班的座位表是个n*m的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了,每个同学对于选择文科与理科有着自己的喜悦值,而一对好朋友如果能同时选文科或者理科 ...

  7. BZOJ 2127: happiness [最小割]

    2127: happiness Time Limit: 51 Sec  Memory Limit: 259 MBSubmit: 1815  Solved: 878[Submit][Status][Di ...

  8. BZOJ 2127: happiness(最小割解决集合划分)

    Time Limit: 51 Sec  Memory Limit: 259 MBSubmit: 2350  Solved: 1138[Submit][Status][Discuss] Descript ...

  9. 【BZOJ2127】happiness(最小割)

    [BZOJ2127]happiness(最小割) 题面 Description 高一一班的座位表是个n*m的矩阵,经过一个学期的相处,每个同学和前后左右相邻的同学互相成为了好朋友.这学期要分文理科了, ...

随机推荐

  1. 如何在创建hive表格的python代码中导入外部文件

    业务场景大概是这样的,我要对用户博文进行分词(这个步骤可以看这篇文章如何在hive调用python的时候使用第三方不存在的库-how to use external python library in ...

  2. windows系统下用VScode配置远程编辑服务器文件的环境!通过Rmate方法

    虽然公司电脑win可以通过Xshell通过SSH远程连接家中内网linux服务器了,但是只能用vim编辑文件有点不爽. 于是上网查询,windows下使用vscode远程编辑服务器文件的办法.参照博文 ...

  3. 顺序链表的C风格实现

    //头文件 #ifndef _SEQLIST_H_ #define _SEQLIST_H_ //定义数据类型 typedef void SeqList; typedef void SeqListNod ...

  4. react+redux状态管理实现排序 合并多个reducer文件

    这个demo只有一个reducer 所以合并reducer这个demo用不到 ,但是我写出来这样大家以后可以用到,很好用,管理多个reducer,因为只要用到redux就不会只有一个reducer所以 ...

  5. python基础之装饰器扩展和迭代器

    wraps模块 让原函数保留原来的说明信息 import time import random from functools import wraps def auth(func): '''auth ...

  6. 一步一步学Linq to sql(一):预备知识

    什么是Linq to sql Linq to sql(或者叫DLINQ)是LINQ(.NET语言集成查询)的一部分,全称基于关系数据的 .NET 语言集成查询,用于以对象形式管理关系数据,并提供了丰富 ...

  7. 第四模块:网络编程进阶&数据库开发 练习

    练习题 基于queue模块实现线程池 import threading from multiprocessing import Queue class A(threading.Thread): def ...

  8. Android 快捷方式的创建与查询 快捷方式问题大全 获取快捷方式在Launcher数据库中的信息 Failed to find provider info for com.android.la

    /** * 创建添加快捷方式 * 其中需要设置的有: * 1. 快捷方式的标题 * 2. 快捷方式的图标 * 3. 点击快捷方式后的跳转 */ public static void createSho ...

  9. CC3200在AP模式的TCP sock作为客户端连接时返回SL_ECONNREFUSED(-111) Connection refused

    1. CC3200处于AP模式(电脑无线连接CC3200的WIFI信号),开启一个TCP socket,这个socket作为TCP客户端去连接TCP服务器端 struct sockaddr_in ad ...

  10. PowerShell脚本授权最佳实践

    [TechTarget中国原创] Windows PowerShell已成为微软在Windows Server上提供的首选管理界面.因为深度整合到Windows Server操作系统,PowerShe ...