背包变形 思路

居民集会

蓝桥村的居民都生活在一条公路的边上，公路的长度为L，每户家庭的位置都用这户家庭到公路的起点的距离来计算，第i户家庭距起点的距离为di。
每年，蓝桥村都要举行一次集会。今年，由于村里的人口太多，村委会决定要在4个地方举行集会，其中3个位于公路中间，1个位最公路的终点。
已知每户家庭都会向着远离公路起点的方向去参加集会，参加集会的路程开销为家庭内的人数ti与距离的乘积。
给定每户家庭的位置di和人数ti，请为村委会寻找最好的集会举办地：p1, p2, p3, p4 (p1<=p2<=p3<=p4=L),使得村内所有人的路程开销和最小。
【输入格式】
输入的第一行包含两个整数n, L，分别表示蓝桥村的家庭数和公路长度。
接下来n行，每行两个整数di, ti，分别表示第i户家庭距离公路起点的距离和家庭中的人数。
【输出格式】
输出一行，包含一个整数，表示村内所有人路程的开销和。
【样例输入】
6 10
1 3
2 2
4 5
5 20
6 5
8 7
【样例输出】
18
【样例说明】
在距起点2, 5, 8, 10这4个地方集会，6个家庭需要的走的距离分别为1, 0, 1, 0, 2, 0，总的路程开销为1*3+0*2+1*5+0*20+2*5+0*7=18。
【数据规模与约定】
对于10%的评测数据，1<=n<=300。
对于30%的评测数据，1<=n<=2000，1<=L<=10000，0<=di<=L，di<=di+1，0<=ti<=20。
对于100%的评测数据，1<=n<=100000，1<=L<=1000000，0<=di<=L，di<=di+1，0<=ti<=1000000。
资源约定：
峰值内存消耗 < 512M
CPU消耗  < 5000ms

参考大佬的思路，是不是大佬都不写注释，想了好久才想明白。。。

思路：把居民看成是排列好的一堆物品，其价值就是，居民的人数到聚会地点的距离。然后开会的地点看成是从这个位置开始将左边的某一段物品全部放进背包，共计拿三次，求出背包的最小值。上式dp[i][j]代表从第i个位置向左第j次拿取物品。

代码：

#include<iostream>  
#include<string>  
#include<cstdio>  
#include<algorithm>  
#include<cmath>  
#include<cstring>  
#include<set>  
using namespace std;  
#define INF 1.0e14;  
const int maxn = 1e5+100;  
typedef long long LL;  
int A[maxn];  
int B[maxn];  
LL dp[maxn][4];  
LL sum[maxn][4];  
int main(){  
    int n,m;  
    scanf("%d%d",&n,&m);  
    int i,j,k;  
    for(i=1;i<=n;i++)  
        scanf("%d%d",&A[i],&B[i]);  
    for(i=1;i<=n;i++){  
        dp[i][0]=dp[i-1][0]+sum[i-1][0]*(A[i]-A[i-1]); 
        sum[i][0]=sum[i-1][0]+B[i];///从第一个点向左拿 
        
        for(j=1;j<=3;j++)  
         if(i==j){///如果是相等则不能再向左边取了。  
          dp[i][j]=dp[i-1][j-1];  
            sum[i][j]=B[i];  
         }  
        else if(i>j){  
            dp[i][j]=dp[i-1][j]+sum[i-1][j]*(A[i]-A[i-1]);  
            sum[i][j]=sum[i-1][j]+B[i];  
            if(dp[i][j]>=dp[i-1][j-1])  
               dp[i][j]=dp[i-1][j-1],sum[i][j]=B[i];  
        }///i<j 没东西拿，直接是零 
    }  
    LL ans=dp[n][0]; ///只有终点一个集会地点，剩余三个都在起点 
    if(n>1) ans=min(ans,dp[n][1]);///如果n==2，最多只需要2个聚会地点  
    if(n>2) ans=min(ans,dp[n][2]);  
    if(n>3) ans=min(ans,dp[n][3]+sum[n][3]*(m-A[n]));  
    cout<<ans<<endl;  
    return 0;  
}  

同时这题要注意的是因为可以在一个点拿多次，这种情况下相当于聚会地点小于三个，所以最后还要再比较一下输出答案。

网上有其他做法，但是感觉对于10%的数据应该会TLE，但是比赛时能写出来得分就很不错了。

#include<vector>  
#include<iostream>  
#include<cstring>  
#include<algorithm>  
using namespace std;  
const int maxn=10000+5;   
int d[maxn],t[maxn],n,L;  
int w[maxn][maxn];  
  
int solve(int x,int y,int cur)  //[x,y] 还能分cur次   
{  
    if(!cur) return w[x][y];  
      
    int ML,MR,M;  
    M=2*w[0][L];  
    for(int m=x+1;m<y;m++)  
    {  
        ML=solve(x,m,cur-1);  //分为两个区间   
        MR=solve(m,y,cur-1);  
        M=min(M,ML+MR);  
    }  
    return M;  
}  
int main()  
{  
    scanf("%d%d",&n,&L);  
    for(int i=0;i<n;i++)  
        scanf("%d%d",&d[i],&t[i]);  
    for(int i=0;i<=L;i++)  //[i,j] w[i][j]为区间i~j所需的开销   
    for(int j=i;j<=L;j++)  
    {  
        for(int k=0;k<n;k++)  
        {  
            if(d[k]>i&&d[k]<j)  
            {  
                w[i][j]+=(j-d[k])*t[k];  
            }  
        }  
    }  
    cout<<solve(0,L,2);  
    return 0;  
}