迭代学习

原题： 给n个项链，每个项链有价值vi，重量wi的属性，
       选出k个，让SUM(vi)/SUM(wi)最大。

 原题地址：http://acm.pku.edu.cn/JudgeOnline/problem?id=3111

采用方法是迭代：
令F（x）= (Vi1+Vi2+...+Vik) - x*(Wi1+Wi2+...+Wik);
F（x）关于x单调，故收敛。
原题在这个式子可以体现成F(x)=0,求x的最大值。 从x=0开始迭代，
按照V - x*W非升序排列，令F(x) =0, 计算x', 如果|x-x'|<一定值，则迭代

结束，否则继续迭代。

做了这个题：感觉好多数学方法人在第一思考中用不到，经常用的往往是最基本的想法，但那样可能就解决不了问题。看来要多多练习呀。
这道题我搞了好长时间，现在搞出来了，很爽，呵呵！

程序如下：   AC过了

#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;

struct Data
{
 int v,w,id;
};
#define MAX 100002
Data d[MAX]; double s;

bool gre(const Data & m1, const Data & m2)
{
 return (double)((double)m1.v - s*(double)m1.w) > (double)

((double)m2.v - s*(double)m2.w) ;
}

int main()
{
  //  freopen("C://ACMData.txt","r",stdin);
 int n,k,i; double st;
 _int64 V,W;
 while( scanf("%d%d",&n,&k)!=EOF)
 {
        i=0;
  while( i<n){
            scanf("%d%d", &d[i].v,&d[i].w);
   d[i].id = i+1;
   i++;
  }
  s =1.0;  st=0.0;

  while( fabs( s-st) >=1e-9){
   s = st;
   sort(d,d+n,gre);
   V=0;W =0;i=0;
      while(i<k){
      V+=d[i].v;
      W+=d[i].w;
      i++;
   }
   st = (double)V/(double)W;
  } 
        i=0;
  while(i<k)
  {
   printf("%d ",d[i].id); i++;
  }
  printf("/n");

 }
 return 0;
}