[cf1219G]Harvester

分类讨论（以下仅考虑行，列的情况）：
1.4行的，求出每一行的和后找到4个最大值即可；
2.3行1列，枚举列，再将每一行最大值减去那一列的值后取3个最大值得和即可；
3.2行2列，发现行和列是等价的，因此可以通过旋转使得$mle n$，又因为$nmle 10^5$，所以$mle sqrt{10^5}$，枚举两列后用同样的做法求出最大值即可
最终复杂度$o(nmcodt min(n,m))$，也就是根号的复杂度

#include<bits/stdc++.h>
using namespace std;
#define N 100005
#define ll long long
int n,m;
ll ans,b[N],s1[N],s2[N];
vector<int>a[N];
bool cmp(ll x,ll y){
    return x>y;
}
ll qu(ll *a,int n,int k){
    if (n<k)return 0;
    ll x;
    for(int i=1;i<=k;i++)b[i]=0;
    for(int i=1;i<=n;i++){
        x=a[i];
        for(int j=1;j<=k;j++)
            if (x>b[j])swap(x,b[j]);
    }
    x=0;
    for(int i=1;i<=k;i++)x+=b[i];
    return x;
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        a[i].push_back(0);
        for(int j=1;j<=m;j++){
            scanf("%lld",&ans);
            a[i].push_back(ans);
            s1[i]+=a[i][j];
            s2[j]+=a[i][j];
        }
    }
    ans=max(qu(s1,n,4),qu(s2,m,4));
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++)s2[j]-=a[i][j];
        ans=max(ans,s1[i]+qu(s2,m,3));
        for(int j=1;j<=m;j++)s2[j]+=a[i][j];
    }
    for(int j=1;j<=m;j++){
        for(int i=1;i<=n;i++)s1[i]-=a[i][j];
        ans=max(ans,s2[j]+qu(s1,n,3));
        for(int i=1;i<=n;i++)s1[i]+=a[i][j];
    }
    if (n<=m)
        for(int i=1;i<n;i++)
            for(int j=i+1;j<=n;j++){
                for(int k=1;k<=m;k++)s2[k]-=a[i][k]+a[j][k];
                ans=max(ans,s1[i]+s1[j]+qu(s2,m,2));
                for(int k=1;k<=m;k++)s2[k]+=a[i][k]+a[j][k];
            }
    else
        for(int i=1;i<m;i++)
            for(int j=i+1;j<=m;j++){
                for(int k=1;k<=n;k++)s1[k]-=a[k][i]+a[k][j];
                ans=max(ans,s2[i]+s2[j]+qu(s1,n,2));
                for(int k=1;k<=n;k++)s1[k]+=a[k][i]+a[k][j];
            }
    printf("%lld",ans);
}