hdu 4412 Sky Soldiers （dp）

具体思路的解释见  http://hi.baidu.com/renxl51/item/d7437a30bca43883f4e4ad0d

f[i][j] 表示用j个点覆盖前i个点。  f[i][j]=min{ f[k][j-1] + cost }  ( j-1 <= k <= i-1 )   就是枚举第j个点覆盖了从第k+1个点到第i个点  ，  其中cost的计算可以利用单调性

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<map>
#define INF ((1<<23)-1)
using namespace std;
map<int ,double> mp;
double f[1010][60];
typedef struct
{
    int x; double p;
}Point;
Point pt[1010];
int kk,m;

void initial()
{
    for(int i=0;i<=1000;i++)
        for(int j=0;j<=50;j++)
            f[i][j]=INF;
    for(int i=0;i<=m;i++)
        for(int j=i;j<=m;j++)
            f[i][j]=0.0;
}

int main()
{
    int k,i,j,l,x; double p;
    while( scanf("%d%d",&kk,&m),kk+m )
    {
        mp.clear();
        initial();
        for(i=1;i<=kk;i++)
        {
            scanf("%d",&l);
            for(j=1;j<=l;j++)
            {
                scanf("%d%lf",&x,&p);
                mp[x]+=p;
            }
        }
        map<int,double>::iterator it;
        j=1;
        for( it=mp.begin(); it!=mp.end();it++,j++)
        {
            pt[j].x=(it->first);
            pt[j].p=(it->second);
        }
        int n=j-1;
        for( i=1;i<=n;i++ )
        {
            for(j=1;j<=m && j<=i; j++ )
            {
                f[i][j]=INF;
                double pre=0.000,suf=0.000;
                int cur=i; double add=0.000;
                for(k=i-1;k>=j-1;k--)
                {
                    if( k+1 != i )
                    {
                        pre+=pt[k+1].p;
                        add+=pt[k+1].p*(pt[cur].x-pt[k+1].x);
                    }
                    while( pt[cur].p+suf < pre )
                    {
                        suf+=pt[cur].p;
                        cur--;
                        pre-=pt[cur].p;
                        add-=pt[cur].p*(pt[cur+1].x-pt[cur].x);
                        add+=pt[cur+1].p*(pt[cur+1].x-pt[cur].x);
                        add-=pre*(pt[cur+1].x-pt[cur].x);
                        add+=(suf-pt[cur+1].p)*(pt[cur+1].x-pt[cur].x);
                    }
                    f[i][j]=min(f[i][j],f[k][j-1]+add);
                }
            }
        }
        printf("%.2lf\n",f[n][m]);
    }
    return 0;
}