[ACM_几何] UVA 11300 Spreading the Wealth [分金币 左右给 最终相等 方程组 中位数]
Problem
A Communist regime is trying to redistribute wealth in a village. They have have decided to sit everyone around a circular table. First, everyone has converted all of their properties to coins of equal value, such that the total number of coins is divisible by the number of people in the village. Finally, each person gives a number of coins to the person on his right and a number coins to the person on his left, such that in the end, everyone has the same number of coins. Given the number of coins of each person, compute the minimum number of coins that must be transferred using this method so that everyone has the same number of coins.
The Input
There is a number of inputs. Each input begins with n(n<1000001), the number of people in the village. n lines follow, giving the number of coins of each person in the village, in counterclockwise order around the table. The total number of coins will fit inside an unsigned 64 bit integer.
The Output
For each input, output the minimum number of coins that must be transferred on a single line.
Sample Input
3 100 100 100 4 1 2 5 4
Sample Output
0 4
Problem setter: Josh Bao
题目大意:n人围绕圆桌,每个人有一定数量金币,金币总数能被n整除。每个人可以给他左右的人一些金币,最终使每个人金币数相等。求被转手金币最小值。
解题思路:
1、列出方程组:
设Ai为初始金币数,xi表示i号给i-1号金币数[可正可负],M为最终每个人的金币数,则:
for the first person: A1-x1+x2=M-->x2=M-A1+x1=x1-C1[规定C1=A1-M]
for the second person: A2-x2+x3=M-->x3=M-A2+x2=2M-A1-A2+x1=x1-C2
for the third person: A3-x3+x4=M-->x4=M-A3+x3=3M-A1-A2-A3+x1=x1-C3
......
for the nth person: An-xn+x1=M[这是一个多余的式子,并不能给我们更多的信息]
2、转换为单变量问题:
我们希望所有xi的绝对值之和最小,即:|x1|+|x1-C1|+|x1-C2|+...+|x1-Cn|最小。
3、转换为几何问题:
注意上式的几何意义就是数轴上给定n个点找出一个到他们之和尽量小的点。
4、递归求最优解:
这个最优的x1就是他们的中位数。
1 #include<iostream> 2 #include<algorithm> 3 #include<cstdio> 4 typedef long long LL; 5 using namespace std; 6 LL A[1000010],C[1000010];//A[i]表示第i人初始金币数 7 int main(){ 8 int n; 9 while(cin>>n){ 10 LL tot=0;//总数 11 for(int i=0;i<n;i++){ 12 cin>>A[i]; 13 tot+=A[i]; 14 } 15 LL M=tot/n;//平均数 16 C[0]=0; 17 for(int i=1;i<n;i++){ 18 C[i]=C[i-1]+A[i]-M;//递推数组C 19 } 20 sort(C,C+n); 21 LL x1=C[n/2],ans=0;//计算x1 22 for(int i=0;i<n;i++){ 23 ans+=abs(x1-C[i]); 24 } 25 cout<<ans<<'\n'; 26 }return 0; 27 }