SPJ有问题趴
查看原帖
SPJ有问题趴
291604
王茗仟楼主2022/11/6 20:52

先给代码

#define int long long 
#define re  register
#define in inline
#define max(a,b) ((a)>(b)?a:b)
#define min(a,b) ((a)<(b)?a:b)
using namespace std;
const int MAXN=1e5+10;
const int MAXM=2e6+19;
const int INF=0x3f3f3f3f3f;

in int read(){
    re int x=0,f=0;re char c=getchar();
    while(!isdigit(c)) f|=(c=='-'),c=getchar();
    while(isdigit(c)) x=(x<<3)+(x<<1)+c-'0',c=getchar();
    return f?-x:x;
}

int n,m,s,t,dis[MAXN],cur[MAXN];

struct edge{
    int u,v,w;
    int nx;
}e[MAXM];
int head[MAXN],tot;

in void  add(int u,int v,int w){
    e[tot].u=u;
    e[tot].v=v;
    e[tot].w=w;
    e[tot].nx=head[u];
    head[u]=tot++;
}

bool bfs(){
//广搜一下测一下剩下点的深度,顺便找出增广路
    queue<int>q;
    memset(dis,-1,sizeof(dis));
    q.push(s);
    dis[s]=0;
    cur[s]=head[s];//cur是当前弧优化的宝贝
    while(q.size()){
        int u=q.front();
        q.pop();
        for(int i=head[u];~i;i=e[i].nx){
            int v=e[i].v;
            cur[v]=head[v];
            if(dis[v]==-1&&e[i].w){
                dis[v]=dis[u]+1;
                if(v==t){
                    return 1;
                }
                q.push(v);
            }
        }
    }
    return 0;//不存在增广路了,可以结束le
}

int dfs(int u,int lim){
// 当前节点     能增加的最大流量
//深搜??也就是增广一下,给下流
    if(u==t){
        return lim;
    }

    int flow=0;//从u点分出去的
    for(int i=cur[u];~i&&flow<lim;i=e[i].nx){
        cur[u]=i;//当前弧优化,很新对你来说,(走过就不再走了)
        int v=e[i].v;
        if(dis[v]==dis[u]+1&&e[i].w){
            int minf=dfs(v,min(e[i].w,lim-flow));//求要给v多少
//limit-flow给这条边能分的最多有多少
            e[i].w-=minf;
            e[i^1].w+=minf;
            flow+=minf;
            if(flow==lim){
            //全分完了直接结束吧
                return flow;
            }
        }   
    }
    return flow;

}

int dinic(){
    int ans=0;
    while(bfs()){
  //    cout<<ans<<endl;
        ans+=dfs(s,INF);
    }
    return ans;
}

signed main(){
    m=read();n=read();
    memset(head,-1,sizeof(head));
    
    s=0;t=n+1;
    for(re int i=1;i>0;i++){
        int u,v,w;
        u=read();v=read();
        if(u==-1&&v==-1) break;
        add(u,v,INF);add(v,u,0);
    }
    for(int i=1;i<=m;i++){
    	add(s,i,1);add(i,s,0);
	}
	for(int i=m+1;i<=n;i++){
		add(i,t,1);add(t,i,0);
	}
    
    cout<<dinic()<<endl;
    
    for(int u=1;u<=m;u++){
    	for(int i=head[u];i;i=e[i].nx){
    		int v=e[i].v;
    		if(v!=s&&v!=t&&e[i].w!=0&&e[i^1].w!=0){
    			cout<<u<<" "<<v<<endl;
			}
		}
	}
    
    return 0;
}

这是样例11的输入

10 55
19 52
20 44
5 55
14 31
29 35
20 37
13 37
4 45
28 34
13 53
21 44
11 30
15 53
3 37
20 50
20 41
10 46
11 36
2 38
23 42
20 42
29 32
11 38
16 50
1 32
2 55
2 42
11 50
7 39
23 47
22 33
10 43
5 40
4 55
7 31
2 52
2 36
11 33
17 38
29 36
7 54
26 52
26 35
5 39
15 53
12 32
9 38
25 43
27 55
16 36
14 51
19 39
21 35
22 44
24 51
16 42
8 52
8 55
19 56
1 53
25 34
4 47
27 37
5 50
21 32
15 45
15 42
27 53
5 35
13 51
27 40
14 44
8 55
9 48
13 44
19 40
14 37
15 51
14 46
9 47
28 33
17 52
19 41
2 50
5 48
3 55
1 38
28 47
18 49
22 49
13 42
5 53
13 34
26 36
23 33
29 39
14 44
26 54
16 33
19 52
1 34
11 32
21 40
9 52
12 39
6 38
22 56
29 42
27 32
15 44
13 39
4 35
20 44
3 42
15 55
10 51
6 53
25 40
22 37
14 49
13 49
18 42
6 42
27 47
3 31
27 43
23 40
5 55
22 55
18 43
-1 -1

这是样例11的输出

1 34
2 42
3 37
4 45
5 35
6 38
7 54
8 55
9 48
10 46
11 30
14 31
15 44
16 50
17 52
18 49
19 39
20 41
21 40
22 56
23 47
24 51
25 43
26 36
27 53
28 33
29 32

这是我的输出

1 34
2 50
3 37
5 35
6 53
7 31
8 52
9 48
11 30
13 39
14 46
15 45
16 36
17 38
18 49
19 41
20 44
21 32
22 56
23 33
24 51
25 40
26 54
27 43
28 47
29 42

我们两个4和13那里有差别,请大佬来看看,我感觉我的也是对的

2022/11/6 20:52
加载中...