题目链接: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);
}