码迷,mamicode.com
首页 > 其他好文 > 详细

POJ之01背包系列

时间:2016-08-17 21:11:59      阅读:130      评论:0      收藏:0      [点我收藏+]

标签:

poj3624 Charm Bracelet

模板题

没有要求填满,所以初始化为0就行

#include<cstdio>
#include<iostream>
using namespace std;
#define N 15010
int n,m,v[N],c[N],f[N];
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d%d",&v[i],&c[i]);
    for(int i=1;i<=n;i++){
        for(int j=m;j>=v[i];j--){
            f[j]=max(f[j],f[j-v[i]]+c[i]);
        }
    }
    printf("%d\n",f[m]);
    return 0;
}

 

poj3628 Bookshelf 2

 

#include<cstdio>
#include<iostream>
using namespace std;
#define N 21000000
int n,m,sum,ans,c[25],f[N];
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&c[i]),sum+=c[i];
    for(int i=1;i<=n;i++){
        for(int j=sum;j>=c[i];j--){
            f[j]=max(f[j],f[j-c[i]]+c[i]);
        }
    }
    for(int i=1;i<=sum;i++){
        if(f[i]>=m){
            ans=f[i]-m;break;
        }
    }
    printf("%d\n",ans);
    return 0;
} 

 

 

 

 

poj1745 Divisibility

 

这道题如果取摸后范围比较小,第二维不大于100,然后记忆化背包(非递归搜索),根据能达到的值推能达到的值,衍生出前n个时的所有状态。

 

#include<cstdio>
#include<iostream>
using namespace std;
#define N 10005
int n,k,c[N];
bool f[N][105];
int main(){
    scanf("%d%d",&n,&k);
    for(int i=0,x;i<n;i++){
        scanf("%d",&x);
        if(x<0) x=-x;
        c[i]=x%k;
    }
    f[0][c[0]]=1;
    for(int i=1;i<n;i++){
        for(int j=0;j<k;j++){
            if(f[i-1][j]){
                f[i][(j+c[i])%k]=1;
                f[i][(j+k-c[i])%k]=1;
            }
        }
    }
    puts(f[n-1][0]?"Divisible":"Not divisible");
    return 0;
} 

 

 

 

 

poj1976  A Mini Locomotive

 

3辆车运货,共有n堆货,每辆可以运连续k堆,求最大运货量

 

保证k*3<=n;也就是说要运货量最大必须堆数为k。

 

f[i][j]为前j次前i堆最大运货量的最大运货量

 

如果i<j*k;那么只能全运

 

第i堆不运:f[i-1][j]

 

第i堆运:f[i-k][j-1]+sum;

 

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define N 50010
int n,k,T,s[N],f[N][4];
int main(){
    scanf("%d",&T);
    while(T--){
        memset(f,0,sizeof f);
        scanf("%d",&n);
        for(int i=1;i<=n;i++) scanf("%d",&s[i]),s[i]+=s[i-1];
        scanf("%d",&k);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=3;j++){
                if(i<j*k)
                    f[i][j]=i;
                else
                    f[i][j]=max(f[i-1][j],f[i-k][j-1]+s[i]-s[i-k]);
            }
        }
        printf("%d\n",f[n][3]);
    }
    return 0;
}

 

 

 

 

poj1837 Balance

 

状态压缩求方案数

注意要平移数组

 

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define p 7500
#define N 21
int n,m,h[N],w[N],f[N][p*2+10];
int main(){
    scanf("%d%d",&m,&n);
    for(int i=0;i<m;i++) scanf("%d",&h[i]);
    for(int i=1;i<=n;i++) scanf("%d",&w[i]);
    f[0][p]=1;
    for(int i=1;i<=n;i++){
        for(int k=0;k<=p<<1;k++){
            if(f[i-1][k]){
                for(int j=0;j<m;j++){
                    f[i][k+h[j]*w[i]]+=f[i-1][k];
                }
            }
        }
    }
    printf("%d\n",f[n][p]);
    return 0;
}

 

 

 

poj1948  Triangular Pastures

 

n个线,组一个三角形,求三角形的最大面积。

 

海伦公式

 

  • 技术分享,p为半周长

 

因为周长已知,知道两条边就能确定面积。

 

设f[i][j](i>j),然后确定每一条边是否加入那个边,由已知状态推出未知状态,随之更新最大面积。

 

每条边不可能超过周长的一半。

 

#include<cstdio>
#include<cmath>
#include<iostream>
using namespace std;
#define M 45
#define N 1605
int n,c,a[M];
bool f[N][N];
int ans=-0x7fffffff;
int check(int x,int y){
    double t1=x;
    double t2=y;
    double t3=c-x-y;
    if(t1+t2<t3||t1+t3<t2||t3+t2<t1) return -1.0;
    double t=c*1.0/2.0;
    double res=sqrt(t*(t-t1)*(t-t2)*(t-t3))*100.0;
    return res;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",a+i),c+=a[i];
    f[0][0]=1;
    for(int i=1;i<=n;i++){
        for(int j=c>>1;j>=0;j--){
            for(int k=c>>1;k>=j;k--){
                if(j>=a[i]&&f[j-a[i]][k]) f[j][k]=1; 
                if(k>=a[i]&&f[j][k-a[i]]) f[j][k]=1; 
                if(f[j][k]) ans=max(ans,check(j,k)); 
            }
        }
    }
    printf("%d\n",ans);
    return 0;
}

 

POJ之01背包系列

标签:

原文地址:http://www.cnblogs.com/shenben/p/5781739.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!