Big Event in HDU（HDU1171）可用背包和母函数求解

Big Event in HDU  HDU1171

就是求一个简单的背包：

题意：就是给出一系列数，求把他们尽可能分成均匀的两堆

如：2 10 1 20 1     结果是：20 10。才最均匀！

三种解法：

多重背包的优化与否：（1031MS）

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;
int dp[250005];
int a[55],b[55];
int main()
{
    int n,s,i,j,k;
    while(scanf("%d",&n)!=EOF)
    {
        if(n<0)
            break;
        s=0;
        for(i=0;i<n;i++)
        {
            scanf("%d%d",&a[i],&b[i]);
            s+=a[i]*b[i];
        }
        memset(dp,0,sizeof(dp));
        for(i=0;i<n;i++)
            for(j=1;j<=b[i];j++)
            for(k=s/2;k>=a[i];k--)
            dp[k]=max(dp[k],dp[k-a[i]]+a[i]);
        if(dp[s/2]>s-dp[s/2])
            printf("%d %d
",dp[s/2],s-dp[s/2]);
        else
            printf("%d %d
",s-dp[s/2],dp[s/2]);
    }
    return 0;
}

二进制优化的：（46MS）

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<string.h>
using namespace std;
int dp[250005],a1[250005];
int main()
{
    int a[5005],b[5005],n,i,j,k,s,cout1;
    while(scanf("%d",&n)!=EOF)
    {
        if(n<0)
            break;
        s=0;
        cout1=0;
        for(i=0;i<n;i++)
        {
            scanf("%d%d",&a[i],&b[i]);
            s+=a[i]*b[i];
            for(k=1;k<=b[i];k<<=1)
            {
                a1[cout1++]=k*a[i];
                b[i]-=k;
            }
            if(b[i]>0)
                a1[cout1++]=b[i]*a[i];
        }
        memset(dp,0,sizeof(dp));
        for(i=0;i<cout1;i++)
            for(j=s/2;j>=a1[i];j--)
            dp[j]=max(dp[j],dp[j-a1[i]]+a1[i]);
        if(dp[s/2]>=s-dp[s/2])
            printf("%d %d
",dp[s/2],s-dp[s/2]);
        else
            printf("%d %d
",s-dp[s/2],dp[s/2]);
    }
    return 0;
}

最后是母函数求解的：（875MS）

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;
int c1[250005],c2[250005];
int main()
{
    int n,i,j,k,s;
    int a[55],b[55];
    while(scanf("%d",&n)!=EOF)
    {
        s=0;
        if(n<0)
            break;
        for(i=1;i<=n;i++)
        {
            scanf("%d%d",&a[i],&b[i]);
            s+=a[i]*b[i];
        }
        for(i=0;i<=s;i++)
        {
            c1[i]=0;
            c2[i]=0;
        }
        c1[0]=1;
        for(i=1;i<=n;i++)
        {
            for(j=0;j<=s;j++)
                for(k=0;k+j<=s,k<=a[i]*b[i];k+=a[i])
                c2[k+j]+=c1[j];
            for(j=0;j<=s;j++)
            {
                c1[j]=c2[j];
                c2[j]=0;
            }
        }
        for(i=s/2;i<=s;i++)
        {
            if(c1[i]!=0)
                break;
        }
        if(i>=s-i)//且要取最大值
        printf("%d %d
",i,s-i);
        else
        printf("%d %d
",s-i,i);
    }
    return 0;
}