Codeforces gym 102511 J

可以发现只有$1\sim 10^9$之间的$l$是有用的。
一个最暴力的想法是对每一个有用的$l$，分别算出每名选手在此时的总得分，然后可以算出每名选手此时的排名并更新对应的答案，可惜$l$的范围太大，这么做不太现实。不过我们可以换一个角度思考，我们只需要知道每名选手的最好排名，而排名的定义是分数不超过他的选手人数，那么我们可以尝试对于任意两名选手$i$和$j$，算出在哪些时刻$j$的得分不超过$i$的得分，这样也容易得到答案。事实上，下面我们将会证明，所有$j$得分不超过$i$得分的时刻，可以写为$\mathcal{O}(h)$个不相交区间的并。
对于第$i$名选手，显然他的得分关于$l$的函数$f(i,l)$是一个分成$\mathcal{O}(h)$段的一次函数，这个分段一次函数的表达式在对$a$数组排序后容易得出。那么对于任意两名选手$i$和$j$，令$g(i,j,l)=f(i,l)-f(j,l)$，则当$l=k$时$j$的得分不超过$i$的得分当且仅当$g(i,j,k)\le 0$。由于$f(i,l)$和$f(j,l)$都是$\mathcal{O}(h)$段的分段一次函数，因此$g(i,j,l)$也是$\mathcal{O}(h)$段的分段一次函数。我们对于$g(i,j,l)$的每一段分别考虑，显然每一段$\le 0$的时刻都形成了一段区间，因此总体来看，$g(i,j,l)\le 0$的$l$可写为$\mathcal{O}(h)$个不相交区间的并。
这也给出了本题的算法：对于每名选手$i$，我们先计算出$f(i,l)$的分段一次函数表达式。接着对于第$i$名选手，枚举$j=1\sim n$，计算出$g(i,j,l)$的分段一次函数表达式，然后就容易将$j$得分不超过$i$得分的时刻写为$\mathcal{O}(h)$个不相交区间的并。对于第$i$名选手，每有一个区间$[L,R]$，意味着当$l\in [L,R]$时第$i$名选手的排名要增加$1$。这样我们可以将一个区间$[L,R]$拆成两个事件$(L,+1)$和$(R+1,-1)$，对于这些事件按发生时间排序后扫描一遍，维护每个时刻第$i$名选手的排名，并更新答案即可。
时间复杂度瓶颈主要在最后的排序部分，总时间复杂度为$\mathcal{O}(p^{2}h\log (ph))$。

#include <bits/stdc++.h>
#define FR first
#define SE second

using namespace std;

typedef long long ll;
typedef pair<int,int> pr;

int val[100005],sz;

int num[505][55];
int pos[105];

void calc(int x,int y,int m) {
  int sz1=0,r1=1,r2=1;
  while (r1<=m||r2<=m) {
  	if (r1<=m&&(r2>m||num[x][r1]<=num[y][r2])) {
  		pos[++sz1]=num[x][r1];
  		r1++;
	  }
  	else {
	    pos[++sz1]=num[y][r2];
	    r2++;
	}
  	while (r1<=m&&num[x][r1]<=pos[sz1]) r1++;
  	while (r2<=m&&num[y][r2]<=pos[sz1]) r2++;
  } 
  pos[++sz1]=1e9+1;
  int v1=m,v2=m;
  r1=r2=1;
  ll s1=0,s2=0;
  for(int i=1;i<=sz1;i++) {
  	while (r1<=m&&num[x][r1]<=pos[i-1]) {
  		r1++;
  		v1--;
	  } 
	while (r2<=m&&num[y][r2]<=pos[i-1]) {
		r2++;
		v2--;
	}
	if (s1<s2) {
		if (s1+(ll)v1*(pos[i]-pos[i-1])>=s2+(ll)v2*(pos[i]-pos[i-1])) {
			int t=((s2-s1)+(ll)(v1-v2)*pos[i-1]-1)/(v1-v2);
			val[++sz]=((t+1)<<1)+1;
		}
	}
	else {
		if (s1+(ll)v1*(pos[i]-pos[i-1])<s2+(ll)v2*(pos[i]-pos[i-1])) {
			int t=((s1-s2)+(ll)(v2-v1)*pos[i-1])/(v2-v1);
			val[++sz]=((t+1)<<1);
		}
	}
	s1+=(ll)v1*(pos[i]-pos[i-1]);
	s2+=(ll)v2*(pos[i]-pos[i-1]);
  }
}

int solve(int x,int n,int m) {
  sz=0;
  for(int i=1;i<=n;i++)
    if (x!=i) calc(x,i,m);
  sort(val+1,val+sz+1);
  int ans=n,s=n;
  for(int i=1,j=1;i<=sz;i=j) {
  	while (j<=sz&&(val[j]>>1)==(val[i]>>1)) j++;
  	for(int k=i;k<j;k++)
  	  if (val[k]&1) s++; else s--;
	ans=min(ans,s); 
  }
  return ans;
}

int main() {
  int n,m;
  scanf("%d%d",&n,&m);
  for(int i=1;i<=n;i++) {
  	for(int j=1;j<=m;j++) scanf("%d",&num[i][j]);
  	sort(num[i]+1,num[i]+m+1);
  }
  for(int i=1;i<=n;i++) {
  	int ans=solve(i,n,m);
  	printf("%d\n",ans);
  }
  return 0;
}