bzoj3308: 九月的咖啡店 数学 玄学 最大流

bzoj3308: 九月的咖啡店

Description

深绘里在九份开了一家咖啡让，如何调配咖啡民了她每天的头等大事
我们假设她有N种原料，第i种原料编号为i,调配一杯咖啡则需要在这
里若干种兑在一起。不过有些原料不能同时在一杯中，如果两个编号
为i,j的原料，当且仅当i与j互质时，才能兑在同一杯中。
现在想知道，如果用这N种原料来调同一杯咖啡，使用的原料编号之和
最大可为多少。

Input

一个数字N

Output

如题

Sample Input

10

Sample Output

30

HINT

1<=N<=200000

分析

网络流&数学神题。
两个结论扔出来的时候已经懵逼了。
对于最后的最大编号之和的数的选择:
1.至多2个质因数存在于一个数中
2.一个大于$\sqrt n$，一个小于$\sqrt n$
不会证不会证不会证。
大概的思路就是如果单独放质因数应该会很优，但是举举例子会发现也有可能存在两个质因数更优的情况。
比如hzwer大神的：
3，5，n=20
9+5<3*5
然后就是费用流匹配小于$\sqrt n$和大于$\sqrt n$的数。
但是会炸。
于是优化，考虑单独选数对答案的贡献，再考虑把两个单独的数修改成匹配状态的收益$V_{ab}-V_{a}-V_{b}$
然后不断spfa增广直到费用小于零。
虚渊玄

代码

/**************************************************************
    Problem: 3308
    User: 2014lvzelong
    Language: C++
    Result: Accepted
    Time:140 ms
    Memory:13540 kb
****************************************************************/

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
const int N = 220000;
const int M = 440000;
const int inf = 1e9;
int read() {
    char ch = getchar();int x = 0, f = 1;
    while(ch < '0' || ch > '9') {if(ch == '-') f = -1; ch = getchar();}
    while(ch >= '0' && ch <= '9') {x = x * 10 + ch - '0'; ch = getchar();}
    return x * f; 
}
int pre[N], head[N], to[M], w[M], f[M], nxt[M], dis[N], Q[M], p[N], top, tot, n, m, ans, s, t;
bool vis[N];
void add(int u, int v, int ww, int ff) {
    to[++top] = v; w[top] = ww; nxt[top] = head[u]; head[u] = top; f[top] = ff;
}
void adds(int u, int v, int w, int f) {add(u, v, w, f); add(v, u, 0, -f);}

void pri(int n) {
    vis[1] = 1;
    for(int i = 2;i <= n; ++i) {
        if(!vis[i]) p[++tot] = i;
        for(int j = 1;j <= tot && i * p[j] <= n; ++j) {
            vis[i * p[j]] = true;
            if(!(i % p[j])) break;
        }
    }
}

int calc(int n, int x) {int t; for(t = x; 1LL * t * x <= n; t = t * x) ;return t;}

void Build() {
    s = tot + 1; t = s + 1; int pos = 0;
    for(int i = 1;i <= tot; ++i) {
        if(p[i] >= n / 2) {ans += p[i]; continue;}
        if(1LL * p[i] * p[i] <= n) adds(s, i, 1, 0), ans += calc(n, p[i]);
        else pos = (!pos) ? i : pos, adds(i, t, 1, 0), ans += p[i];
    }
    for(int i = 1;i < pos; ++i)
        for(int j = pos; j <= tot; ++j) {
            if(1LL * p[i] * p[j] > n) break;
            int f = calc(n / p[j], p[i]) * p[j] - calc(n, p[i]) - p[j];
            if(f > 0) adds(i, j, 1, f);
        }
}

bool spfa() {
    for(int i = 0;i <= t; ++i) dis[i] = -inf, vis[i] = 0;
    int L = 0, R; Q[R = 1] = s; dis[s] = 0; vis[s] = true;
    while(L < R) {
        int u = Q[++L]; vis[u] = false;
        for(int i = head[u]; i; i = nxt[i]) 
            if(w[i] && dis[to[i]] < dis[u] + f[i]) {
                dis[to[i]] = dis[u] + f[i];
                pre[to[i]] = i;
                if(!vis[to[i]]) {
                    vis[to[i]] = true;
                    Q[++R] = to[i];
                }
            }
    }
    return dis[t] >= 0;
}

int augment() {
    int flow = inf;
    for(int i = t; i != s; i = to[pre[i] ^ 1]) flow = min(flow, w[pre[i]]);
    for(int i = t; i != s; i = to[pre[i] ^ 1]) w[pre[i]] -= flow, w[pre[i] ^ 1] += flow;
    return flow * dis[t];
}

void MCMF() {
    int mincost;
    while(spfa())
        mincost += augment();
    ans += mincost;
}

int main() {
    n = read(); 
    pri(n); top = 1;
    Build();
    MCMF();
    printf("%d\n", ans + 1);
    return 0;
}