zoukankan      html  css  js  c++  java
  • Codeforces 1107 F (二分图最大权匹配)

    传送门

    题意:

    nn个贷款,你在月初可以申请任意一个贷款。对于第ii个贷款,你可以获得aia_i元,但是你要在未来的kik_i个月底还bib_i元。现在每个贷款只能申请一次,问在某个时刻,你能够获得的最多的钱。

    题目分析:

    我们思考这样的一个问题,假设我们选取了第ii个贷款,并经过了jj个月,那么对于选取第ii个贷款这个行为,它对答案的贡献为:max(aibi×(min(ki,j1)),0)max(a_i-b_i imes(min(k_i,j-1)),0)而显然,对于每一个贷款,都会有一个这样的跟月份有关的方程。

    同时我们发现,每一个贷款都是是相互独立的;同时,对于每一个月份也是相互独立的(即它们之间没有相互的制约关系)。而唯一存在制约关系的是贷款与月份。因此,我们可以把贷款和月份的关系看作一张二分图。其中第ii个贷款到第jj个月份的权值即为max(aibi×(min(ki,j1)),0)max(a_i-b_i imes(min(k_i,j-1)),0)

    因此,要求出某个时刻的最大值,我们即求出这张二分图的最大权匹配即可。这个我们只需要用KM算法即可。时间复杂度:O(n3)mathcal{O}(n^3)

    这个题还可以用时间复杂度为O(n2)mathcal{O}(n^2)dpdp去解决,dpdp的方法有待更新dots

    代码:

    #include <bits/stdc++.h>
    #define maxn 505
    using namespace std;
    typedef long long ll;
    const ll INF=0x3f3f3f3f3f3f3f3f;
    ll g[maxn][maxn];
    int a[maxn],b[maxn],k[maxn];
    int n;
    int nx,ny;
    ll linker[maxn],lx[maxn],ly[maxn];
    ll slack[maxn];
    bool visx[maxn],visy[maxn];
    //KM板子
    bool dfs(int x){
        visx[x]=true;
        for(int y=0;y<ny;y++){
            if(visy[y]) continue;
            int tmp=lx[x]+ly[y]-g[x][y];
            if(tmp==0){
                visy[y]=true;
                if(linker[y]==-1||dfs(linker[y])){
                    linker[y]=x;
                    return true;
                }
            }
            else if(slack[y]>tmp)
                slack[y]=tmp;
        }
        return false;
    }
    ll KM(){
        memset(linker,-1,sizeof(linker));
        memset(ly,0,sizeof(ly));
        for(int i=0;i<nx;i++){
            lx[i]=-INF;
            for(int j=0;j<ny;j++)
                if(g[i][j]>lx[i])
                    lx[i]=g[i][j];
        }
        for(int x=0;x<nx;x++){
            for(int i=0;i<ny;i++)
                slack[i]=INF;
            while(true){
                memset(visx,false,sizeof(visx));
                memset(visy,false,sizeof(visy));
                if(dfs(x)) break;
                int d=INF;
                for(int i=0;i<ny;i++){
                    if(!visy[i]&&d>slack[i])
                        d=slack[i];
                }
                for(int i=0;i<nx;i++)
                    if(visx[i])
                        lx[i]-=d;
                for(int i=0;i<ny;i++){
                    if(visy[i]) ly[i]+=d;
                    else slack[i]-=d;
                }
            }
        }
        ll res=0;
        for(int i=0;i<ny;i++){
            if(linker[i]!=-1)
                res+=g[linker[i]][i];
        }
        return res;
    }
    int main()
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++){
            scanf("%d%d%d",&a[i],&b[i],&k[i]);
        }
        for(int i=0;i<n;i++){
            for(int j=0;j<n;j++){//对于每个贷款与月份建边
                g[i][j]=max(a[i]-1ll*b[i]*(min(k[i],j)),0ll);
            }
        }
        ll res=0;
        nx=ny=n;
        res=KM();
        cout<<res<<endl;
    }
    
  • 相关阅读:
    Linux中怎么通过PID号找到对应的进程名及所在目录
    MYSQL 1093 之You can't specify target table for update in FROM clause解决办法
    Spring注解@Resource和@Autowired区别对比
    Java数据类型和MySql数据类型对应一览
    java高分局之jstat命令使用(转)
    为python安装matplotlib模块
    Python中的文件IO操作(读写文件、追加文件)
    Python 3语法小记(九) 异常 Exception
    SpringBoot下的Job定时任务
    linux的top命令参数详解
  • 原文地址:https://www.cnblogs.com/Chen-Jr/p/11007163.html
Copyright © 2011-2022 走看看