poj 3111 迭代法

给N个珠宝，每个珠宝有v与w两个参量，从中取K个使得∑v/∑w 最大。

利用迭代法，首先任取K个最为初始值，计算s，然后对每个珠宝计算其y=v-s*w的值，根据y值对其排序，然后取前K个作为本次选取的K个，直到选取的K个珠宝计算出的s值不再增加为止。明显只需要证明s是单调递增的即可。

对于当前K个珠宝计算值s有∑v-s*∑w=0，即∑y=0，但是根据y排序后选取的前K个必有∑y>=0, 即∑v/∑w>=s,即s是单调递增的。

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<iostream>
#include<math.h>
using namespace std;
const double esp=1e-8;
int N,K;
struct Je
{
    int v,w,id;
    double y;
}je[100005];
int cmp(Je je1,Je je2)
{
    return je1.y>je2.y;
}
int main()
{
//    freopen("test.txt","r",stdin);
    cin>>N>>K;
    for(int i=0;i<N;i++)
    {
        scanf("%d%d",&je[i].v,&je[i].w);
        je[i].y=0;je[i].id=i+1;
    }
    double ps=1,s=0;
    while(fabs(s-ps)>esp)
    {
        for(int i=0;i<N;i++)
            je[i].y=je[i].v-s*je[i].w;
        sort(je,je+N,cmp);
        double sv=0,sw=0;
        for(int i=0;i<K;i++)
        {
            sv+=je[i].v;sw+=je[i].w;
        }
        ps=s;s=sv/sw;
    }
    for(int i=0;i<K;i++)
        printf("%d ",je[i].id);
    printf("\n");
    return 0;
}