KM模板

var
 n,m,i,j:longint;
 ans:int64;
 sel,lx,ly,slack:array[..] of int64;
 a:array[..,..] of int64;
 visx,visy:array[..] of boolean;

 function dfs(po:longint):boolean;
 var
  i:longint;
  t:int64;
   begin
     visx[po]:=true;

     for i:= to m do
       if visy[i]=false then
         begin
           t:=lx[po]+ly[i]-a[po,i];
           if t= then
             begin
               visy[i]:=true;
               if sel[i]= then
                begin
                  sel[i]:=po;
                  exit(true);
                end else
               if dfs(sel[i]) then
                 begin
                   sel[i]:=po;
                   exit(true);
                 end;
             end else
             if slack[i]>t then slack[i]:=t;
         end;
     exit(false);
   end;

 procedure km;
 var
  i,j:longint;
  d:int64;
   begin
     for i:= to n do
       lx[i]:=-*maxlongint;

     for i:= to n do
       for j:= to m do
         if a[i,j]>lx[i] then lx[i]:=a[i,j];

     for i:= to n do
       begin
         while true do
           begin
             fillchar(visx,sizeof(visx),);
             fillchar(visy,sizeof(visy),);
             for j:= to m do
             slack[j]:=*maxlongint;

             if dfs(i) then break;

             d:=*maxlongint;

             for j:= to  m do
               if (visy[j]=false) and (slack[j]<d) then d:=slack[j];

             for j:= to n do
               if visx[j] then dec(lx[j],d);

             for j:= to m do
               if visy[j] then inc(ly[j],d);
           end;
       end;
   end;

  begin
    read(n,m);

    for i:= to n do
      for j:= to m do
       begin
         read(a[i,j]);
       end;

    km;

    for i:= to n do ans:=ans+lx[i];

    for i:= to  m do ans:=ans+ly[i];

  writeln(ans);
  end.
#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;

  int m,n,k;
  int visx[],visy[],sel[];
  long long lx[],ly[],a[][],slack[];

  int dfs(int po){
      visx[po]=;
      for (int i=;i<=m;i++)
          if (visy[i]==){
              long long  t=lx[po]+ly[i]-a[po][i];
              if  (t==){
                  visy[i]=;
                  if ((sel[i]==-)||(dfs(sel[i]))){
                      sel[i]=po;
                      return();
                  }
              }else slack[i]=min(slack[i],t);
          }
      return();
  }

  void km(){
      for (int i=;i<=n;i++) lx[i]=-1e10;
      for (int i=;i<=m;i++) ly[i]=;
      for (int i=;i<=m;i++) sel[i]=-;
      for (int i=;i<=n;i++)
          for(int j=;j<=m;j++)
               lx[i]=max(lx[i],a[i][j]);    

      for (int i=;i<=n;i++)
           while (){
               for (int j=;j<=n;j++) visx[j]=;
               for (int j=;j<=m;j++) visy[j]=,slack[j]=1e10;

               if (dfs(i)) break;    

               long long d=1e10;
               for (int j=;j<=m;j++)
                   if (visy[j]==)
                       d=min(d,slack[j]);

               for (int j=;j<=n;j++)
                    if (visx[j])  lx[j]-=d;

               for (int j=;j<=m;j++)
                   if (visy[j])  ly[j]+=d;
           }
  }

  int main(){
      int T;
      scanf("%d",&T);
      while (T--){
          scanf("%d%d%d",&n,&m,&k);

          for (int i=;i<=n;i++)
              for (int j=;j<=m;j++)
                   a[i][j]=;

          for (int i=;i<=k;i++){
              int x,l,p,r;
              scanf("%d%d%d%d",&x,&l,&r,&p);
              for (int j=l;j<=r;j++) a[x][j]=p;
          }

          km();
          long long  ans1=;
          for (int i=;i<=n;i++) ans1+=lx[i];
          for (int i=;i<=m;i++) ans1+=ly[i];

          for (int i=;i<=n;i++)
               for (int j=;j<=m;j++)
                    a[i][j]*=-;

          km();
          long long ans2=;
          for (int i=;i<=n;i++) ans2-=lx[i];
          for (int i=;i<=m;i++) ans2-=ly[i];

          if (ans1==ans2) printf("YES\n");else printf("NO\n");
      }
  }

求最小值所有边取相反数后KM，输出有答案的相反数

求最小乘积先求对数再KM