洛谷 [P3973] 线性代数
最大权闭合子图,神题
这不是线性代数,这是网络流。
我们看见这是一堆矩阵的运算,而且最后变成了一个数,那么我们就想到,把这个矩阵乘法的过程用具体的数字推出来
我们发现,a是一个01矩阵,然后其实就可以化成这么一个问题:
有n个东西,选了i,j两件东西能得到b[i,j]的价值,然而选i需要c[i]的花费,选j需要c[j]的花费……
这是一个经典的最小割模型,最大权闭合子图,详见胡伯涛论文。
建立S,T。
S连(i,j)边,边权为b[i,j],(i,j)连i、连j边,边权均为∞,i向T连边,边权为c[i]。
然后求最小割,最后答案就是
sum(b[i][j])-最小割答案 (i∈[1..n],j∈[1..n])
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <bitset>
#include <queue>
#define inf 0x3f3f3f3f
using namespace std;
const int MAXN = 300005;
int n, b[505][505], c[505], s, t, head[MAXN], nume, MaxFlow, ans, dep[MAXN], cur[MAXN];
int init() {
int rv = 0, fh = 1;
char c = getchar();
while(c < '0' || c > '9'){
if(c == '-') fh = -1;
c = getchar();
}
while(c >= '0' && c <= '9'){
rv = (rv<<1) + (rv<<3) + c -'0';
c = getchar();
}
return fh * rv;
}
struct edge{
int to, nxt, flow, cap;
}e[MAXN<<4];
void adde(int from, int to, int cap){
e[++nume].to = to;
e[nume].nxt = head[from];
e[nume].cap = cap;
head[from] = nume;
}
queue <int> q;
bool bfs(){
q.push(s);
memset(dep,0,sizeof(dep));
dep[s]=1;
while(!q.empty()){
int u = q.front();q.pop();
for(int i = head[u] ; i ; i = e[i].nxt){
int v = e[i].to;
if(!dep[v]&&e[i].flow < e[i].cap){
dep[v] = dep[u] + 1;
q.push(v);
}
}
}
return dep[t];
}
int dfs(int u, int flow) {
if(u == t) return flow;
int tot = 0;
for(int &i = cur[u] ; i&&tot < flow ; i = e[i].nxt) {
int v = e[i].to;
if(dep[v] == dep[u] + 1&&e[i].flow < e[i].cap) {
if(int t = dfs(v, min(e[i].cap - e[i].flow, flow - tot))) {
tot += t;
e[i].flow += t;
e[((i-1) ^ 1 ) + 1].flow -= t;
}
}
}
return tot;
}
void dinic(){
while(bfs()) {
memcpy(cur,head,t*4+4);
MaxFlow+=dfs(s, 0x3f3f3f3f);
}
}
int main() {
freopen("in.txt", "r", stdin);
n=init();
for(int i = 1 ; i <= n ; i++) {
for(int j = 1; j <= n; j++) {
b[i][j] = init();
ans += b[i][j];
}
}
for(int i = 1;i <= n;i++){
c[i]=init();
}
s=0;t=n*n+n+1;
for(int i = 1;i <= n ;i++){
for(int j = 1;j <= n;j++){
adde(s,i*n+j-n,b[i][j]);
adde(i*n+j-n,s,0);
}
}
for(int i = 1;i <= n;i++) {
adde(n*n+i, t, c[i]);
adde(t, n*n+i, 0);
}
int kkk=n*n;
for(int i = 1;i <= n ;i++) {
int ttt=i*n-n;
for(int j=1 ;j <=n ;j++) {
adde(ttt+j,kkk+i,inf);
adde(kkk+i,ttt+j,0);
adde(ttt+j,kkk+j,inf);
adde(kkk+j,ttt+j,0);
}
}
dinic();
cout<<ans-MaxFlow<<endl;
fclose(stdin);
return 0;
}
洛谷 [P3973] 线性代数的更多相关文章
- 洛谷P3973 - [TJOI2015]线性代数
Portal Description 给定一个\(n\times n\)的矩阵\(B\)和一个\(1×n\)的矩阵\(C\).求出一个\(1×n\)的01矩阵\(A\),使得\(D=(A×B-C)×A ...
- 【洛谷P3973】[TJOI2015]线性代数(最小割)
洛谷 题意: 给出一个\(n*n\)的矩阵\(B\),再给出一个\(1*n\)的矩阵\(C\). 求一个\(1*n\)的\(01\)矩阵\(A\),使得\(D=(A\cdot B-C)\cdot A^ ...
- 洛谷1640 bzoj1854游戏 匈牙利就是又短又快
bzoj炸了,靠离线版题目做了两道(过过样例什么的还是轻松的)但是交不了,正巧洛谷有个"大牛分站",就转回洛谷做题了 水题先行,一道傻逼匈牙利 其实本来的思路是搜索然后发现写出来类 ...
- 洛谷P1352 codevs1380 没有上司的舞会——S.B.S.
没有上司的舞会 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题目描述 Description Ural大学有N个职员,编号为1~N.他们有 ...
- 洛谷P1108 低价购买[DP | LIS方案数]
题目描述 “低价购买”这条建议是在奶牛股票市场取得成功的一半规则.要想被认为是伟大的投资者,你必须遵循以下的问题建议:“低价购买:再低价购买”.每次你购买一支股票,你必须用低于你上次购买它的价格购买它 ...
- 洛谷 P2701 [USACO5.3]巨大的牛棚Big Barn Label:二维数组前缀和 你够了 这次我用DP
题目背景 (USACO 5.3.4) 题目描述 农夫约翰想要在他的正方形农场上建造一座正方形大牛棚.他讨厌在他的农场中砍树,想找一个能够让他在空旷无树的地方修建牛棚的地方.我们假定,他的农场划分成 N ...
- 洛谷P1710 地铁涨价
P1710 地铁涨价 51通过 339提交 题目提供者洛谷OnlineJudge 标签O2优化云端评测2 难度提高+/省选- 提交 讨论 题解 最新讨论 求教:为什么只有40分 数组大小一定要开够 ...
- 洛谷P1371 NOI元丹
P1371 NOI元丹 71通过 394提交 题目提供者洛谷OnlineJudge 标签云端评测 难度普及/提高- 提交 讨论 题解 最新讨论 我觉得不需要讨论O long long 不够 没有取 ...
- 洛谷P1538迎春舞会之数字舞蹈
题目背景 HNSDFZ的同学们为了庆祝春节,准备排练一场舞会. 题目描述 在越来越讲究合作的时代,人们注意的更多的不是个人物的舞姿,而是集体的排列. 为了配合每年的倒计时,同学们决定排出——“数字舞蹈 ...
随机推荐
- c语言中标识符的作用域
1.代码块作用域(block scope) 位于一对花括号之间的所有语句称为一个代码块,在代码块的开始位置声明的标识符具有代码块作用域,表示它们可以被这个代码中的所有语句访问.函数定义的形式参数在函数 ...
- spring中用到的设计模式
http://www.cnblogs.com/pengmengnan/p/6717766.html 一 : 工厂模式工厂模式主要是为创建对象提供过度接口,以便将创建对象的具体 过程屏蔽隔离起来,达到提 ...
- tinyxml的封装与使用(转载)
tinyxml是个高效精简的xml解析开源代码. 针对tinyxml直接使用对于对xml不是很熟悉的入门新手来说,有些概念难以理解,因此我将其封装后,供大家使用. 头文件: #include &quo ...
- C#面试题整理(1)
最近在看CLR VIA C#,发现了一些案例很适合来做面试题.特此整理: 1,System.Object里的GetType方法是否为虚函数?说出理由. 答案:不是,因为C#是一种类型安全的语言,如果覆 ...
- seaJS 模块加载过程分析
先看一个seajs的官方example, 以下以seajs.use('main')为例, 解析加载mod main的过程 //app.html seajs.use("main") ...
- webzip怎么用 如何用webzip下载整个网站?
相信很多站长对webzip这款软件都并不感到陌生,它功能强大,能够完整下载网站的内容,或者你也可以选择自行设置下载的层数.文件类型.网页与媒体文件的定位等等.具体详情你可以在百度上去搜一下.由于web ...
- Struts2 05---拦截器
一.Struts2拦截器原理: Struts2拦截器的实现原理相对简单,使用了aop思想和责任链模式,当请求struts2的action时,Struts 2会查找配置文件,并根据其配置实例化相对的拦截 ...
- dedecms环境优化
路径:dedecms/dede/templates/index_body.htm <script type="text/javascript">function sho ...
- Redis的事务功能详解
Redis的事务功能详解 MULTI.EXEC.DISCARD和WATCH命令是Redis事务功能的基础.Redis事务允许在一次单独的步骤中执行一组命令,并且可以保证如下两个重要事项: >Re ...
- eclipse导入web项目变成java项目解决办法
右键工程,properties-> Project Facets-> 点convert to faceted..连接 -> 把Dynamic Web Moudle勾上