sgu 210


KM ,把王子放在二分图 X 侧,妹子放在Y侧,边权为 (Ai)2 KM 即可,
具体实现时我把王子放在二分图 Y 侧,妹子放在X侧,这样就更容易输出答案了。


#include<cstdio>
#include<cstdlib>
#include<ctime>
#include<cmath>
#include<cstring>
#include<map>
#include<string>
#include<stack>
#include<queue>
#include<utility>
#include<iostream>
#include<algorithm>

const int maxn = 405, INF = 0x3f3f3f3f;

int n, a[maxn];
int w[maxn][maxn];
bool visx[maxn], visy[maxn];
int slack[maxn], lx[maxn], ly[maxn];
int linky[maxn];

void init()
{
    std::cin >> n;
    for(int i = 1; i <= n; i++)
        std::cin >> a[i];
    for(int i = 1, num; i <= n; i++)
    {
        std::cin >> num;
        for(int j = 1, to; j <= num; j++)
            std::cin >> to, w[to][i] = a[i]*a[i];
    }
}
bool find(int x)
{
    visx[x] = true;

    for(int y = 1; y <= n; y++)
        if(!visy[y])
        {
            int dec = lx[x] + ly[y] - w[x][y];

            if(!dec)
            {
                visy[y] = true;
                if(!linky[y] || find(linky[y]))
                {
                    linky[y] = x;
                    return true;
                }
            }
            slack[y] = std::min(slack[y], dec);
        }
    return false;   
}
void KM()
{
    memset(lx,0,sizeof(lx));
    memset(ly,0,sizeof(ly));
    memset(linky,0,sizeof(linky));

    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= n; j++)
            lx[i] = std::max(lx[i],w[i][j]);

    for(int i = 1; i <= n; i++)
    {
        memset(slack,INF,sizeof(slack));

        while(true)
        {
            memset(visx,false,sizeof(visx));
            memset(visy,false,sizeof(visy));

            if(find(i)) break;

            int d = INF;

            for(int i = 1; i <= n; i++)
                if(!visy[i]) d = std::min(slack[i], d);
            for(int i = 1; i <= n; i++)
            {
                if(visx[i]) lx[i] -= d; 
                if(visy[i]) ly[i] += d;
                else slack[i] -= d;
            }
        }
    }
}
void prt()
{
    for(int i = 1; i <= n; i++)
        std::cout << (w[linky[i]][i]?linky[i]:0) << ' ';
}


int main()
{
#ifndef ONLINE_JUDGE
    freopen("sgu210.in","r",stdin);
    freopen("sgu210.out","w",stdout);
#endif

    init(), KM(), prt();

#ifndef ONLINE_JUDGE
    fclose(stdin);
    fclose(stdout);
#endif
    return 0;
}

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值