枚举 + 三分 (游标)
The swimming area of Berhattan's city beach is marked out with n buoys. The buoys form a straight line. When the buoys were being put into the water, nobody cared to observe the same distance between each pair of adjacent buoys.
Now the beach keeper wants the distance between any two adjacent buoys to be the same. He plans to shift some or all of the buoys without changing their respective order. To facilitate the task, he wants the total length of all shifts to be as small as possible.
Given coordinates of the buoys, you should find the minimum possible length of all shifts, as well as new coordinates of the buoys.
Input
The first line of input contains a single integer n (2 ≤ n ≤ 400), n — the number of buoys. The second line contains buoys' integer coordinates x1, x2, ..., xn ( - 10000 ≤ xi ≤ 10000). No two given buoys will share the same place. The coordinates are given in strictly increasing order.
To the first line print a real number t — the minimum possible total length of required shifts. Output this value with at least 4 digits after the decimal point.
To the second line print n numbers — new coordinates of the buoys. The new coordinates should be printed in strictly increasing order with at least 7 digits after the decimal point. If there are several optimal ways to shift the buoys, you may output any of them.
4
-2 2 6 9
1.0000
-2.0000000000 1.6666666667 5.3333333333 9.0000000000
All buoys are located on the Ox axis. You may move buoys only along the Ox axis.
题意 : 在 X 轴上有一串浮标,移动一些浮标,让他们之间的距离相等,问最小移动的总距离之和是多少,并输出此时的浮标位置
思路 : 对于两个浮标间的距离 X ,当X过大或过小的时候都会使答案很大,因此中间存在一个最小值,这里三分就可以,然后既然要使距离最小,那我们就先确保一个浮标不动去移动其他的,这里枚举就可以。
代码示例:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 | #define ll long long const int maxn = 1e6+5; const double pi = acos (-1.0); const int inf = 0x3f3f3f3f; double pre[405]; double ans[405]; double sum = 1.0*inf, lenn; int n, pos; double fun( double len, int base){ double x = pre[base]; double s = 0; for ( int i = base-1; i >= 1; i--){ x -= len; s += fabs (x-pre[i]); } x = pre[base]; for ( int i = base+1; i <= n; i++){ x += len; s += fabs (x-pre[i]); } if (s < sum) { sum = s; lenn = len; pos = base; } return s; } int main() { freopen ( "input.txt" , "r" ,stdin); freopen ( "output.txt" , "w" ,stdout); cin >> n; for ( int i = 1; i <= n; i++){ scanf ( "%lf" , &pre[i]); } for ( int i = 1; i <= n; i++){ double l = 0, r = 200000; for ( int j = 1; j <= 100; j++){ double lm = l + (r-l)/3; double rm = r - (r-l)/3; double lf = fun(lm, i); double rf = fun(rm, i); if (lf > rf) l = lm; else r = rm; } } ans[pos] = pre[pos]; for ( int i = pos+1; i <= n; i++) ans[i] = ans[i-1] + lenn; for ( int i = pos-1; i >= 1; i--) ans[i] = ans[i+1] - lenn; printf ( "%.4f\n" , sum); for ( int i = 1; i <= n; i++) printf ( "%.10f%c" , ans[i], i==n? '\n' : ' ' ); return 0; } |
作者:静默虚空
欢迎任何形式的转载,但请务必注明出处。
限于本人水平,如果文章和代码有表述不当之处,还请不吝赐教。
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步