这就很离谱
查看原帖
这就很离谱
169928
YRZ001030楼主2020/9/29 15:47

以下代码的主函数建图代码,我通过copy别人的网络流跑出来了,主函数建图代码是正确的; 我的网络流最大流板子是蓝书上面的,并且经过落谷上的最大流模板检测也是正确的;可是为啥一结合就错呢?

///链接:https://www.luogu.com.cn/problem/P3355
/*题意: 一个n*n的棋盘,马走日的规则互相攻击,棋盘有一些点不可放置马,问在都不互相攻击的前提下,能放多少马。
  题解:有点类似网格取数的最小割。
*/

#include"stdio.h"
#include"string.h"
#include"stack"
#include"map"
#include"math.h"
#include"iostream"
#include"vector"
#include"queue"
#include"algorithm"
using namespace std;
#define OK printf("\n");
#define Debug printf("this_ok\n");
#define INF 1e18
typedef long long ll;
#define scanll(a,b) scanf("%lld%lld",&a,&b);
#define scanl(a) scanf("%lld",&a);
#define printl(a,b) if(b == 0) printf("%lld ",a); else printf("%lld\n",a);
#define print_int(a,b) if(b == 0) printf("%d ",a); else printf("%d\n",a);
typedef pair<int,int> PII;

inline int read(){
    int s = 0, w = 1; char ch = getchar();
    while(ch < '0' || ch > '9')   { if(ch == '-') w = -1; ch = getchar(); }
    while(ch >= '0' && ch <= '9') { s = (s << 3) + (s << 1) + (ch ^ 48); ch = getchar(); }
    return s * w;
}
const ll mod = 998244353;
const int N = 500010,M = 500010;
const  double pi = acos(-1);
const int inf = 1 << 29;
const int dirx[9] = {-2,-2,2,2,1,1,-1,-1};
const int diry[9] = {-1,1,-1,1,-2,2,-2,2};
int n,m,t,s,tot;
int maxflow,sum;
int head[N],ver[M],Next[M],edge[M],d[M];
queue<int> q;
int id[210],tot1;
int vis[210],to[210];

void add(int x,int y,int z){
    ver[++ tot] = y; Next[tot] = head[x];  edge[tot] = z; head[x] = tot;
    ver[++ tot] = x; edge[tot] = 0; Next[tot] = head[y]; head[y] = tot;
}

bool bfs(){
    memset(d,0,sizeof(d));
    while(q.size())q.pop();
    q.push(s); d[s] = 1;
    while(q.size()){
        int x = q.front(); q.pop();
        for(int i = head[x]; i!=-1; i = Next[i])
        if(edge[i] && !d[ver[i]]){
            q.push(ver[i]); d[ver[i]] = d[x] + 1;
            if(ver[i] == t) return 1;
        }
    }
    return 0;
}

int dinic(int x,ll flow){
    if(x == t) return flow;
    ll rest = flow,k;
    for(int i = head[x]; i!=-1 && rest; i = Next[i]){
         if(edge[i] && d[ver[i]] == d[x] + 1){
            k = dinic(ver[i],min(rest,(ll)edge[i]));
            if(!k) d[ver[i]] = 0;
            edge[i] -= k;
            edge[i ^ 1] += k;
            rest -= k;
         }
    }
    return flow - rest;
}
int main(){
    n = read(),s = 0; t = 50003;
    for(int i = 0; i <= 1010; i ++) head[i] = -1;
    tot = 1; tot1=  0;
    int flow = 0; int num = 0;
    for(int i = 1; ; i ++){
       add(s,i * 2,1); add(i * 2 + 1,t,1);
       /*for(int j = 1; j < i; j ++){
          int sum = i + j;
          if(sqrt(sum) * sqrt(sum) == sum){
            add(j * 2,i * 2 + 1,1);
          }
       }*/
       num = i;
        for(int j=sqrt(num)+1;j*j<2*num;++j)
            add((j*j-num)<<1,(num<<1)|1,1);
       flow = 0;
       while(bfs())
         while(flow = dinic(s,inf)) maxflow += flow;
       if(!flow) {
          id[++ tot1] = i;
       }
       if(tot1 > n) {
         num = i;    break;
       }
    }
    printf("%d\n",num);
    for(int i = 1; i < num; i ++){
        for(int j = head[i * 2]; j;j = Next[j]){
            if(edge[j] == 0){
                to[i] = ver[j] / 2; break;
            }
        }
    }
    for(int i = 1; i <= n; i ++){
        int k = id[i];
        while(to[k] != 0){
            printf("%d ",k); k = to[k];
        }
        printf("%d\n",k);
    }
}
/*
3
1 2 3
2
2 6
*/


2020/9/29 15:47
加载中...