1634: [Usaco2007 Jan]Protecting the Flowers 护花

1634: [Usaco2007 Jan]Protecting the Flowers 护花

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 493  Solved: 310
[Submit][Status]

Description

Farmer John went to cut some wood and left N (2 <= N <= 100,000) cows eating the grass, as usual. When he returned, he found to his horror that the cows were in his garden eating his beautiful flowers. Wanting to minimize the subsequent damage, FJ decided to take immediate action and transport the cows back to their barn. Each cow i is at a location that is Ti minutes (1 <= Ti <= 2,000,000) away from the barn. Furthermore, while waiting for transport, she destroys Di (1 <= Di <= 100) flowers per minute. No matter how hard he tries,FJ can only transport one cow at a time back to the barn. Moving cow i to the barn requires 2*Ti minutes (Ti to get there and Ti to return). Write a program to determine the order in which FJ should pick up the cows so that the total number of flowers destroyed is minimized.

   约翰留下他的N只奶牛上山采木．他离开的时候，她们像往常一样悠闲地在草场里吃草．可是，当他回来的时候，他看到了一幕惨剧：牛们正躲在他的花园里，啃食着他心爱的美丽花朵！为了使接下来花朵的损失最小，约翰赶紧采取行动，把牛们送回牛棚． 牛们从1到N编号．第i只牛所在的位置距离牛棚Ti(1≤Ti《2000000)分钟的路程，而在约翰开始送她回牛棚之前，她每分钟会啃食Di(1≤Di≤100)朵鲜花．无论多么努力，约翰一次只能送一只牛回棚．而运送第第i只牛事实上需要2Ti分钟，因为来回都需要时间．    写一个程序来决定约翰运送奶牛的顺序，使最终被吞食的花朵数量最小．

Input

* Line 1: A single integer

N * Lines 2..N+1: Each line contains two space-separated integers, Ti and Di, that describe a single cow's characteristics

 

    第1行输入N，之后N行每行输入两个整数Ti和Di．

 

Output

* Line 1: A single integer that is the minimum number of destroyed flowers

    一个整数，表示最小数量的花朵被吞食．

Sample Input

6
3 1
2 5
2 3
3 2
4 1
1 6

Sample Output

86

HINT

 

    约翰用6，2，3，4，1，5的顺序来运送他的奶牛．

 

Source

Silver

 题解：呵呵呵，又忘开int64了，结果WA了一次（长点记性吧QAQ）。。。我们可以这样分析此题——如题，假如按照T_i,D_i,T_i+1,D_i+1的顺序为最优解的话，则此时消耗的花为2T_iD_i+1,假如将这两个的顺序反过来的话，则消耗的花为2T_i+1D_i，由于前者为最优解，所以T_iD_i+1<T_i+1D_i，即(T_i/D_i)<(T_i+1/D_i+1)，则可以用冒泡排序说明顺序中各个牛的(T_i/D_i)必然递增（虽是这么分析的，但还是不建议先做商再排序，小心卡精度）。。。That's all...

 

var
   i,k,m,n:longint;
   j,l:int64;
   a,b:array[0..200000] of int64;
procedure swap(var x,y:int64);
          var z:int64;
          begin
               z:=x;x:=y;y:=z;
          end;
procedure sort(l,r:longint);
          var
             i,j:longint;
             x,y:int64;
          begin
               i:=l;j:=r;
               x:=a[(l+r) div 2];
               y:=b[(l+r) div 2];
               repeat
                     while (a[i]*y)<(x*b[i]) do inc(i);
                     while (a[j]*y)>(x*b[j]) do dec(j);
                     if i<=j then
                        begin
                             swap(a[i],a[j]);
                             swap(b[i],b[j]);
                             inc(i);dec(j);
                        end;
               until i>j;
               if l<j then sort(l,j);
               if i<r then sort(i,r);
          end;
begin
     readln(n);
     for i:=1 to n do
         begin
              readln(a[i],b[i]);
              a[i]:=a[i]*2;
         end;
     sort(1,n);
     j:=0;l:=0;
     for i:=1 to n-1 do
         begin
              j:=j+a[i];
              l:=l+j*b[i+1];
         end;
     writeln(l);
end.