Codeforces 394D Physical Education and Buns

题目链接：http://codeforces.com/problemset/problem/394/D

题目大意：

完成一些操作，使得给出的n(n<=1000)个数A1,A2,A3..An(doesn't exceed 104.)（可随意调换位置）构成一个不降的等差数列。

可以对任意一个数加、减1，令第i个数的加减次数为c[i]。

寻找一个等差数列使得每个数加减次数最大的最小(  Max{c[i]} (1<=i<=n) 最小 )，输出这个最小值和该等差数列的首项和公差。

无视上面，语死早就不说了。。

先肯定升序排个序。。

由于是让最大的最小，可以考虑二分答案，正好也有二分性质。lol~

二分最小值k，然后每个数的变动范围不超过k，数i的范围就可以确定 [ a[i] - k,a[i] + k ]

看成二维平面上竖着的一些区间，找一条y=kx+b k,b都为整数的直线穿过去，如果找得到就是有，找不到就没有。

斜率一定的时候，把该直线向上一直移动，直到碰到一个区间的上端点，  <- 这只是说明若找得到直线就一定有一条经过区间上端点，

例如 K = 1的情况下：

枚举区间上端点，维护一个斜率（视角）范围，到最后看范围内有没有整数 

这样就是 O(n^2 logn ) 

。。。。这是蒟蒻的方法大家可以无视掉。。。。。

其实可以枚举斜率，

令直线是 y=d*(x-1)+b 等差数列那个东西  其中b是未知数

直线到每一个点的有向距离就是 deltay[i]=(i-1)*d+b-a[i]

现在我们的目的就是使Max(deltay[i]) 最小，

deltay[i]=(i-1)*d-a[i]+b 调整b使得这个东西最小，

那么算出(i-1)*d-a[i]的最大值Max 和最小值 Min ,  MIN(Max(deltay))=(Max-Min+1)/2 （整除。加1只是为了四舍五入），就是数轴上有一个长为Max-Min的区间，选一个点到区间端点的最大值最小，当然就是中点

得到MINMAXdeltay了，首项？

Max(deltay)=b-min

deltay[i]=(i-1)*d-a[i] +b

当(i-1)*d-a[i]=Min 的时候deltay最大

所以b=Max(deltay)+min

y=d*(x-1)+b

首项就是y=b

公差是枚举的，所以就在O(n*公比d)的时间解决，每个数不超过10000，所以Max公比d=20000

不过对于枚举公比只能呵呵了。。。

感觉这个Minmax和公比(斜率)成一个凹函数的关系，是否可以考虑三分啊。（你去写吧）

#include <cstdio>
#include <algorithm>
#define rep(i,l,r) for (int i=l;i<=r;++i)
int getx(){
	char c;int x;bool pd=false;
	for (c=getchar();c!='-'&&(c<'0'||c>'9');c=getchar());
	if (c=='-') c=getchar(),pd=true;
	for (x=0;c>='0'&&c<='9';c=getchar())
		x=(x<<3)+(x<<1)+c-'0';
	return pd?-x:x;
}
bool upmin(int &a,const int &b){return a>b?a=b,1:0;}
bool upmax(int &a,const int &b){return a<b?a=b,1:0;}
const int MAX_N=1050;
const int INF=~0U>>2;
int a[MAX_N];
int n;
int minN=INF,minD,minF;
int main(){
	n=getx();
	rep(i,1,n) a[i]=getx();
	std::sort(a+1,a+n+1);
	rep(d,0,20000){
		int min=INF,max=-INF;
		rep(i,1,n){
			upmin(min,a[i]-(i-1)*d);
			upmax(max,a[i]-(i-1)*d);
			}
		int now=(max-min+1)/2;
		if (upmin(minN,now)){
			minD=d;
			minF=now+min;
			}
		}
	printf("%d\n%d %d\n",minN,minF,minD);
}