Assignment 【HDU - 2853】【KM算法+思维】

题目链接

---

  这个思维，当真是有点想不到，一开始净是往link_x[]这个角度去想了，但是这样的话，又会影响KM形成的交错树，显然不是这么解的。

  题意：有N个公司，还有M个任务，任务数大于等于公司数，然后每家公司做每个任务都有对应的贡献，于是乎，初始在加上每个公司有原本应该做的任务，就有一个初始贡献，现在我们想要最少的改变公司和对应的做的任务，想知道最大可以改变的值是多少？

  思路：首先，改KM板子的话，我自认为没这个实力，然后的话，这题的突破口在于是否要修改这个问题上面。然后，初略看了他人的解法主要是有个放大。为什么要放大？我想了想，一共我们要取N个任务，那么，如果给每个值都放大成K倍(K > N)的话，我们再考虑每个原来的值（假定给原来的每个初定任务都“+1”贡献），我们如果全部选择的是初始条件的话，它除以K之后的答案是不变的。并且，如果有换的话，证明一点，换的方向肯定是比它更优的，于是这个K倍就引刃而解了！

  说实话，这是个天才的游戏——K倍，真的很有趣。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <unordered_map>
#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
#define MP(a, b) make_pair(a, b)
#define Min3(a, b, c) min(a, min(b, c))
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 55;
int N, M, cop[maxN], K, old_val;
struct KM
{
    int n, mp[maxN][maxN], link_x[maxN], link_y[maxN];
    bool vis_x[maxN], vis_y[maxN];
    int que[maxN], top, fail, pre[maxN];
    ll hx[maxN], hy[maxN], slk[maxN];
    bool check(int i)
    {
        vis_x[i] = true;
        if(link_x[i])
        {
            que[fail++] = link_x[i];
            vis_y[link_x[i]] = true;
            return true;
        }
        while(i)
        {
            link_x[i] = pre[i];
            swap(i, link_y[pre[i]]);
        }
        return false;
    }
    void bfs(int S)
    {
        for(int i=1; i<=n; i++)
        {
            slk[i] = INF;
            vis_x[i] = vis_y[i] = false;
        }
        top = 0; fail = 1; que[0] = S;
        vis_y[S] = true;
        while(true)
        {
            ll d = 0;
            while(top < fail)
            {
                for(int i=1, j = que[top++]; i<=n; i++)
                {
                    if(!vis_x[i] && slk[i] >= (d = hx[i] + hy[j] - mp[i][j]))
                    {
                        pre[i] = j;
                        if(d) slk[i] = d;
                        else if(!check(i)) return;
                    }
                }
            }
            d = INF;
            for(int i=1; i<=n; i++)
            {
                if(!vis_x[i] && d > slk[i]) d = slk[i];
            }
            for(int i=1; i<=n; i++)
            {
                if(vis_x[i]) hx[i] += d;
                else slk[i] -= d;
                if(vis_y[i]) hy[i] -= d;
            }
            for(int i=1; i<=n; i++)
            {
                if(!vis_x[i] && !slk[i] && !check(i)) return;
            }
        }
    }
    void Init()
    {
        for(int i=1; i<=n; i++)
        {
            link_x[i] = link_y[i] = 0;
            hy[i] = slk[i] = 0;
            vis_y[i] = false;
        }
        for(int i=1; i<=n; i++)
        {
            hx[i] = 0;
            for(int j=1; j<=n; j++)
            {
                if(hx[i] < mp[i][j]) hx[i] = mp[i][j];
            }
        }
    }
    void solve()
    {
        n = M;
        Init();
        for(int i=1; i<=n; i++) bfs(i);
        ll ans = 0;
        for(int i=1; i<=n; i++) ans += hx[i] + hy[i];
        printf("%lld %lld\n", N - (ans % K), ans / K - old_val);
    }
}km;
int main()
{
    while(scanf("%d%d", &N, &M) != EOF)
    {
        old_val = 0;
        K = N + 1;  //放大倍数
        for(int i=1; i<=N; i++)
        {
            for(int j=1; j<=M; j++)
            {
                scanf("%d", &km.mp[i][j]);
                km.mp[i][j] *= K;
            }
        }
        for(int i = N + 1; i<=M; i++)
        {
            for(int j=1; j<=M; j++) km.mp[i][j] = 0;
        }
        for(int i=1, id; i<=N; i++)
        {
            scanf("%d", &id);
            old_val += km.mp[i][id] / K;
            km.mp[i][id] ++;
        }
        km.solve();
    }
    return 0;
}