NEERC 17 G.The Great Wall

最新推荐文章于 2021-05-22 21:16:35 发布

All_fade_away

最新推荐文章于 2021-05-22 21:16:35 发布

阅读量292

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.csdn.net/xmr_pursue_dreams/article/details/111381713

博客围绕NEERC 17 G.The Great Wall题目展开，先通过削减一维简化问题，确定a、b、c集合关系。接着令bi′、ci′，二分答案k，分x+r−1≥y和x+r−1<y两种情况用二维数点解决，最终得出时间复杂度为O(nlg2n)。

NEERC 17 G.The Great Wall

Solution

这题的第一步 $t r i c k$ 是：我们注意到取 $a, b, c$ 的集合两两不交且并集为 $U$ ，因此确定了 $b, c$ 之后可以简单地唯一确定 $a$ ，于是我们通过让 $bi−=ai,ci−=ai,Ans+=∑aib_i-=a_i,c_i-=a_i,Ans+=\sum a_i$ 来削减一维。

然后令 $bi′=∑j=1ibib'_i=\sum_{j=1}^ib_i$ ， $ci′=∑j=1icic'_i=\sum_{j=1}^ic_i$ ，二分答案 $k$ 。

当 $x+r−1≥yx+r-1\geq y$ ， $sum=cx+r−1′−cy−1′+by−1′−bx−1′+by+r−1′−bx+r−1′=(cx+r−1′−bx−1′−bx+r−1′)+(−cy−1′+by−1′+by+r−1′)≤ksum=c'_{x+r-1}-c'_{y-1}+b'_{y-1}-b'_{x-1}+b'_{y+r-1}-b'_{x+r-1}=(c'_{x+r-1}-b'_{x-1}-b'_{x+r-1})+(-c'_{y-1}+b'_{y-1}+b'_{y+r-1})\leq k$ ，可以简单地用二维数点解决。
当 $x + r - 1 < y$ ， $sum=(bx+r−1′−bx−1′)+(by+r−1′−by−1′)≤ksum=(b'_{x+r-1}-b'_{x-1})+(b'_{y+r-1}-b'_{y-1})\leq k$ ，同样可以简单地用二维数点解决。

时间复杂度 $O(n lg^2n)$ 。

Code

#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <cstring>
#include <ctime>
#include <cassert>
#include <string.h>
//#include <unordered_set>
//#include <unordered_map>
//#include <bits/stdc++.h>

#define MP(A,B) make_pair(A,B)
#define PB(A) push_back(A)
#define SIZE(A) ((int)A.size())
#define LEN(A) ((int)A.length())
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define fi first
#define se second

using namespace std;

template<typename T>inline bool upmin(T &x,T y) { return y<x?x=y,1:0; }
template<typename T>inline bool upmax(T &x,T y) { return x<y?x=y,1:0; }

typedef long long ll;
typedef unsigned long long ull;
typedef long double lod;
typedef pair<int,int> PR;
typedef vector<int> VI;

const lod eps=1e-11;
const lod pi=acos(-1);
const int oo=1<<30;
const ll loo=1ll<<62;
const int mods=998244353;
const int MAXN=600005;
const int INF=0x3f3f3f3f;//1061109567
/*--------------------------------------------------------------------*/
inline ll read()
{
	ll f=1,x=0; char c=getchar();
	while (c<'0'||c>'9') { if (c=='-') f=-1; c=getchar(); }
	while (c>='0'&&c<='9') { x=(x<<3)+(x<<1)+(c^48); c=getchar(); }
	return x*f;
}
int n,r,s[MAXN];
ll p[MAXN],b[MAXN],c[MAXN],B[MAXN],a[MAXN],k;
void add(int x,int y) { if (!x) return; for (;x<=n-r+1;x+=x&(-x)) s[x]+=y; }
int query(int x) { upmin(x,n-r+1); int ans=0; for (;x;x-=x&(-x)) ans+=s[x]; return ans; }
ll solve1(ll x)
{
	ll ans=0;
	for (int i=1;i<=n-r+1;i++) B[i]=p[i]=-c[i-1]+b[i-1]+b[i+r-1];
	sort(p+1,p+n-r+2);
	for (int i=1;i<=n-r+1;i++) B[i]=lower_bound(p+1,p+n-r+2,B[i])-p;
	for (int i=1;i<=r-1;i++) add(B[i],1);
	for (int i=1;i<=n-r+1;i++)
	{
		add(B[i],-1);
		if (i<=n-r-r+2) add(B[i+r-1],1);
		ans+=query(upper_bound(p+1,p+n-r+2,x-(c[i+r-1]-b[i-1]-b[i+r-1]))-p-1);
	}
	return ans;
}
ll solve2(ll x)
{
	ll ans=0;
	for (int i=1;i<=n-r+1;i++) B[i]=p[i]=b[i+r-1]-b[i-1];
	sort(p+1,p+n-r+2);
	for (int i=1;i<=n-r+1;i++) B[i]=lower_bound(p+1,p+n-r+2,B[i])-p;
	for (int i=r;i<=n-r+1;i++) add(B[i],1);
	for (int i=1;i<=n-r+1;i++)
	{
		if (i<=n-r-r+2) add(B[i+r-1],-1);
		ans+=query(upper_bound(p+1,p+n-r+2,x-(b[i+r-1]-b[i-1]))-p-1);
	}
	return ans;
}
signed main()
{
	ll sum=0;
	n=read(),r=read(),k=read();
	for (int i=1;i<=n;i++) a[i]=read(),sum+=a[i];
	for (int i=1;i<=n;i++) b[i]=read()-a[i]+b[i-1];
	for (int i=1;i<=n;i++) c[i]=read()-a[i]+c[i-1];
	
	ll l=0,r=1000000ll*n;
	while (l<r)
	{
		ll mid=(l+r)>>1,t=solve1(mid)+solve2(mid);
		if (t>=k) r=mid;
		else l=mid+1;
	}
	printf("%lld\n",r+sum);
	return 0;
}