萌新求助网络流 73 pts
查看原帖
萌新求助网络流 73 pts
220824
yyz1005楼主2024/1/19 11:06

具体的想法是把图搞成下面这样

女生那一边对称过去

从小到大枚举答案,每次给所有 L[i] 连一条流量为 11 的边到超级源点,给所有 R[i] 连一条流量为 11 的边到超级汇点,看新增的流量是否为 nn。

不知道为什么 Wa 了 3 个,都是大了一点

#include<bits/stdc++.h>
/*
#pragma GCC optimize(3)
#pragma GCC optimize("Ofast,no-stack-protector")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("inline")
#pragma GCC optimize("-fgcse")
#pragma GCC optimize("-fgcse-lm")
#pragma GCC optimize("-fipa-sra")
#pragma GCC optimize("-ftree-pre")
#pragma GCC optimize("-ftree-vrp")
#pragma GCC optimize("-fpeephole2")
#pragma GCC optimize("-ffast-math")
#pragma GCC optimize("-fsched-spec")
#pragma GCC optimize("unroll-loops")
#pragma GCC optimize("-falign-jumps")
#pragma GCC optimize("-falign-loops")
#pragma GCC optimize("-falign-labels")
#pragma GCC optimize("-fdevirtualize")
#pragma GCC optimize("-fcaller-saves")
#pragma GCC optimize("-fcrossjumping")
#pragma GCC optimize("-fthread-jumps")
#pragma GCC optimize("-funroll-loops")
#pragma GCC optimize("-fwhole-program")
#pragma GCC optimize("-freorder-blocks")
#pragma GCC optimize("-fschedule-insns")
#pragma GCC optimize("inline-functions")
#pragma GCC optimize("-ftree-tail-merge")
#pragma GCC optimize("-fschedule-insns2")
#pragma GCC optimize("-fstrict-aliasing")
#pragma GCC optimize("-fstrict-overflow")
#pragma GCC optimize("-falign-functions")
#pragma GCC optimize("-fcse-skip-blocks")
#pragma GCC optimize("-fcse-follow-jumps")
#pragma GCC optimize("-fsched-interblock")
#pragma GCC optimize("-fpartial-inlining")
#pragma GCC optimize("no-stack-protector")
#pragma GCC optimize("-freorder-functions")
#pragma GCC optimize("-findirect-inlining")
#pragma GCC optimize("-fhoist-adjacent-loads")
#pragma GCC optimize("-frerun-cse-after-loop")
#pragma GCC optimize("inline-small-functions")
#pragma GCC optimize("-finline-small-functions")
#pragma GCC optimize("-ftree-switch-conversion")
#pragma GCC optimize("-foptimize-sibling-calls")
#pragma GCC optimize("-fexpensive-optimizations")
#pragma GCC optimize("-funsafe-loop-optimizations")
#pragma GCC optimize("inline-functions-called-once")
#pragma GCC optimize("-fdelete-null-pointer-checks")
#pragma GCC optimize("-fdelete-null-pointer-checks,inline-functions-called-once,-funsafe-loop-optimizations,-fexpensive-optimizations,-foptimize-sibling-calls,-ftree-switch-conversion,-finline-small-functions,inline-small-functions,-frerun-cse-after-loop,-fhoist-adjacent-loads,-findirect-inlining,-freorder-functions,no-stack-protector,-fpartial-inlining,-fsched-interblock,-fcse-follow-jumps,-fcse-skip-blocks,-falign-functions,-fstrict-overflow,-fstrict-aliasing,-fschedule-insns2,-ftree-tail-merge,inline-functions,-fschedule-insns,-freorder-blocks,-fwhole-program,-funroll-loops,-fthread-jumps,-fcrossjumping,-fcaller-saves,-fdevirtualize,-falign-labels,-falign-loops,-falign-jumps,unroll-loops,-fsched-spec,-ffast-math,Ofast,inline,-fgcse,-fgcse-lm,-fipa-sra,-ftree-pre,-ftree-vrp,-fpeephole2",3)
*/
using namespace std;
typedef long long ll;
namespace IO{
	char buf[1 << 20], *p1 = buf, *p2 = buf;
	#define getchar() (p1 == p2 && (p2 = (p1 = buf) + \
		fread(buf, 1, 1 << 20, stdin), p1 == p2) ? EOF : *p1++)
	inline ll read(){
		ll x = 0,h = 1;
		char ch = getchar();
		while(!(ch=='-'||('0'<=ch&&ch<='9'))) ch = getchar();
		if(ch=='-') h = -1;
		else x = ch-'0';
		ch = getchar();
		while('0'<=ch&&ch<='9'){
			x = x*10+ch-'0';
			ch = getchar();
		}
		return x*h;
	}
	inline ll digread(){
		char ch = getchar();
		while(!('0'<=ch&&ch<='9')) ch = getchar();
		return ch-'0';
	}
	inline char readCH(){
		char ch = getchar();
		while(!(ch=='Y'||ch=='N')) ch = getchar();
		return ch;
	}
	inline void print(ll x){
		if(x>9) print(x/10);
		putchar(x%10+'0');
	}
	inline void write(ll x){
		if(x<0) putchar('-'),x = -x;
		print(x);
	}
}
using namespace IO;
#define File(name)                     \
    freopen(#name ".in", "r", stdin);  \
    freopen(#name ".out", "w", stdout);
#define ERR assert(0)
#define RunTime ((double)clock()/CLOCKS_PER_SEC)
/*
有时候程序没反应就把 IO 里面重定义 getchar 删掉 
-std=c++11  -Wall -Wextra -Wl,--stack=102400000000 -fno-ms-extensions
找不出来哪错了要检查 RE 问题,比如数组开小了,越界,函数返回值
用 long long TLE 可以尝试换成 int
对 vector 进行修改时,要 for(auto &v : vec) v++。必须加 &
n 行 m 列的矩阵先从左到右在从上到下编号,第 i 行 j 列是 (i-1)*m+j
网络流建图正边反边一定连着建
*/
const ll N = 100010;
const ll INF = 1e16;
struct Edge{
    ll from,to,val,nxt;
} edge[N];
ll D = 0;
ll s,t;
ll head[N],cnt = 1;
ll dis[N],cur[N];
ll dfs(ll id,ll flw){
    if(id==t) return flw;
    ll ret = 0;
    for(ll i = cur[id]; i&&flw>0; i = edge[i].nxt){
        cur[id] = i;
        ll v = edge[i].to;
        if(edge[i].val>0&&dis[v]==dis[id]+1){
            ll diff = dfs(v,min(flw,edge[i].val));
            if(diff==0) dis[v] = INF;
            edge[i].val-=diff;
            edge[i^1].val+=diff;
            ret+=diff;
            flw-=diff;
        }
    }
    return ret;
}
bool bfs(ll id){
    for(ll i = 1; i <= D; i++) dis[i] = INF;
    queue<ll> q;
    q.push(id);
    dis[id] = 0;
    cur[id] = head[id];
    while(!q.empty()){
        ll fir = q.front();
        q.pop();
        for(ll i = head[fir]; i; i = edge[i].nxt){
            ll v = edge[i].to;
            if(edge[i].val>0&&dis[v]==INF){
                q.push(v);
                cur[v] = head[v];
                dis[v] = dis[fir]+1;
                if(v==t) return 1;
            }
        }
    }
    return 0;
}
ll Dinic(){
    ll ans = 0;
    while(bfs(s)) ans+=dfs(s,INF);
    return ans;
}
void AddEdge(ll u,ll v,ll w){
    cnt++;
    edge[cnt].from = u,edge[cnt].to = v,edge[cnt].val = w;
    edge[cnt].nxt = head[u],head[u] = cnt;
}
void add(ll u,ll v,ll w){
	AddEdge(u,v,w);
	AddEdge(v,u,0);
}
ll n,K;
ll L1[N],L2[N],L3[N],L4[N],L[N];
ll R1[N],R2[N],R3[N],R4[N],R[N];
int main(){
	n = read(),K = read();
	for(ll i = 1; i <= n; i++) L1[i] = ++D,R1[i] = ++D;
	for(ll i = 1; i <= n; i++) L2[i] = ++D,R2[i] = ++D;
	for(ll i = 1; i <= n; i++) L3[i] = ++D,R3[i] = ++D;
	for(ll i = 1; i <= n; i++) L4[i] = ++D,R4[i] = ++D;
	for(ll i = 1; i <= n; i++) L[i] = ++D,R[i] = ++D;
	for(ll i = 1; i <= n; i++) add(L[i],L1[i],INF),add(L[i],L3[i],INF);
	for(ll i = 1; i <= n; i++) add(R2[i],R[i],INF),add(R4[i],R[i],INF);
	for(ll i = 1; i <= n; i++) add(L1[i],L2[i],K),add(R1[i],R2[i],K);
	for(ll i = 1; i <= n; i++) add(L3[i],L4[i],INF),add(R3[i],R4[i],INF);
	for(ll i = 1; i <= n; i++){
		for(ll j = 1; j <= n; j++){
			bool x = (readCH()=='Y');
			if(x) add(L4[i],R3[i],1);
			else add(L2[i],R1[i],1);
		}
	}
	s = ++D,t = ++D;
	for(ll res = 1; ; res++){
		for(ll i = 1; i <= n; i++) add(s,L[i],1),add(R[i],t,1);
		ll ans = Dinic();
		if(ans!=n){
			printf("%lld\n",res-1);
			return 0;
		}
	}
	fprintf(stderr, "TIME:%.3lf s\n", RunTime);
	return 0;
}
/*

*/
2024/1/19 11:06
加载中...