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