题目描述
一条东西走向的穆西河将巴邻旁市一分为二，分割成了区域 $A$ 和区域 $B$。
每一块区域沿着河岸都建了恰好 $1000000001$ 栋的建筑，每条岸边的建筑都从 $0$ 编号到 $1000000000$。相邻的每对建筑相隔 $1$ 个单位距离，河的宽度也是 $1$ 个单位长度。区域 $A$ 中的 $i$ 号建筑物恰好与区域 $B$ 中的 $i$ 号建筑物隔河相对。
城市中有 $N$ 个居民。第 $i$ 个居民的房子在区域 $P_i$ 的 $S_i$ 号建筑上，同时他的办公室坐落在 $Q_i$ 区域的 $T_i$ 号建筑上。一个居民的房子和办公室可能分布在河的两岸，这样他就必须要搭乘船只才能从家中去往办公室，这种情况让很多人都觉得不方便。为了使居民们可以开车去工作，政府决定建造不超过 $K$ 座横跨河流的大桥。
由于技术上的原因，每一座桥必须刚好连接河的两岸，桥梁必须严格垂直于河流，并且桥与桥之间不能相交。
当政府建造最多 $K$ 座桥之后，设 $D_i$ 表示第 $i$ 个居民此时开车从家里到办公室的最短距离。请帮助政府建造桥梁，使得 $D_1 + D_2 + \dots + D_N$ 最小。
输入格式
输入的第一行包含两个正整数 $K$ 和 $N$，分别表示桥的上限数量和居民的数量。
接下来 $N$ 行，每一行包含四个参数：$P_i, S_i, Q_i$ 和 $T_i$，表示第 $i$ 个居民的房子在区域 $P_i$ 的 $S_i$ 号建筑上，且他的办公室位于 $Q_i$ 区域的 $T_i$ 号建筑上。
输出格式
输出仅为一行，包含一个整数，表示 $D_1 + D_2 + \dots + D_N$ 的最小值。
子任务
所有数据都保证：$P_i$ 和 $Q_i$ 为字符 “<samp>A</samp>” 和 “<samp>B</samp>” 中的一个， $0 \leq S_i, T_i \leq 1000000000$，同一栋建筑内可能有超过 $1$ 间房子或办公室（或二者的组合，即房子或办公室的数量同时大于等于 $1$）。
- 子任务 1 （8 分）
- $K = 1$
- $1 \leq N \leq 1000$
- 子任务 2 （14 分）
- $K = 1$
- $1 \leq N \leq 100000$
- 子任务 3 （9 分）
- $K = 2$
- $1 \leq N \leq 100$
- 子任务 4 （32 分）
- $K = 2$
- $1 \leq N \leq 1000$
- 子任务 5 （37 分）
- $K = 2$
- $1 \leq N \leq 100000$
分析1
我们先考虑\(m=1\)的时候的做法
对于子任务1，首先把所有的坐标离散出来，最多有\(2n\)个
然后暴力枚举\(2n\)的坐标的位置造桥，再\(O(n)\)计算代价
时间复杂度\(O(2n^2)\)，期望得分8分
紧接着考虑优化，我们发现，不需要每次花费\(O(n)\)来计算代价
我们只需要把所有的区间离散到坐标上，然后顺着扫一遍，直接记录一下，每次就可以\(O(1)\)修改得到代价
时间复杂度\(O(2n)\)，期望得分22分
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<string>
#include<climits>
#include<vector>
#include<cmath>
#include<map>
#include<queue>
#define LL long long
#define pii pair<int,int>
#define mp make_pair
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
return *p1++;
}
inline void read(int &x){
char c=nc();int b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
inline void read(LL &x){
char c=nc();LL b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
inline int read(char *s)
{
char c=nc();int len=1;
for(;!(c>='a' && c<='z');c=nc()) if (c==EOF) return 0;
for(;(c>='a' && c<='z');s[len++]=c,c=nc());
s[len++]='\0';
return len;
}
inline void read(char &x){
for (x=nc();!(x>='A' || x<='B');x=nc());
}
int wt,ss[19];
inline void print(int x){
if (x<0) x=-x,putchar('-');
if (!x) putchar(48); else {
for (wt=0;x;ss[++wt]=x%10,x/=10);
for (;wt;putchar(ss[wt]+48),wt--);}
}
inline void print(LL x){
if (x<0) x=-x,putchar('-');
if (!x) putchar(48); else {for (wt=0;x;ss[++wt]=x%10,x/=10);for (;wt;putchar(ss[wt]+48),wt--);}
}
int n,m,k;
LL b[2000010];
LL c[2000010],d[2000010];
struct data
{
int xa,xb;
LL ya,yb;
}a[1000010];
int Hash(int x)
{
return lower_bound(b+1,b+1+k,x)-b;
}
int main()
{
read(m);read(n);
char ch;k=0;
for (int i=1;i<=n;i++)
{
read(ch),read(a[i].ya),a[i].xa=ch-'A'+1,b[++k]=a[i].ya,
read(ch),read(a[i].yb),a[i].xb=ch-'A'+1,b[++k]=a[i].yb;
if (a[i].ya>a[i].yb) swap(a[i].ya,a[i].yb),swap(a[i].xa,a[i].xb);
}
sort(b+1,b+1+k);
k=unique(b+1,b+1+k)-b-1;
if (n<=1000 && m==1)
{
LL res=(LL)(1e17),s;
for (int i=1;i<=k;i++)
{
s=0;
for (int j=1;j<=n;j++)
if (a[j].xa==a[j].xb) s+=a[j].yb-a[j].ya;
else if (a[j].ya<=b[i] && a[j].yb>=b[i]) s+=a[j].yb-a[j].ya+1LL;
else if (b[i]<a[j].ya) s+=a[j].yb-a[j].ya+1+(a[j].ya-b[i])*2LL;
else s+=a[j].yb-a[j].ya+1LL+(b[i]-a[j].yb)*2LL;
res=min(res,s);
}
print(res),puts("");
return 0;
}
if (m==1 && n<=100000)
{
LL s=0,p=0,q=0;
memset(c,0,sizeof(c));
memset(d,0,sizeof(d));
for (int i=1;i<=n;i++)
if (a[i].xa==a[i].xb) s+=a[i].yb-a[i].ya;
else p++,c[Hash(a[i].ya)]++,d[Hash(a[i].yb)]++,s+=a[i].yb-a[i].ya+1+(a[i].ya-b[1])*2;
LL res=s;p-=c[1],q+=d[1];
for (int i=2;i<=k;i++)
{
s-=2*p*(b[i]-b[i-1]);p-=c[i];
s+=2*q*(b[i]-b[i-1]);q+=d[i];
res=min(res,s);
}
print(res),puts("");
return 0;
}
return 0;
}分析2
然后在思考\(m=2\)做法的时候，猛然醒悟...自己原来是个***...
我们在回过头去看\(m=1\)的时候
我们发现，答案由两部分组成，当两个建筑位于同一侧的时候\(ans1=|S_i-T_i|\)，显然\(ans1\)是固定的
而第二部分的答案\(ans2=\sum |P_i-x|+|T_i-x|+1\)，\(x\)表示桥的位置，显然当\(x\)集合\({S_i,T_i}\)的中位数的时候\(ans2\)取得最小值
时间复杂度\(O(nlogn)\)，期望得分22
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<string>
#include<climits>
#include<vector>
#include<cmath>
#include<map>
#include<queue>
#define LL long long
#define pii pair<int,int>
#define mp make_pair
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
return *p1++;
}
inline void read(int &x){
char c=nc();int b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
inline void read(LL &x){
char c=nc();LL b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
inline int read(char *s)
{
char c=nc();int len=1;
for(;!(c>='a' && c<='z');c=nc()) if (c==EOF) return 0;
for(;(c>='a' && c<='z');s[len++]=c,c=nc());
s[len++]='\0';
return len;
}
inline void read(char &x){
for (x=nc();!(x>='A' || x<='B');x=nc());
}
int wt,ss[19];
inline void print(int x){
if (x<0) x=-x,putchar('-');
if (!x) putchar(48); else {
for (wt=0;x;ss[++wt]=x%10,x/=10);
for (;wt;putchar(ss[wt]+48),wt--);}
}
inline void print(LL x){
if (x<0) x=-x,putchar('-');
if (!x) putchar(48); else {for (wt=0;x;ss[++wt]=x%10,x/=10);for (;wt;putchar(ss[wt]+48),wt--);}
}
int n,m,k;
LL b[2000010];
LL c[2000010],d[2000010];
struct data
{
int xa,xb;
LL ya,yb;
}a[1000010];
int main()
{
read(m);read(n);
char ch;k=0;
for (int i=1;i<=n;i++)
{
read(ch),read(a[i].ya),a[i].xa=ch-'A'+1;
read(ch),read(a[i].yb),a[i].xb=ch-'A'+1;
if (a[i].xa!=a[i].xb) b[++k]=a[i].ya,b[++k]=a[i].yb;
if (a[i].ya>a[i].yb) swap(a[i].ya,a[i].yb),swap(a[i].xa,a[i].xb);
}
sort(b+1,b+1+k);
if (m==1)
{
LL x=b[k/2+1],res=0;
for (int i=1;i<=n;i++)
if (a[i].xa==a[i].xb) res+=a[i].yb-a[i].ya;
else if (a[i].ya<=x && a[i].yb>=x) res+=a[i].yb-a[i].ya+1;
else if (a[i].ya>x) res+=a[i].yb-a[i].ya+1+2*(a[i].ya-x);
else res+=a[i].yb-a[i].ya+1+2*(x-a[i].yb);
print(res),puts("");
}
return 0;
}分析3
考虑\(m=2\)的做法
对于子任务3，同上述一样离散以后，暴力枚举两座桥的位置，然后花费\(O(n)\)计算代价
时间复杂度\(O(4n^3)\)，期望得分9分
继续考虑，我们可以发现，当有两座桥的时候，所有的居民唯一需要做的就是，找一座代价小的桥
而对于每一个居民而言，代价\(ans=|P_i-x|+|T_i-x|\)，\(x\)为桥的位置，即如下图
那么，增加相同的距离，也就是说，我们需要比较\(\frac {P_i+T_i}{2}\)距离哪个比较近即可
这样，我们如果把所有的居民按照\(\frac {P_i+T_i}{2}\)排序的话，我们可以把所有的居民分成连续的两段，前一段使用前一座桥，后一段时候后一座桥
那么，我们把\(m=2\)问题，规约成了两个\(m=1\)问题
枚举断点，暴力的做\(m=1\)问题即可
时间复杂度\(O(n^2logn)\)，期望得分41分，合并第一部分，可以得到63分，是一个很可观的分数了
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<string>
#include<climits>
#include<vector>
#include<cmath>
#include<map>
#include<queue>
#define LL long long
#define pii pair<int,int>
#define mp make_pair
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
return *p1++;
}
inline void read(int &x){
char c=nc();int b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
inline void read(LL &x){
char c=nc();LL b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
inline int read(char *s)
{
char c=nc();int len=1;
for(;!(c>='a' && c<='z');c=nc()) if (c==EOF) return 0;
for(;(c>='a' && c<='z');s[len++]=c,c=nc());
s[len++]='\0';
return len;
}
inline void read(char &x){
for (x=nc();!(x>='A' || x<='B');x=nc());
}
int wt,ss[19];
inline void print(int x){
if (x<0) x=-x,putchar('-');
if (!x) putchar(48); else {
for (wt=0;x;ss[++wt]=x%10,x/=10);
for (;wt;putchar(ss[wt]+48),wt--);}
}
inline void print(LL x){
if (x<0) x=-x,putchar('-');
if (!x) putchar(48); else {for (wt=0;x;ss[++wt]=x%10,x/=10);for (;wt;putchar(ss[wt]+48),wt--);}
}
int n,m,k;
LL b[2000010];
vector<int> e;
LL c[2000010],d[2000010];
struct data
{
int xa,xb;
LL ya,yb;
}a[1000010];
bool cmp(data x,data y)
{
return x.ya+x.yb<y.ya+y.yb;
}
int main()
{
read(m);read(n);
char ch;k=0;LL S=0;int x=0;
for (int i=1;i<=n;i++)
{
x++;
read(ch),read(a[x].ya),a[x].xa=ch-'A'+1;
read(ch),read(a[x].yb),a[x].xb=ch-'A'+1;
if (a[x].xa==a[x].xb) {S+=abs(a[x].ya-a[x].yb),x--;continue;}
b[++k]=a[x].ya,b[++k]=a[x].yb;
if (a[x].ya>a[x].yb) swap(a[x].ya,a[x].yb),swap(a[x].xa,a[x].xb);
}
n=x;
sort(b+1,b+1+k);
if (m==2 && n<=1000)
{
sort(a+1,a+1+n,cmp);
LL res=(LL)1e17,s,y;
for (int x=1;x<=n+1;x++)
{
s=0;
if (x>1)
{
e.clear();
for (int i=1;i<x;i++)
e.push_back(a[i].ya),e.push_back(a[i].yb);
sort(e.begin(),e.end());
y=e[e.size()/2];
for (int i=1;i<x;i++)
s+=abs(a[i].ya-y)+abs(a[i].yb-y)+1;
}
if (x<=n)
{
e.clear();
for (int i=x;i<=n;i++)
e.push_back(a[i].ya),e.push_back(a[i].yb);
sort(e.begin(),e.end());
y=e[e.size()/2];
for (int i=x;i<=n;i++)
s+=abs(a[i].ya-y)+abs(a[i].yb-y)+1;
}
res=min(res,s);
}
print(res+S),puts("");
}
return 0;
}分析4
继续考虑\(m=2\)的满分做法
我们得想办法把找中位数和算\(\sum |P_i-x|+|T_i-x|+1\)加速
我们每次会加入两个数字，或者删除两个数字，需要动态的维护中位数和上面的绝对值之和
用线段树来维护，很容易维护出中位数
那么绝对值怎么处理？显然，把绝对值拆开，因为，维护中位数的时候，顺便记录下中位数之前的和，然后就可以做啦
时间复杂度\(O(nlogn)\)，期望得分100
满分程序
我的线段树写的比较丑QAQ
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<string>
#include<climits>
#include<vector>
#include<cmath>
#include<map>
#include<queue>
#define LL long long
#define pii pair<int,int>
#define mp make_pair
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
return *p1++;
}
inline void read(int &x){
char c=nc();int b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
inline void read(LL &x){
char c=nc();LL b=1;
for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
inline int read(char *s)
{
char c=nc();int len=1;
for(;!(c>='a' && c<='z');c=nc()) if (c==EOF) return 0;
for(;(c>='a' && c<='z');s[len++]=c,c=nc());
s[len++]='\0';
return len;
}
inline void read(char &x){
for (x=nc();!(x>='A' || x<='B');x=nc());
}
int wt,ss[19];
inline void print(int x){
if (x<0) x=-x,putchar('-');
if (!x) putchar(48); else {
for (wt=0;x;ss[++wt]=x%10,x/=10);
for (;wt;putchar(ss[wt]+48),wt--);}
}
inline void print(LL x){
if (x<0) x=-x,putchar('-');
if (!x) putchar(48); else {for (wt=0;x;ss[++wt]=x%10,x/=10);for (;wt;putchar(ss[wt]+48),wt--);}
}
int n,m,k;
LL b[200010];
LL c[200010],d[200010];
struct data
{
int xa,xb;
LL ya,yb;
}a[100010];
struct ST
{
int sum;LL x;
}s1[800010],s2[800010];
int Hash(int x)
{
return lower_bound(b+1,b+1+k,x)-b;
}
void change1(int q,int l,int r,int x,int y,LL z)
{
if (l==r) {s1[x].sum+=y;s1[x].x+=z;return ;}
int mid=l+r>>1;
if (q<=mid) change1(q,l,mid,x<<1,y,z);else change1(q,mid+1,r,x<<1|1,y,z);
s1[x].sum=s1[x<<1].sum+s1[x<<1|1].sum;
s1[x].x=s1[x<<1].x+s1[x<<1|1].x;
}
void change2(int q,int l,int r,int x,int y,LL z)
{
if (l==r) {s2[x].sum+=y;s2[x].x+=z;return ;}
int mid=l+r>>1;
if (q<=mid) change2(q,l,mid,x<<1,y,z);else change2(q,mid+1,r,x<<1|1,y,z);
s2[x].sum=s2[x<<1].sum+s2[x<<1|1].sum;
s2[x].x=s2[x<<1].x+s2[x<<1|1].x;
}
int query1(int q,int l,int r,int x)
{
if (l==r) return b[l];
int mid=l+r>>1;
if (s1[x<<1].sum>=q) return query1(q,l,mid,x<<1);
return query1(q-s1[x<<1].sum,mid+1,r,x<<1|1);
}
int query2(int q,int l,int r,int x)
{
if (l==r) return b[l];
int mid=l+r>>1;
if (s2[x<<1].sum>=q) return query2(q,l,mid,x<<1);
return query2(q-s2[x<<1].sum,mid+1,r,x<<1|1);
}
LL sum1(int q,int l,int r,int x)
{
if (l==r) return s1[x].x;
int mid=l+r>>1;
if (s1[x<<1].sum>=q) return sum1(q,l,mid,x<<1);
return s1[x<<1].x+sum1(q-s1[x<<1].sum,mid+1,r,x<<1|1);
}
LL sum2(int q,int l,int r,int x)
{
if (l==r) return s2[x].x;
int mid=l+r>>1;
if (s2[x<<1].sum>=q) return sum2(q,l,mid,x<<1);
return s2[x<<1].x+sum2(q-s2[x<<1].sum,mid+1,r,x<<1|1);
}
LL num1(int q,int l,int r,int x)
{
if (l==r) return s1[x].sum;
int mid=l+r>>1;
if (s1[x<<1].sum>=q) return num1(q,l,mid,x<<1);
return s1[x<<1].sum+num1(q-s1[x<<1].sum,mid+1,r,x<<1|1);
}
LL num2(int q,int l,int r,int x)
{
if (l==r) return s2[x].sum;
int mid=l+r>>1;
if (s2[x<<1].sum>=q) return num2(q,l,mid,x<<1);
return s2[x<<1].sum+num2(q-s2[x<<1].sum,mid+1,r,x<<1|1);
}
bool cmp(data x,data y)
{
return x.ya+x.yb<y.ya+y.yb;
}
int main()
{
read(m);read(n);
char ch;k=0;LL S=0;int x=0;
for (int i=1;i<=n;i++)
{
x++;
read(ch),read(a[x].ya),a[x].xa=ch-'A'+1;
read(ch),read(a[x].yb),a[x].xb=ch-'A'+1;
if (a[x].xa==a[x].xb) {S+=abs(a[x].ya-a[x].yb),x--;continue;}
b[++k]=a[x].ya,b[++k]=a[x].yb;
if (a[x].ya>a[x].yb) swap(a[x].ya,a[x].yb),swap(a[x].xa,a[x].xb);
}
n=x;
sort(b+1,b+1+k);
if (m==1)
{
LL x=b[k/2+1],res=0;
for (int i=1;i<=n;i++)
if (a[i].ya<=x && a[i].yb>=x) res+=a[i].yb-a[i].ya+1;
else if (a[i].ya>x) res+=a[i].yb-a[i].ya+1+2*(a[i].ya-x);
else res+=a[i].yb-a[i].ya+1+2*(x-a[i].yb);
print(res+S),puts("");
}
if (m==2)
{
memset(s1,0,sizeof(s1));
memset(s2,0,sizeof(s2));
sort(a+1,a+1+n,cmp);
LL res=(LL)1e17,s,y,z,t,p=0,q=0;
for (int i=1;i<=n;i++)
change2(Hash(a[i].ya),1,k,1,1,a[i].ya),change2(Hash(a[i].yb),1,k,1,1,a[i].yb),q+=2LL;
for (int x=1;x<=n+1;x++)
{
s=0;
if (x>1)
{
y=query1(p/2+1,1,k,1);
z=sum1(p/2+1,1,k,1);
t=num1(p/2+1,1,k,1);
s+=y*t-z+(s1[1].x-z)-y*(p-t)+(LL)(x-1);
}
if (x<=n)
{
y=query2(q/2+1,1,k,1);
z=sum2(q/2+1,1,k,1);
t=num2(q/2+1,1,k,1);
s+=y*t-z+(s2[1].x-z)-y*(q-t)+(LL)(n-x+1);
}
if (x<=n)
{
change1(Hash(a[x].ya),1,k,1,1,a[x].ya),change1(Hash(a[x].yb),1,k,1,1,a[x].yb);p+=2LL;
change2(Hash(a[x].ya),1,k,1,-1,-a[x].ya),change2(Hash(a[x].yb),1,k,1,-1,-a[x].yb);q-=2LL;
}
res=min(res,s);
}
print(res+S),puts("");
}
return 0;
}
京公网安备 11010502036488号