Asignación de asiento D ZOJ-3997 (flujo máximo)

Título:
Hay 10 tipos de personas, $yo\; la$ gente solo quiere hacer$Tengo$ varios asientos numerados. Todo tipo de persona tiene$a\left[i\right]$ uds.
Un total de$m$ asientos, ¿cuántas personas pueden sentarse?

Idea:
un número puede ser un múltiplo de varios números del 1 al 10, así que si pones $Si\; es$ un múltiplo de$i$ , también puede afectar a otros múltiplos. Esto tiene que ponderarse.

Primero averigüe el múltiplo común de cualquier número del 1 al 10, después de eliminar la duplicación, quedan 48. Supuesto $i\; los$ números son$num\left[i\right]$， 则num [i] num [i]El número de$n$$u$$m$$[$$i$$]$ es$cnt\left[i\right]=\frac{m}{num\left[i\right]}$。

---

A continuación, tenemos que eliminar los duplicados (equivalente a la exclusión), recorrer num num una vez de atrás hacia $num$ matriz, para$num\left[i\right]$ este número, encuentra$El\; número\; de\; la$ matriz$n$$u$$m$ que es mayor y un múltiplo de él$num\left[j\right]$ , luego$cnt\left[i\right]=cnt\left[i\right]-cnt\left[j\right]$。

Entonces, para cada número $num\left[i\right]$ Asigne tantos asientos como pueda.

---

El siguiente paso es construir un mapa.
Establecer el punto de origen en 0 y el punto de sumidero en 100. Luego, el punto de origen construye un mapa para cada tipo de persona, y la capacidad es el número de tales personas $A\left[I\right]$ ,
cada uno de los cualesnum numLos múltiplos de la matriz$n$$u$$m$ se utilizan para crear un gráfico, la capacidad es$cnt\left[i\right]$，
$num$ matriz para construir un mapa para cada dirección numérica al punto de hundimiento, la capacidad es$cnt\left[i\right]$ , el
resultado es el caudal máximo.

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
typedef long long ll;

//网络流开始
const int inf = 1 << 30, N = 2005, M = 2005;
int head[N], ver[M], edge[M], Next[M], d[N];
int n, m, s, t, tot;
int maxflow;
queue<int> q;
void add(int x, int y, int z) {
    
    
    ver[++tot] = y, edge[tot] = z, Next[tot] = head[x], 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; 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, int flow) {
    
     // 在当前分层图上增广
    if (x == t) return flow;
    int rest = flow, k;
    for (int i = head[x]; i && rest; i = Next[i])
        if (edge[i] && d[ver[i]] == d[x] + 1) {
    
    
            k = dinic(ver[i], min(rest, edge[i]));
            if (!k) d[ver[i]] = 0; // 剪枝，去掉增广完毕的点
            edge[i] -= k;
            edge[i ^ 1] += k;
            rest -= k;
        }
    return flow - rest;
}


int cal() {
    
    
    int flow = 0;
    while (bfs())
        while (flow = dinic(s, inf)) maxflow += flow;
    return maxflow;
}

//网络流结束

int num[2005],cnt;

int gcd(int n,int m) {
    
    
    return m == 0 ? n : gcd(m,n % m);
}

void dfs(int x,int now) {
    
    
    if(x == 11) {
    
    
        num[++cnt] = now;
        return;
    }
    int tmp = gcd(now,x);
    dfs(x + 1,now * x / tmp); //算公倍数
    dfs(x + 1,now); //不算
}

void Pre() {
    
    
    dfs(1,1);
    sort(num + 1,num + 1 + cnt);
    cnt = unique(num + 1,num + 1 + cnt) - (num + 1);
}

void init() {
    
    
    s = 0;t = 100;
    tot = 1;
    maxflow = 0;
    memset(head,0,sizeof(head));
}

int cnt_p[20]; //输入的1~10
int cnt_n[N]; //每个数字对应的个数

int main() {
    
    
    Pre();
    int T;scanf("%d",&T);
    while(T--) {
    
    
        scanf("%d",&n);
        init();
        for(int i = 1;i <= 10;i++) {
    
    
            scanf("%d",&cnt_p[i]);
        }
        for(int i = 1;i <= cnt;i++) {
    
    
            cnt_n[i] = n / num[i];
        }
        for(int i = cnt;i >= 1;i--) {
    
    
            for(int j = i + 1;j <= cnt;j++) {
    
    
                if(num[j] % num[i] == 0) {
    
    
                    cnt_n[i] -= cnt_n[j];
                }
            }
            if(cnt_n[i] != 0) {
    
    
                add(i + 10,100,cnt_n[i]);
            }
        }
        for(int i = 1;i <= 10;i++) {
    
    
            if(cnt_p[i] == 0) continue;
            add(0,i,cnt_p[i]);
            for(int j = 1;j <= cnt;j++) {
    
    
                if(num[j] > n) break;
                if(num[j] % i == 0 && cnt_n[j] != 0) {
    
    
                    add(i,j + 10,cnt_n[j]);
                }
            }
        }
        printf("%d\n",cal());
    }
    return 0;
}