jzoj5236. 【NOIP2017模拟8.7A组】利普希茨

Description

Input

输入文件名为lipschitz.in。
第一行一个整数n。
接下来一行n个整数，描述序列A。
第三行一个数q 。
接下来q行，每行三个整数。其中第一个整数type表示操作的类型。 type=0对应修改操作， type=1对应查询操作。

Output

输出文件名为lipschitz.out。
对于每个查询，给出f(A[l…r]) 。

Sample Input

输入1：
6
90 50 78 0 96 20
6
0 1 35
1 1 4
0 1 67
0 4 11
0 3 96
1 3 5

输入2：
50
544 944 200 704 400 150 8 964 666 596 850 608 452 103 988 760 370 723 350 862 856 0 724 544 668 891 575 448 16 613 952 745 990 459 740 960 752 194 335 575 525 12 618 80 618 224 240 600 562 283
10
1 6 6
1 1 3
0 11 78279
0 33 42738
0 45 67270
1 1 26
1 19 24
1 37 39
1 8 13
0 7 64428

Sample Output

输出1：
78
85

输出2：
0
744
77683
856
558
77683

Data Constraint

对于30%的数据，n,q<=500
对于60%的数据，n,q<=5000
对于100%的数据，n,q<=100000,0<=ai,val<=10^9

题外话

大爷省选加油啊。
我这个菜逼就在图书馆当个加油的吧。

题解

啊，一道包着大数据结构的sb结论题。
本着想到一个神奇的错误结论能过掉的，然后就没有然后了。
这结论也是真的妙，自己还是太菜了，没想到。

我们可以考虑那条神奇的式子：
$$\lceil \frac{|a_i-a_j|}{|i-j|}\rceil$$
那么我们考虑把上面这个丑陋的式子换个方式看：
那么我们考虑把每个点都变成一个点：$(i,a_i)$
然后上面的式子不就变成了求斜率吗？

然后画个图就可以发现，我们对于每三个点，其中必定有一个相邻的点的斜率最小。

线段树即可。

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
const int maxn=100010;

int n,m,a[maxn],ans,xx,yy,ki,x,y,jia,answer,wz;
struct node{
	int ma,mi,wz1,wz2;
}tree[maxn*4];

void findmax(int x,int l,int r,int st,int en)
{
	if (l==st && r==en)
	{
		if (tree[x].ma>answer) wz=tree[x].wz1;
		answer=max(answer,tree[x].ma);
	}
	else
	{
		int mid=(l+r)/2;
		if (en<=mid) findmax(x*2,l,mid,st,en);
		else if (st>mid) findmax(x*2+1,mid+1,r,st,en);
		else
		{
			findmax(x*2,l,mid,st,mid);
			findmax(x*2+1,mid+1,r,mid+1,en);
		}
	}
}

void findmin(int x,int l,int r,int st,int en)
{
	if (l==st && r==en)
	{
		if (tree[x].mi<answer) wz=tree[x].wz2;
		answer=min(answer,tree[x].mi);
	}
	else
	{
		int mid=(l+r)/2;
		if (en<=mid) findmin(x*2,l,mid,st,en);
		else if (st>mid) findmin(x*2+1,mid+1,r,st,en);
		else
		{
			findmin(x*2,l,mid,st,mid);
			findmin(x*2+1,mid+1,r,mid+1,en);
		}
	}
}

void modify(int x,int l,int r,int st,int en)
{
	if (l==r)
	{
		tree[x].ma=en;tree[x].mi=en;
		tree[x].wz1=l;tree[x].wz2=l;
	}
	else
	{
		int mid=(l+r)/2;
		if (st<=mid) modify(x*2,l,mid,st,en);
		else modify(x*2+1,mid+1,r,st,en);
		if (tree[x*2].ma>tree[x*2+1].ma) tree[x].wz1=tree[x*2].wz1;else tree[x].wz1=tree[x*2+1].wz1;
		if (tree[x*2].mi<tree[x*2+1].mi) tree[x].wz2=tree[x*2].wz2;else tree[x].wz2=tree[x*2+1].wz2;
		tree[x].ma=max(tree[x*2].ma,tree[x*2+1].ma);
		tree[x].mi=min(tree[x*2].mi,tree[x*2+1].mi);
	}
}

__attribute__((optimize("-O3")))
int main()
{
	freopen("lipschitz.in","r",stdin);
	freopen("lipschitz.out","w",stdout);
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
	{
		scanf("%d",&a[i]);
		if (i>1)
		modify(1,1,n-1,i-1,abs(a[i]-a[i-1]));
	}
	scanf("%d",&m);
	while (m>0)
	{
		m--;
		scanf("%d%d%d",&ki,&x,&y);
		if (ki==0)
		{
			a[x]=y; 
			modify(1,1,n-1,x-1,abs(a[x]-a[x-1]));
			modify(1,1,n-1,x,abs(a[x]-a[x+1]));
		}
		else
		{
			if (x==y) printf("0\n");
			else
			{
				answer=0;
				wz=0;
				findmax(1,1,n-1,x,y-1);
				printf("%d\n",answer);
			}
		}
	}
}