https://www.luogu.com.cn/problem/AT_arc142_e
他的充要条件是是什么:
- a i , a j ≥ m i n ( b i , b j ) a_i,a_j\ge min(b_i,b_j) ai,aj≥min(bi,bj)
- 存在 a i ≥ m a x ( b i , b j ) a_i\ge max(b_i,b_j) ai≥max(bi,bj)
第一个条件直接预处理一下。
第二个条件,同时有非常多限制交叉,我们就非常想网络流建模来搞。(而且数据范围明细提示着我们)但是有个尴尬的问题,哪些点连源点,哪些连汇点。
然后就是此题最巧妙的地方了,可以构造二分图!
我们发现 ( a i , b i ) (a_i,b_i) (ai,bi) 是成对的,这似乎提示着我们什么。考虑对其分类
对于一对限制 ( i , j ) (i,j) (i,j),如果目前只满足条件1,不满足条件2 , a i ≥ b i a_i\ge b_i ai≥bi 和 a j ≥ b j a_j\ge b_j aj≥bj 必然是一个成立,一个不成立。同时,这里 i , j i,j i,j 独立,所以这就变成了一个二分图了!
最后是建模。我们这里发现是左边和右边至少满足一个,典型的最小割模型。
#include<bits/stdc++.h>
using namespace std;
#define int long long
inline int read(){int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;
ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<1)+
(x<<3)+(ch^48);ch=getchar();}return x*f;}
#define Z(x) (x)*(x)
#define pb push_back
//#define M
//#define mo
#define N 100010
using namespace atcoder;
struct node {
int u, v;
};
int n, m, i, j, k, T;
int ans, a[N], b[N], u, v, S;
vector<node>G;
signed main()
{
// freopen("in.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
// srand(time(NULL));
// T=read();
// while(T--) {
//
// }
mf_graph<int>F(200010);
n=read(); S=(n+5)*100+1; T=(n+5)*100+2;
for(i=1; i<=n; ++i) a[i]=read(), b[i]=read();
m=read();
for(i=1; i<=m; ++i) {
u=read(); v=read();
k=min(b[u], b[v]);
if(a[u]<k) ans+=k-a[u], a[u]=k;
if(a[v]<k) ans+=k-a[v], a[v]=k;
G.pb({u, v});
}
// printf("%lld\n", ans);
for(i=1; i<=n; ++i) {
// printf("%lld %lld\n", a[i], b[i]);
if(a[i]>=b[i]) {
// printf("> %lld\n", i);
for(j=a[i]+1; j<=100; ++j)
F.add_edge(S, i*100+j, 1);
for(j=1; j<100; ++j)
F.add_edge(i*100+j, i*100+j+1, 1e9);
}
else {
for(j=a[i]+1; j<=100; ++j)
F.add_edge(i*100+j, T, 1);
for(j=1; j<100; ++j)
F.add_edge(i*100+j+1, i*100+j, 1e9);
}
}
for(auto t : G) {
u=t.u; v=t.v;
k=max(b[u], b[v]);
if(a[u]>=b[u] && a[v]>=b[v]) continue;
if(a[u]>=b[v] && a[v]>=b[u]) continue;
// if(k==13) continue;
if(a[u]>=b[u]) {
// printf("%d -> %d\n", u*100+k, v*100+k);
F.add_edge(u*100+k, v*100+k, 1e9);
}
else {
// printf("%d -> %d\n", v*100+k, u*100+k);
F.add_edge(v*100+k, u*100+k, 1e9);
}
}
ans=ans+F.flow(S, T);
printf("%lld", ans);
return 0;
}