最小割小记

参考博客：最小割浅谈

关于最小割

1. 常用描述
   表述一：删去若干条边使得源点到汇点不连通，求删边的权值和的最小可能值。
   表述二：将点集分为$(S,T)$，记所有从$S$中出发到$T$中的边的权值和为$c(S,T)$，求$c(S,T)$的最小值。

2. 求最小割
   a. 以权值为容量，该网络最大流的值即为最小割的值
   b. 在残量网络中，从源点出发进行一次增广BFS，即得到一个分割。该分割是一个最小割。

题型1：对求最小割原理的理解

[AHOI2009]最小割

摘自此Blog

[SDOI2014]LIS

摘自此Blog

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#define re register
#define maxn 1505
#define LL long long
#define inf 999999999
inline int max(int a,int b){
   return a>b?a:b;}
inline int min(int a,int b){
   return a<b?a:b;}
inline int read()
{
   
    char c=getchar();int x=0;while(c<'0'||c>'9') c=getchar();
    while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
struct node{
   int a,b,c,rk;}g[maxn];
inline int cmp(node A,node B) {
   return A.c<B.c;}
struct E{
   int v,nxt,w,f;}e[maxn*maxn];
int n,num=1,S,T,Test;
int ans[maxn],cnt,id[maxn],vis[maxn];
int d[maxn],dp[maxn],head[maxn],cur[maxn],in[maxn],out[maxn];
inline void add(int x,int y,int z) {
   e[++num].v=y;e[num].nxt=head[x];head[x]=num;e[num].w=z;}
inline void C(int x,int y,int z) {
   add(x,y,z),add(y,x,0);}
inline int check(int s,int t)
{
   
    std::queue<int> q;
    memset(vis,0,sizeof(vis));
    q.push(s);vis[s]=1;
    while(!q.empty())
    {
   
        int k=q.front();q.pop();
        if(k==t) return 1;
        for(re int i=head[k];i;i=e[i].nxt)
        if(!vis[e[i].v]&&e[i].w>e[i].f) q.push(e[i].v),vis[e[i].v]=1;
    }
    return 0;
}
inline int BFS(int s,int t)
{
   
    std::queue<int> q;
    memcpy(cur,head,sizeof(head));
    memset(d,0,sizeof(d));
    d[s]=1,q.push(s);
    while(!q.empty())
    {
   
        int k=q.front();q.pop();
        for(re int i=head[k];i;i=e[i].nxt)
        if(!d[e[i].v]&&e[i].w>e[i].f) 
        {
   
            d[e[i].v]=d[k]+1;
            if(e[i].v==t) return 1;
            q.push(e[i].v);
        }
    }
    return d[t];
}
int dfs(int x,int now,int t)
{
   
    if(x==t||!now) return now;
    int flow=0,ff;
    for(re int& i=cur[x];i;i=e[i].nxt)
    if(d[e[i].v]==d[x]+1)
    {
   
        ff=dfs(e[i].v,min(now,e[i].w-e[i].f),t);
        if(ff<=0) continue;
        now-=ff,flow+=ff;
        e[i].f+=ff,e[i^1].f-=ff;
        if(!now) break;
    }
    return flow;
}
int main()
{
   
    Test=read();
    while(Test--)
    {
   
        n=read();cnt=0