分支定界法 0-1多背包问题 0-1背包问题的多种解法代码(动态规划、贪心法、回溯法、分支...

\u7528\u5206\u652f\u9650\u754c\u6cd5\u6c42\u89e30/1\u80cc\u5305\u95ee\u9898

1.\u95ee\u9898\u63cf\u8ff0\uff1a\u5df2\u77e5\u6709N\u4e2a\u7269\u54c1\u548c\u4e00\u4e2a\u53ef\u4ee5\u5bb9\u7eb3M\u91cd\u91cf\u7684\u80cc\u5305\uff0c\u6bcf\u79cd\u7269\u54c1I\u7684\u91cd\u91cf\u4e3aWEIGHT\uff0c\u4e00\u4e2a\u53ea\u80fd\u5168\u653e\u5165\u6216\u8005\u4e0d\u653e\u5165\uff0c\u6c42\u89e3\u5982\u4f55\u653e\u5165\u7269\u54c1\uff0c\u53ef\u4ee5\u4f7f\u80cc\u5305\u91cc\u7684\u7269\u54c1\u7684\u603b\u6548\u76ca\u6700\u5927\u3002

2.\u8bbe\u8ba1\u601d\u60f3\u4e0e\u5206\u6790\uff1a\u5bf9\u7269\u54c1\u7684\u9009\u53d6\u4e0e\u5426\u6784\u6210\u4e00\u68f5\u89e3\u6811\uff0c\u5de6\u5b50\u6811\u8868\u793a\u4e0d\u88c5\u5165\uff0c\u53f3\u8868\u793a\u88c5\u5165\uff0c\u901a\u8fc7\u68c0\u7d22\u95ee\u9898\u7684\u89e3\u6811\u5f97\u51fa\u6700\u4f18\u89e3\uff0c\u5e76\u7528\u7ed3\u70b9\u4e0a\u754c\u6740\u6b7b\u4e0d\u7b26\u5408\u8981\u6c42\u7684\u7ed3\u70b9\u3002

#include

struct good
{
int weight;
int benefit;
int flag;//\u662f\u5426\u53ef\u4ee5\u88c5\u5165\u6807\u8bb0
};

int number=0;//\u7269\u54c1\u6570\u91cf
int upbound=0;
int curp=0, curw=0;//\u5f53\u524d\u6548\u76ca\u503c\u4e0e\u91cd\u91cf
int maxweight=0;
good *bag=NULL;

void Init_good()
{
bag=new good [number];

for(int i=0; i {
cout>bag[i].weight;
cout>bag[i].benefit;
bag[i].flag=0;//\u521d\u59cb\u6807\u5fd7\u4e3a\u4e0d\u88c5\u5165\u80cc\u5305
cout< }

}

int getbound(int num, int *bound_u)//\u8fd4\u56de\u672c\u7ed3\u70b9\u7684c\u9650\u754c\u548cu\u9650\u754c
{
for(int w=curw, p=curp; num {
w=w+bag[num].weight;
p=w+bag[num].benefit;
}

*bound_u=p+bag[num].benefit;
return ( p+bag[num].benefit*((maxweight-w)/bag[num].weight) );
}

void LCbag()
{
int bound_u=0, bound_c=0;//\u5f53\u524d\u7ed3\u70b9\u7684c\u9650\u754c\u548cu\u9650\u754c

for(int i=0; i {
if( ( bound_c=getbound(i+1, &bound_u) )>upbound )//\u904d\u5386\u5de6\u5b50\u6811
upbound=bound_u;//\u66f4\u6539\u5df2\u6709u\u9650\u754c,\u4e0d\u66f4\u6539\u6807\u5fd7

if( getbound(i, &bound_u)>bound_c )//\u904d\u5386\u53f3\u5b50\u6811
//\u82e5\u88c5\u5165\uff0c\u5224\u65ad\u53f3\u5b50\u6811\u7684c\u9650\u754c\u662f\u5426\u5927\u4e8e\u5de6\u5b50\u6811\u6839\u7684c\u9650\u754c\uff0c\u662f\u5219\u88c5\u5165
{
upbound=bound_u;//\u66f4\u6539\u5df2\u6709u\u9650\u754c
curp=curp+bag[i].benefit;
curw=curw+bag[i].weight;//\u4ece\u5df2\u6709\u91cd\u91cf\u548c\u6548\u76ca\u52a0\u4e0a\u65b0\u7269\u54c1
bag[i].flag=1;//\u6807\u8bb0\u4e3a\u88c5\u5165
}
}

}

void Display()
{

cout<<"\u53ef\u4ee5\u653e\u5165\u80cc\u5305\u7684\u7269\u54c1\u7684\u7f16\u53f7\u4e3a\uff1a";
for(int i=0; iif(bag[i].flag>0)

cout<
\uff5d


\u53c2\u8003\uff1a
http://cache.baidu.com/c?word=%B7%D6%D6%A7%3B%CF%DE%BD%E7%3B%B7%A8%3B%C7%F3%BD%E2%3B0%2C1%3B%B1%B3%B0%FC%3B%CE%CA%CC%E2&url=http%3A//www%2Edaxiongonline%2Ecom/dvbbs/dispbbs%2Easp%3FboardID%3D31%26ID%3D1016&b=10&a=0&user=baidu

\u4e00.\u52a8\u6001\u89c4\u5212\u6c42\u89e30-1\u80cc\u5305\u95ee\u9898
/************************************************************************/
/* 0-1\u80cc\u5305\u95ee\u9898\uff1a
/* \u7ed9\u5b9an\u79cd\u7269\u54c1\u548c\u4e00\u4e2a\u80cc\u5305
/* \u7269\u54c1i\u7684\u91cd\u91cf\u4e3awi\uff0c\u5176\u4ef7\u503c\u4e3avi
/* \u80cc\u5305\u7684\u5bb9\u91cf\u4e3ac
/* \u5e94\u5982\u4f55\u9009\u62e9\u88c5\u5165\u80cc\u5305\u7684\u7269\u54c1\uff0c\u4f7f\u5f97\u88c5\u5165\u80cc\u5305\u4e2d\u7684\u7269\u54c1
/* \u7684\u603b\u4ef7\u503c\u6700\u5927\uff1f
/* \u6ce8\uff1a\u5728\u9009\u62e9\u88c5\u5165\u80cc\u5305\u7684\u7269\u54c1\u65f6\uff0c\u5bf9\u7269\u54c1i\u53ea\u6709\u4e24\u79cd\u9009\u62e9\uff0c
/* \u5373\u88c5\u5165\u6216\u4e0d\u88c5\u5165\u80cc\u5305\u3002\u4e0d\u80fd\u5c06\u7269\u54c1i\u88c5\u5165\u591a\u6b21\uff0c\u4e5f
/* \u4e0d\u80fd\u53ea\u88c5\u5165\u90e8\u5206\u7684\u7269\u54c1i\u3002
/*
/* 1. 0-1\u80cc\u5305\u95ee\u9898\u7684\u5f62\u5f0f\u5316\u63cf\u8ff0\uff1a
/* \u7ed9\u5b9ac>0, wi>0, vi>0, 0<=i<=n\uff0c\u8981\u6c42\u627e\u5230\u4e00\u4e2an\u5143\u7684
/* 0-1\u5411\u91cf(x1, x2, ..., xn), \u4f7f\u5f97\uff1a
/* max sum_{i=1 to n} (vi*xi),\u4e14\u6ee1\u8db3\u5982\u4e0b\u7ea6\u675f\uff1a
/* (1) sum_{i=1 to n} (wi*xi) <= c
/* (2) xi\u2208{0, 1}, 1<=i<=n
/*
/* 2. 0-1\u80cc\u5305\u95ee\u9898\u7684\u6c42\u89e3
/* 0-1\u80cc\u5305\u95ee\u9898\u5177\u6709\u6700\u4f18\u5b50\u7ed3\u6784\u6027\u8d28\u548c\u5b50\u95ee\u9898\u91cd\u53e0\u6027\u8d28\uff0c\u9002\u4e8e
/* \u91c7\u7528\u52a8\u6001\u89c4\u5212\u65b9\u6cd5\u6c42\u89e3
/*
/* 2.1 \u6700\u4f18\u5b50\u7ed3\u6784\u6027\u8d28
/* \u8bbe(y1,y2,...,yn)\u662f\u7ed9\u5b9a0-1\u80cc\u5305\u95ee\u9898\u7684\u4e00\u4e2a\u6700\u4f18\u89e3\uff0c\u5219\u5fc5\u6709
/* \u7ed3\u8bba\uff0c(y2,y3,...,yn)\u662f\u5982\u4e0b\u5b50\u95ee\u9898\u7684\u4e00\u4e2a\u6700\u4f18\u89e3\uff1a
/* max sum_{i=2 to n} (vi*xi)
/* (1) sum_{i=2 to n} (wi*xi) <= c - w1*y1
/* (2) xi\u2208{0, 1}, 2<=i<=n
/* \u56e0\u4e3a\u5982\u82e5\u4e0d\u7136\uff0c\u5219\u8be5\u5b50\u95ee\u9898\u5b58\u5728\u4e00\u4e2a\u6700\u4f18\u89e3(z2,z3,...,zn)\uff0c
/* \u800c(y2,y3,...,yn)\u4e0d\u662f\u5176\u6700\u4f18\u89e3\u3002\u90a3\u4e48\u6709\uff1a
/* sum_{i=2 to n} (vi*zi) > sum_{i=2 to n} (vi*yi)
/* \u4e14\uff0cw1*y1 + sum_{i=2 to n} (wi*zi) <= c
/* \u8fdb\u4e00\u6b65\u6709\uff1a
/* v1*y1 + sum_{i=2 to n} (vi*zi) > sum_{i=1 to n} (vi*yi)
/* w1*y1 + sum_{i=2 to n} (wi*zi) <= c
/* \u8fd9\u8bf4\u660e\uff1a(y1,z2,z3,...zn)\u662f\u6240\u7ed90-1\u80cc\u5305\u95ee\u9898\u7684\u66f4\u4f18\u89e3\uff0c\u90a3\u4e48
/* \u8bf4\u660e(y1,y2,...,yn)\u4e0d\u662f\u95ee\u9898\u7684\u6700\u4f18\u89e3\uff0c\u4e0e\u524d\u63d0\u77db\u76fe\uff0c\u6240\u4ee5\u6700\u4f18
/* \u5b50\u7ed3\u6784\u6027\u8d28\u6210\u7acb\u3002
/*
/* 2.2 \u5b50\u95ee\u9898\u91cd\u53e0\u6027\u8d28
/* \u8bbe\u6240\u7ed90-1\u80cc\u5305\u95ee\u9898\u7684\u5b50\u95ee\u9898 P(i,j)\u4e3a\uff1a
/* max sum_{k=i to n} (vk*xk)
/* (1) sum_{k=i to n} (wk*xk) <= j
/* (2) xk\u2208{0, 1}, i<=k<=n
/* \u95ee\u9898P(i,j)\u662f\u80cc\u5305\u5bb9\u91cf\u4e3aj\u3001\u53ef\u9009\u7269\u54c1\u4e3ai,i+1,...,n\u65f6\u7684\u5b50\u95ee\u9898
/* \u8bbem(i,j)\u662f\u5b50\u95ee\u9898P(i,j)\u7684\u6700\u4f18\u503c\uff0c\u5373\u6700\u5927\u603b\u4ef7\u503c\u3002\u5219\u6839\u636e\u6700\u4f18
/* \u5b50\u7ed3\u6784\u6027\u8d28\uff0c\u53ef\u4ee5\u5efa\u7acbm(i,j)\u7684\u9012\u5f52\u5f0f\uff1a
/* a. \u9012\u5f52\u521d\u59cb m(n,j)
/* //\u80cc\u5305\u5bb9\u91cf\u4e3aj\u3001\u53ef\u9009\u7269\u54c1\u53ea\u6709n\uff0c\u82e5\u80cc\u5305\u5bb9\u91cfj\u5927\u4e8e\u7269\u54c1n\u7684
/* //\u91cd\u91cf\uff0c\u5219\u76f4\u63a5\u88c5\u5165\uff1b\u5426\u5219\u65e0\u6cd5\u88c5\u5165\u3002
/* m(n,j) = vn, j>=wn
/* m(n,j) = 0, 0<=j<wn
/* b. \u9012\u5f52\u5f0f m(i,j)
/* //\u80cc\u5305\u5bb9\u91cf\u4e3aj\u3001\u53ef\u9009\u7269\u54c1\u4e3ai,i+1,...,n
/* //\u5982\u679c\u80cc\u5305\u5bb9\u91cfj<wi\uff0c\u5219\u6839\u672c\u88c5\u4e0d\u8fdb\u7269\u54c1i\uff0c\u6240\u4ee5\u6709\uff1a
/* m(i,j) = m(i+1,j), 0<=j<wi
/* //\u5982\u679cj>=wi\uff0c\u5219\u5728\u4e0d\u88c5\u7269\u54c1i\u548c\u88c5\u5165\u7269\u54c1i\u4e4b\u95f4\u505a\u51fa\u9009\u62e9
/* \u4e0d\u88c5\u7269\u54c1i\u7684\u6700\u4f18\u503c\uff1am(i+1,j)
/* \u88c5\u5165\u7269\u54c1i\u7684\u6700\u4f18\u503c\uff1am(i+1, j-wi) + vi
/* \u6240\u4ee5\uff1a
/* m(i,j) = max {m(i+1,j), m(i+1, j-wi) + vi}, j>=wi
/*
/************************************************************************/



#define max(a,b) (((a) > (b)) ? (a) : (b))
#define min(a,b) (((a) < (b)) ? (a) : (b))
template
void Knapsack(Type* v, int *w, int c, int n, Type **m)
{
//\u9012\u5f52\u521d\u59cb\u6761\u4ef6
int jMax = min(w[n] - 1, c);
for (int j=0; j<=jMax; j++) {
m[n][j] = 0;
}

for (j=w[n]; j<=c; j++) {
m[n][j] = v[n];
}

//i\u4ece2\u5230n-1\uff0c\u5206\u522b\u5bf9j>=wi\u548c0<=j<wi\u5373\u4f7fm(i,j)
for (int i=n-1; i>1; i--) {
jMax = min(w[i] - 1, c);
for (int j=0; j<=jMax; j++) {
m[i][j] = m[i+1][j];
}
for (j=w[i]; j<=c; j++) {
m[i][j] = max(m[i+1][j], m[i+1][j-w[i]]+v[i]);
}
}

m[1][c] = m[2][c];
if (c >= w[1]) {
m[1][c] = max(m[1][c], m[2][c-w[1]]+v[1]);
}

}

template
void TraceBack(Type **m, int *w, int c, int n, int* x)
{
for (int i=1; i<n; i++) {
if(m[i][c] == m[i+1][c]) x[i] = 0;
else {
x[i] = 1;
c -= w[i];
}
}
x[n] = (m[n][c])? 1:0;
}


int main(int argc, char* argv[])
{
int n = 5;
int w[6] = {-1, 2, 2, 6, 5, 4};
int v[6] = {-1, 6, 3, 5, 4, 6};
int c = 10;

int **ppm = new int*[n+1];
for (int i=0; i<n+1; i++) {
ppm[i] = new int[c+1];
}

int x[6];

Knapsack(v, w, c, n, ppm);
TraceBack(ppm, w, c, n, x);

return 0;
}
\u4e8c.\u8d2a\u5fc3\u7b97\u6cd5\u6c42\u89e30-1\u80cc\u5305\u95ee\u9898
1.\u8d2a\u5fc3\u6cd5\u7684\u57fa\u672c\u601d\u8def\uff1a
\u2014\u2014\u4ece\u95ee\u9898\u7684\u67d0\u4e00\u4e2a\u521d\u59cb\u89e3\u51fa\u53d1\u9010\u6b65\u903c\u8fd1\u7ed9\u5b9a\u7684\u76ee\u6807\uff0c\u4ee5\u5c3d\u53ef\u80fd\u5feb\u7684\u5730\u6c42\u5f97\u66f4\u597d\u7684\u89e3\u3002\u5f53\u8fbe\u5230\u67d0\u7b97\u6cd5\u4e2d\u7684\u67d0\u4e00\u6b65\u4e0d\u80fd\u518d\u7ee7\u7eed\u524d\u8fdb\u65f6\uff0c\u7b97\u6cd5\u505c\u6b62\u3002
\u8be5\u7b97\u6cd5\u5b58\u5728\u95ee\u9898\uff1a
1).\u4e0d\u80fd\u4fdd\u8bc1\u6c42\u5f97\u7684\u6700\u540e\u89e3\u662f\u6700\u4f73\u7684\uff1b
2).\u4e0d\u80fd\u7528\u6765\u6c42\u6700\u5927\u6216\u6700\u5c0f\u89e3\u95ee\u9898\uff1b
3).\u53ea\u80fd\u6c42\u6ee1\u8db3\u67d0\u4e9b\u7ea6\u675f\u6761\u4ef6\u7684\u53ef\u884c\u89e3\u7684\u8303\u56f4\u3002

\u5b9e\u73b0\u8be5\u7b97\u6cd5\u7684\u8fc7\u7a0b\uff1a
\u4ece\u95ee\u9898\u7684\u67d0\u4e00\u521d\u59cb\u89e3\u51fa\u53d1\uff1b
while \u80fd\u671d\u7ed9\u5b9a\u603b\u76ee\u6807\u524d\u8fdb\u4e00\u6b65 do
\u3000\u3000 \u6c42\u51fa\u53ef\u884c\u89e3\u7684\u4e00\u4e2a\u89e3\u5143\u7d20\uff1b
\u7531\u6240\u6709\u89e3\u5143\u7d20\u7ec4\u5408\u6210\u95ee\u9898\u7684\u4e00\u4e2a\u53ef\u884c\u89e3\uff1b

2.\u4f8b\u9898\u5206\u6790

1).[\u80cc\u5305\u95ee\u9898]\u6709\u4e00\u4e2a\u80cc\u5305\uff0c\u80cc\u5305\u5bb9\u91cf\u662fM=150\u3002\u67097\u4e2a\u7269\u54c1\uff0c\u7269\u54c1\u53ef\u4ee5\u5206\u5272\u6210\u4efb\u610f\u5927\u5c0f\u3002
\u8981\u6c42\u5c3d\u53ef\u80fd\u8ba9\u88c5\u5165\u80cc\u5305\u4e2d\u7684\u7269\u54c1\u603b\u4ef7\u503c\u6700\u5927\uff0c\u4f46\u4e0d\u80fd\u8d85\u8fc7\u603b\u5bb9\u91cf\u3002

\u7269\u54c1 A B C D E F G
\u91cd\u91cf 35 30 60 50 40 10 25
\u4ef7\u503c 10 40 30 50 35 40 30

\u5206\u6790\uff1a
\u76ee\u6807\u51fd\u6570\uff1a \u2211pi\u6700\u5927
\u7ea6\u675f\u6761\u4ef6\u662f\u88c5\u5165\u7684\u7269\u54c1\u603b\u91cd\u91cf\u4e0d\u8d85\u8fc7\u80cc\u5305\u5bb9\u91cf\uff1a\u2211wi<=M( M=150)
\uff081\uff09\u6839\u636e\u8d2a\u5fc3\u7684\u7b56\u7565\uff0c\u6bcf\u6b21\u6311\u9009\u4ef7\u503c\u6700\u5927\u7684\u7269\u54c1\u88c5\u5165\u80cc\u5305\uff0c\u5f97\u5230\u7684\u7ed3\u679c\u662f\u5426\u6700\u4f18\uff1f
\uff082\uff09\u6bcf\u6b21\u6311\u9009\u6240\u5360\u7a7a\u95f4\u6700\u5c0f\u7684\u7269\u54c1\u88c5\u5165\u662f\u5426\u80fd\u5f97\u5230\u6700\u4f18\u89e3\uff1f
\uff083\uff09\u6bcf\u6b21\u9009\u53d6\u5355\u4f4d\u5bb9\u91cf\u4ef7\u503c\u6700\u5927\u7684\u7269\u54c1\uff0c\u6210\u4e3a\u89e3\u672c\u9898\u7684\u7b56\u7565\u3002

(\u73af\u5883:c++)
#include
#define max 100 //\u6700\u591a\u7269\u54c1\u6570
void sort (int n,float a[max],float b[max]) //\u6309\u4ef7\u503c\u5bc6\u5ea6\u6392\u5e8f
{
int j,h,k;
float t1,t2,t3,c[max];
for(k=1;k<=n;k++)
c[k]=a[k]/b[k];
for(h=1;h<n;h++)
for(j=1;j<=n-h;j++)
if(c[j]<c[j+1])
{t1=a[j];a[j]=a[j+1];a[j+1]=t1;
t2=b[j];b[j]=b[j+1];b[j+1]=t2;
t3=c[j];c[j]=c[j+1];c[j+1]=t3;
}
}
void knapsack(int n,float limitw,float v[max],float w[max],int x[max])
{float c1; //c1\u4e3a\u80cc\u5305\u5269\u4f59\u53ef\u88c5\u8f7d\u91cd\u91cf
int i;
sort(n,v,w); //\u7269\u54c1\u6309\u4ef7\u503c\u5bc6\u5ea6\u6392\u5e8f
c1=limitw;
for(i=1;i<=n;i++)
{
if(w[i]>c1)break;
x[i]=1; //x[i]\u4e3a1\u65f6\uff0c\u7269\u54c1i\u5728\u89e3\u4e2d
c1=c1-w[i];
}
}
void main()
{int n,i,x[max];
float v[max],w[max],totalv=0,totalw=0,limitw;
cout<<"\u8bf7\u8f93\u5165n\u548climitw:";
cin>>n >>limitw;
for(i=1;i<=n;i++)
x[i]=0; //\u7269\u54c1\u9009\u62e9\u60c5\u51b5\u8868\u521d\u59cb\u5316\u4e3a0
cout<<"\u8bf7\u4f9d\u6b21\u8f93\u5165\u7269\u54c1\u7684\u4ef7\u503c\uff1a"<<endl;
for(i=1;i<=n;i++)
cin>>v[i];
cout<<endl;
cout<<"\u8bf7\u4f9d\u6b21\u8f93\u5165\u7269\u54c1\u7684\u91cd\u91cf\uff1a"<<endl;
for(i=1;i<=n;i++)
cin>>w[i];
cout<<endl;
knapsack (n,limitw,v,w,x);
cout<<"the selection is:";
for(i=1;i<=n;i++)
{
cout<<x[i];
if(x[i]==1)
totalw=totalw+w[i];
}
cout<<endl;
cout<<"\u80cc\u5305\u7684\u603b\u91cd\u91cf\u4e3a\uff1a"<<totalw<<endl; //\u80cc\u5305\u6240\u88c5\u8f7d\u603b\u91cd\u91cf
cout<<"\u80cc\u5305\u7684\u603b\u4ef7\u503c\u4e3a\uff1a"<<totalv<<endl; //\u80cc\u5305\u7684\u603b\u4ef7\u503c
}
\u4e09.\u56de\u6eaf\u7b97\u6cd5\u6c42\u89e30-1\u80cc\u5305\u95ee\u9898
1.0-l\u80cc\u5305\u95ee\u9898\u662f\u5b50\u96c6\u9009\u53d6\u95ee\u9898\u3002
\u4e00\u822c\u60c5\u51b5\u4e0b\uff0c0-1\u80cc\u5305\u95ee\u9898\u662fNP\u96be\u9898\u30020-1\u80cc\u5305
\u95ee\u9898\u7684\u89e3\u7a7a\u95f4\u53ef\u7528\u5b50\u96c6\u6811\u8868\u793a\u3002\u89e30-1\u80cc\u5305\u95ee\u9898\u7684\u56de\u6eaf\u6cd5\u4e0e\u88c5\u8f7d\u95ee\u9898\u7684\u56de\u6eaf\u6cd5\u5341\u5206\u7c7b
\u4f3c\u3002\u5728\u641c\u7d22\u89e3\u7a7a\u95f4\u6811\u65f6\uff0c\u53ea\u8981\u5176\u5de6\u513f\u5b50\u7ed3\u70b9\u662f\u4e00\u4e2a\u53ef\u884c\u7ed3\u70b9\uff0c\u641c\u7d22\u5c31\u8fdb\u5165\u5176\u5de6\u5b50\u6811\u3002\u5f53
\u53f3\u5b50\u6811\u6709\u53ef\u80fd\u5305\u542b\u6700\u4f18\u89e3\u65f6\u624d\u8fdb\u5165\u53f3\u5b50\u6811\u641c\u7d22\u3002\u5426\u5219\u5c06\u53f3\u5b50\u6811\u526a\u53bb\u3002\u8bber\u662f\u5f53\u524d\u5269\u4f59
\u7269\u54c1\u4ef7\u503c\u603b\u548c\uff1bcp\u662f\u5f53\u524d\u4ef7\u503c\uff1bbestp\u662f\u5f53\u524d\u6700\u4f18\u4ef7\u503c\u3002\u5f53cp+r\u2264bestp\u65f6\uff0c\u53ef\u526a\u53bb\u53f3
\u5b50\u6811\u3002\u8ba1\u7b97\u53f3\u5b50\u6811\u4e2d\u89e3\u7684\u4e0a\u754c\u7684\u66f4\u597d\u65b9\u6cd5\u662f\u5c06\u5269\u4f59\u7269\u54c1\u4f9d\u5176\u5355\u4f4d\u91cd\u91cf\u4ef7\u503c\u6392\u5e8f\uff0c\u7136\u540e
\u4f9d\u6b21\u88c5\u5165\u7269\u54c1\uff0c\u76f4\u81f3\u88c5\u4e0d\u4e0b\u65f6\uff0c\u518d\u88c5\u5165\u8be5\u7269\u54c1\u7684\u4e00\u90e8\u5206\u800c\u88c5\u6ee1\u80cc\u5305\u3002\u7531\u6b64\u5f97\u5230\u7684\u4ef7\u503c\u662f
\u53f3\u5b50\u6811\u4e2d\u89e3\u7684\u4e0a\u754c\u3002
2.\u89e3\u51b3\u529e\u6cd5\u601d\u8def\uff1a
\u4e3a\u4e86\u4fbf\u4e8e\u8ba1\u7b97\u4e0a\u754c\uff0c\u53ef\u5148\u5c06\u7269\u54c1\u4f9d\u5176\u5355\u4f4d\u91cd\u91cf\u4ef7\u503c\u4ece\u5927\u5230\u5c0f\u6392\u5e8f\uff0c\u6b64\u540e\u53ea\u8981\u987a\u5e8f\u8003
\u5bdf\u5404\u7269\u54c1\u5373\u53ef\u3002\u5728\u5b9e\u73b0\u65f6\uff0c\u7531bound\u8ba1\u7b97\u5f53\u524d\u7ed3\u70b9\u5904\u7684\u4e0a\u754c\u3002\u5728\u641c\u7d22\u89e3\u7a7a\u95f4\u6811\u65f6\uff0c\u53ea\u8981\u5176\u5de6\u513f\u5b50\u8282\u70b9\u662f\u4e00\u4e2a\u53ef\u884c\u7ed3\u70b9\uff0c\u641c\u7d22\u5c31\u8fdb\u5165\u5de6\u5b50\u6811\uff0c\u5728\u53f3\u5b50\u6811\u4e2d\u6709\u53ef\u80fd\u5305\u542b\u6700\u4f18\u89e3\u662f\u624d\u8fdb\u5165\u53f3\u5b50\u6811\u641c\u7d22\u3002\u5426\u5219\u5c06\u53f3\u5b50\u6811\u526a\u53bb\u3002

\u56de\u6eaf\u6cd5\u662f\u4e00\u4e2a\u65e2\u5e26\u6709\u7cfb\u7edf\u6027\u53c8\u5e26\u6709\u8df3\u8dc3\u6027\u7684\u7684\u641c\u7d22\u7b97\u6cd5\u3002\u5b83\u5728\u5305\u542b\u95ee\u9898\u7684\u6240\u6709\u89e3\u7684\u89e3\u7a7a\u95f4\u6811\u4e2d\uff0c\u6309\u7167\u6df1\u5ea6\u4f18\u5148\u7684\u7b56\u7565\uff0c\u4ece\u6839\u7ed3\u70b9\u51fa\u53d1\u641c\u7d22\u89e3\u7a7a\u95f4\u6811\u3002\u7b97\u6cd5\u641c\u7d22\u81f3\u89e3\u7a7a\u95f4\u6811\u7684\u4efb\u4e00\u7ed3\u70b9\u65f6\uff0c\u603b\u662f\u5148\u5224\u65ad\u8be5\u7ed3\u70b9\u662f\u5426\u80af\u5b9a\u4e0d\u5305\u542b\u95ee\u9898\u7684\u89e3\u3002\u5982\u679c\u80af\u5b9a\u4e0d\u5305\u542b\uff0c\u5219\u8df3\u8fc7\u5bf9\u4ee5\u8be5\u7ed3\u70b9\u4e3a\u6839\u7684\u5b50\u6811\u7684\u7cfb\u7edf\u641c\u7d22\uff0c\u9010\u5c42\u5411\u5176\u7956\u5148\u7ed3\u70b9\u56de\u6eaf\u3002\u5426\u5219\uff0c\u8fdb\u5165\u8be5\u5b50\u6811\uff0c\u7ee7\u7eed\u6309\u6df1\u5ea6\u4f18\u5148\u7684\u7b56\u7565\u8fdb\u884c\u641c\u7d22\u3002\u56de\u6eaf\u6cd5\u5728\u7528\u6765\u6c42\u95ee\u9898\u7684\u6240\u6709\u89e3\u65f6\uff0c\u8981\u56de\u6eaf\u5230\u6839\uff0c\u4e14\u6839\u7ed3\u70b9\u7684\u6240\u6709\u5b50\u6811\u90fd\u5df2\u88ab\u641c\u7d22\u904d\u624d\u7ed3\u675f\u3002\u800c\u56de\u6eaf\u6cd5\u5728\u7528\u6765\u6c42\u95ee\u9898\u7684\u4efb\u4e00\u89e3\u65f6\uff0c\u53ea\u8981\u641c\u7d22\u5230\u95ee\u9898\u7684\u4e00\u4e2a\u89e3\u5c31\u53ef\u4ee5\u7ed3\u675f\u3002\u8fd9\u79cd\u4ee5\u6df1\u5ea6\u4f18\u5148\u7684\u65b9\u5f0f\u7cfb\u7edf\u5730\u641c\u7d22\u95ee\u9898\u7684\u89e3\u7684\u7b97\u6cd5\u79f0\u4e3a\u56de\u6eaf\u6cd5\uff0c\u5b83\u9002\u7528\u4e8e\u89e3\u4e00\u4e9b\u7ec4\u5408\u6570\u8f83\u5927\u7684\u95ee\u9898\u3002
2.\u7b97\u6cd5\u6846\u67b6\uff1a
a.\u95ee\u9898\u7684\u89e3\u7a7a\u95f4\uff1a\u5e94\u7528\u56de\u6eaf\u6cd5\u89e3\u95ee\u9898\u65f6\uff0c\u9996\u5148\u5e94\u660e\u786e\u5b9a\u4e49\u95ee\u9898\u7684\u89e3\u7a7a\u95f4\u3002\u95ee\u9898\u7684\u89e3\u7a7a\u95f4\u5e94\u5230\u5c11\u5305\u542b\u95ee\u9898\u7684\u4e00\u4e2a\uff08\u6700\u4f18\uff09\u89e3\u3002
b.\u56de\u6eaf\u6cd5\u7684\u57fa\u672c\u601d\u60f3\uff1a\u786e\u5b9a\u4e86\u89e3\u7a7a\u95f4\u7684\u7ec4\u7ec7\u7ed3\u6784\u540e\uff0c\u56de\u6eaf\u6cd5\u5c31\u4ece\u5f00\u59cb\u7ed3\u70b9\uff08\u6839\u7ed3\u70b9\uff09\u51fa\u53d1\uff0c\u4ee5\u6df1\u5ea6\u4f18\u5148\u7684\u65b9\u5f0f\u641c\u7d22\u6574\u4e2a\u89e3\u7a7a\u95f4\u3002\u8fd9\u4e2a\u5f00\u59cb\u7ed3\u70b9\u5c31\u6210\u4e3a\u4e00\u4e2a\u6d3b\u7ed3\u70b9\uff0c\u540c\u65f6\u4e5f\u6210\u4e3a\u5f53\u524d\u7684\u6269\u5c55\u7ed3\u70b9\u3002\u5728\u5f53\u524d\u7684\u6269\u5c55\u7ed3\u70b9\u5904\uff0c\u641c\u7d22\u5411\u7eb5\u6df1\u65b9\u5411\u79fb\u81f3\u4e00\u4e2a\u65b0\u7ed3\u70b9\u3002\u8fd9\u4e2a\u65b0\u7ed3\u70b9\u5c31\u6210\u4e3a\u4e00\u4e2a\u65b0\u7684\u6d3b\u7ed3\u70b9\uff0c\u5e76\u6210\u4e3a\u5f53\u524d\u6269\u5c55\u7ed3\u70b9\u3002\u5982\u679c\u5728\u5f53\u524d\u7684\u6269\u5c55\u7ed3\u70b9\u5904\u4e0d\u80fd\u518d\u5411\u7eb5\u6df1\u65b9\u5411\u79fb\u52a8\uff0c\u5219\u5f53\u524d\u6269\u5c55\u7ed3\u70b9\u5c31\u6210\u4e3a\u6b7b\u7ed3\u70b9\u3002\u6362\u53e5\u8bdd\u8bf4\uff0c\u8fd9\u4e2a\u7ed3\u70b9\u4e0d\u518d\u662f\u4e00\u4e2a\u6d3b\u7ed3\u70b9\u3002\u6b64\u65f6\uff0c\u5e94\u5f80\u56de\u79fb\u52a8\uff08\u56de\u6eaf\uff09\u81f3\u6700\u8fd1\u7684\u4e00\u4e2a\u6d3b\u7ed3\u70b9\u5904\uff0c\u5e76\u4f7f\u8fd9\u4e2a\u6d3b\u7ed3\u70b9\u6210\u4e3a\u5f53\u524d\u7684\u6269\u5c55\u7ed3\u70b9\u3002\u56de\u6eaf\u6cd5\u5373\u4ee5\u8fd9\u79cd\u5de5\u4f5c\u65b9\u5f0f\u9012\u5f52\u5730\u5728\u89e3\u7a7a\u95f4\u4e2d\u641c\u7d22\uff0c\u76f4\u81f3\u627e\u5230\u6240\u8981\u6c42\u7684\u89e3\u6216\u89e3\u7a7a\u95f4\u4e2d\u5df2\u6ca1\u6709\u6d3b\u7ed3\u70b9\u65f6\u4e3a\u6b62\u3002
3.\u8fd0\u7528\u56de\u6eaf\u6cd5\u89e3\u9898\u901a\u5e38\u5305\u542b\u4ee5\u4e0b\u4e09\u4e2a\u6b65\u9aa4\uff1a
a.\u9488\u5bf9\u6240\u7ed9\u95ee\u9898\uff0c\u5b9a\u4e49\u95ee\u9898\u7684\u89e3\u7a7a\u95f4\uff1b
b.\u786e\u5b9a\u6613\u4e8e\u641c\u7d22\u7684\u89e3\u7a7a\u95f4\u7ed3\u6784\uff1b
c.\u4ee5\u6df1\u5ea6\u4f18\u5148\u7684\u65b9\u5f0f\u641c\u7d22\u89e3\u7a7a\u95f4\uff0c\u5e76\u4e14\u5728\u641c\u7d22\u8fc7\u7a0b\u4e2d\u7528\u526a\u679d\u51fd\u6570\u907f\u514d\u65e0\u6548\u641c\u7d22\uff1b
#include

using namespace std;


class Knap
{
friend int Knapsack(int p[],int w[],int c,int n );

public:
void print()
{

for(int m=1;m<=n;m++)
{
cout<<bestx[m]<<" ";
}
cout<<endl;
};

private:
int Bound(int i);
void Backtrack(int i);

int c;//\u80cc\u5305\u5bb9\u91cf
int n; //\u7269\u54c1\u6570
int *w;//\u7269\u54c1\u91cd\u91cf\u6570\u7ec4
int *p;//\u7269\u54c1\u4ef7\u503c\u6570\u7ec4
int cw;//\u5f53\u524d\u91cd\u91cf
int cp;//\u5f53\u524d\u4ef7\u503c
int bestp;//\u5f53\u524d\u6700\u4f18\u503c
int *bestx;//\u5f53\u524d\u6700\u4f18\u89e3
int *x;//\u5f53\u524d\u89e3

};


int Knap::Bound(int i)
{
//\u8ba1\u7b97\u4e0a\u754c
int cleft=c-cw;//\u5269\u4f59\u5bb9\u91cf
int b=cp;
//\u4ee5\u7269\u54c1\u5355\u4f4d\u91cd\u91cf\u4ef7\u503c\u9012\u51cf\u5e8f\u88c5\u5165\u7269\u54c1
while(i<=n&&w[i]<=cleft)
{
cleft-=w[i];
b+=p[i];
i++;
}
//\u88c5\u6ee1\u80cc\u5305
if(i<=n)
b+=p[i]/w[i]*cleft;
return b;
}


void Knap::Backtrack(int i)
{
if(i>n)
{
if(bestp<cp)
{
for(int j=1;j<=n;j++)
bestx[j]=x[j];
bestp=cp;
}
return;
}
if(cw+w[i]<=c) //\u641c\u7d22\u5de6\u5b50\u6811
{
x[i]=1;
cw+=w[i];
cp+=p[i];
Backtrack(i+1);
cw-=w[i];
cp-=p[i];
}
if(Bound(i+1)>bestp)//\u641c\u7d22\u53f3\u5b50\u6811
{
x[i]=0;
Backtrack(i+1);
}

}


class Object
{
friend int Knapsack(int p[],int w[],int c,int n);
public:
int operator<=(Object a)const
{
return (d>=a.d);
}

private:
int ID;
float d;
};


int Knapsack(int p[],int w[],int c,int n)
{
//\u4e3aKnap::Backtrack\u521d\u59cb\u5316
int W=0;
int P=0;
int i=1;
Object *Q=new Object[n];
for(i=1;i<=n;i++)
{
Q[i-1].ID=i;
Q[i-1].d=1.0*p[i]/w[i];
P+=p[i];
W+=w[i];
}
if(W<=c)
return P;//\u88c5\u5165\u6240\u6709\u7269\u54c1
//\u4f9d\u7269\u54c1\u5355\u4f4d\u91cd\u91cf\u6392\u5e8f
float f;
for( i=0;i<n;i++)
for(int j=i;j<n;j++)
{
if(Q[i].d<Q[j].d)
{
f=Q[i].d;
Q[i].d=Q[j].d;
Q[j].d=f;
}


}

Knap K;
K.p = new int[n+1];
K.w = new int[n+1];
K.x = new int[n+1];
K.bestx = new int[n+1];
K.x[0]=0;
K.bestx[0]=0;
for( i=1;i<=n;i++)
{
K.p[i]=p[Q[i-1].ID];
K.w[i]=w[Q[i-1].ID];
}
K.cp=0;
K.cw=0;
K.c=c;
K.n=n;
K.bestp=0;
//\u56de\u6eaf\u641c\u7d22
K.Backtrack(1);
K.print();
delete [] Q;
delete [] K.w;
delete [] K.p;
return K.bestp;

}

void main()
{
int *p;
int *w;
int c=0;
int n=0;
int i=0;
char k;
cout<<"0-1\u80cc\u5305\u95ee\u9898\u2014\u2014\u56de\u6eaf\u6cd5 "<<endl;
cout<<" by zbqplayer "<<endl;
while(k)
{
cout<<"\u8bf7\u8f93\u5165\u80cc\u5305\u5bb9\u91cf(c)\uff1a"<<endl;
cin>>c;
cout<<"\u8bf7\u8f93\u5165\u7269\u54c1\u7684\u4e2a\u6570(n)\uff1a"<<endl;
cin>>n;
p=new int[n+1];
w=new int[n+1];
p[0]=0;
w[0]=0;

cout<<"\u8bf7\u8f93\u5165\u7269\u54c1\u7684\u4ef7\u503c(p)\uff1a"<<endl;
for(i=1;i<=n;i++)
cin>>p[i];

cout<<"\u8bf7\u8f93\u5165\u7269\u54c1\u7684\u91cd\u91cf(w)\uff1a"<<endl;
for(i=1;i<=n;i++)
cin>>w[i];

cout<<"\u6700\u4f18\u89e3\u4e3a(bestx)\uff1a"<<endl;
cout<<"\u6700\u4f18\u503c\u4e3a(bestp)\uff1a"<<endl;
cout<<Knapsack(p,w,c,n)<<endl;
cout<<"[s] \u91cd\u65b0\u5f00\u59cb"<<endl;
cout<<"[q] \u9000\u51fa"<<endl;
cin>>k;
}
\u56db.\u5206\u652f\u9650\u754c\u6cd5\u6c42\u89e30-1\u80cc\u5305\u95ee\u9898
1.\u95ee\u9898\u63cf\u8ff0\uff1a\u5df2\u77e5\u6709N\u4e2a\u7269\u54c1\u548c\u4e00\u4e2a\u53ef\u4ee5\u5bb9\u7eb3M\u91cd\u91cf\u7684\u80cc\u5305\uff0c\u6bcf\u79cd\u7269\u54c1I\u7684\u91cd\u91cf\u4e3aWEIGHT\uff0c\u4e00\u4e2a\u53ea\u80fd\u5168\u653e\u5165\u6216\u8005\u4e0d\u653e\u5165\uff0c\u6c42\u89e3\u5982\u4f55\u653e\u5165\u7269\u54c1\uff0c\u53ef\u4ee5\u4f7f\u80cc\u5305\u91cc\u7684\u7269\u54c1\u7684\u603b\u6548\u76ca\u6700\u5927\u3002

2.\u8bbe\u8ba1\u601d\u60f3\u4e0e\u5206\u6790\uff1a\u5bf9\u7269\u54c1\u7684\u9009\u53d6\u4e0e\u5426\u6784\u6210\u4e00\u68f5\u89e3\u6811\uff0c\u5de6\u5b50\u6811\u8868\u793a\u4e0d\u88c5\u5165\uff0c\u53f3\u8868\u793a\u88c5\u5165\uff0c\u901a\u8fc7\u68c0\u7d22\u95ee\u9898\u7684\u89e3\u6811\u5f97\u51fa\u6700\u4f18\u89e3\uff0c\u5e76\u7528\u7ed3\u70b9\u4e0a\u754c\u6740\u6b7b\u4e0d\u7b26\u5408\u8981\u6c42\u7684\u7ed3\u70b9\u3002

#include

struct good
{
int weight;
int benefit;
int flag;//\u662f\u5426\u53ef\u4ee5\u88c5\u5165\u6807\u8bb0
};

int number=0;//\u7269\u54c1\u6570\u91cf
int upbound=0;
int curp=0, curw=0;//\u5f53\u524d\u6548\u76ca\u503c\u4e0e\u91cd\u91cf
int maxweight=0;
good *bag=NULL;

void Init_good()
{
bag=new good [number];

for(int i=0; i<number; i++)
{
cout<<"\u8bf7\u8f93\u5165\u7b2c\u4ef6"<<i+1<<"\u7269\u54c1\u7684\u91cd\u91cf:";
cin>>bag[i].weight;
cout<<"\u8bf7\u8f93\u5165\u7b2c\u4ef6"<<i+1<<"\u7269\u54c1\u7684\u6548\u76ca:";
cin>>bag[i].benefit;
bag[i].flag=0;//\u521d\u59cb\u6807\u5fd7\u4e3a\u4e0d\u88c5\u5165\u80cc\u5305
cout<<endl;
}

}

int getbound(int num, int *bound_u)//\u8fd4\u56de\u672c\u7ed3\u70b9\u7684c\u9650\u754c\u548cu\u9650\u754c
{
for(int w=curw, p=curp; num<number && (w+bag[num].weight)<=maxweight; num++)
{
w=w+bag[num].weight;
p=w+bag[num].benefit;
}

*bound_u=p+bag[num].benefit;
return ( p+bag[num].benefit*((maxweight-w)/bag[num].weight) );
}

void LCbag()
{
int bound_u=0, bound_c=0;//\u5f53\u524d\u7ed3\u70b9\u7684c\u9650\u754c\u548cu\u9650\u754c

for(int i=0; i<number; i++)//\u9010\u5c42\u904d\u5386\u89e3\u6811\u51b3\u5b9a\u662f\u5426\u88c5\u5165\u5404\u4e2a\u7269\u54c1
{
if( ( bound_c=getbound(i+1, &bound_u) )>upbound )//\u904d\u5386\u5de6\u5b50\u6811
upbound=bound_u;//\u66f4\u6539\u5df2\u6709u\u9650\u754c,\u4e0d\u66f4\u6539\u6807\u5fd7

if( getbound(i, &bound_u)>bound_c )//\u904d\u5386\u53f3\u5b50\u6811
//\u82e5\u88c5\u5165\uff0c\u5224\u65ad\u53f3\u5b50\u6811\u7684c\u9650\u754c\u662f\u5426\u5927\u4e8e\u5de6\u5b50\u6811\u6839\u7684c\u9650\u754c\uff0c\u662f\u5219\u88c5\u5165
{
upbound=bound_u;//\u66f4\u6539\u5df2\u6709u\u9650\u754c
curp=curp+bag[i].benefit;
curw=curw+bag[i].weight;//\u4ece\u5df2\u6709\u91cd\u91cf\u548c\u6548\u76ca\u52a0\u4e0a\u65b0\u7269\u54c1
bag[i].flag=1;//\u6807\u8bb0\u4e3a\u88c5\u5165
}
}

}

void Display()
{

cout<<"\u53ef\u4ee5\u653e\u5165\u80cc\u5305\u7684\u7269\u54c1\u7684\u7f16\u53f7\u4e3a\uff1a";
for(int i=0; i<number; i++)
if(bag[i].flag>0)
cout<<i+1<<" ";
cout<<endl;
delete []bag;
}

的动态规划0-1背包问题
/ ************************************* *********************************** /
/ * 0-1背包问题: /> / * n种物品和一个背包
/ *项目无线网络,我的体重VI的价值
/ *背包的容量为c
/ *应该如何选择项目加载背包使物品装入背包
/ *总最看重的吗?
/ *注:选择项目装入背包的物品我只有两个选择
/ *加载或不装入背包。我不能加载多个
/ *不能只是加载的项目。
/ *

/ * 1。 0-1背包问题的形式化描述:
/ * C> 0,无线网络> 0,六0,0 <= I <= N要求找到一个n
/ * 0-1向量( X1,X2,...,XN),使得:
/ *最大sum_ {= 1到n}(ⅵ*ⅹⅰ),并满足下面的约束:
/ *(1)sum_ {I = 1到n}(WI *十一)<= C

/ *(2)十一∈{0,1} 1 <= I <= N
/ * ...... / a>
/ * 2。
/ * 0-1背包问题0-1背包问题的解决有重叠的性质的最优子结构性质和子适合
/ *使用动态规划方法解决
/

/ * 2.1最优子结构性质
/ *设置(Y1,Y2,...,yn)的0-1背包问题,最佳的解决方案,它肯定
> / *结论(Y2,Y3,...,yn)的是下面的子问题的最优解:
/ *最大sum_ {i = 2到n}(VI *十一)
/ *(1)sum_ {i = 2到N}(WI *十一)<= C - W1 * Y1
/ *(2)十一∈{0,1},2 <= I <= N /> / *否则,子问题有一个最优的解决方案(Z2,Z3,...,zn的),
/ *(Y2,Y3,...,yn)的比的最优解。然后:
/ * sum_ {i = 2到n}(VI *子)sum_ {i = 2到n}(VI *易)
/ *,W1 * Y1 + sum_ {i = 2到n}(WI *子)<= C
/ *进一步指出:
/ * V1 * Y1 + sum_ {i = 2到n}(VI *子)> sum_ {i = 1 N}(VI *易)
/ * W1 * Y1 + sum_ {i = 2到n}(无线ZI)<= C
/ *这说明:(Y1,Z2,Z3, ,...,ZN)0-1背包问题是一个更好的解决方案,然后
/ *说明(Y1,Y2,...,yn)的是不矛盾的前提下,最佳的解决方案,因此,最好
/ *建立的子结构的性质。
/ *
/ * 2.2子重叠性质
/ *设置给定的0-1背包问题的子P(I,J):
/ *最大sum_ {K = I N}(VK * XK)
/ *(1)sum_ {K = I为n}(周* XK)<= J
/ *(2)XK∈{0,1}, I <= K <= N

/ * P(I,J)的问题是背包容量为j,可选的项目我,我+1,...,n的子的问题
/ *设M(I,J)是子P(I,J),最大总价值的最优值。在最佳
/ *子结构的性质可以创建递归米的(I,J):
/ *。递归的初始M(N,J)
/ * / /背包容量为j,可选的项目只有n,,如果背包容量j是大于项目?
/ / /重量直接加载;无法加载。
/ * M(N,J)= VN,J> = WN
/ * M(N,J)= 0,0 <= J <WN
/ * B。递推公式M(I,J)
/ * / /背包容量?,可选的项目我,我+1,...,N
/ * / /背包容量J <无线网络,后来干脆不安装成的文章:
/ * M(I,J)= M(I +1,J),0 <= j的无线
/ * / / J> =无线网络,文章中,我加载的项目之间我选择
/ *没有我的最佳值:M(i +1,j)的安装项目
/ *加载项,我的最优值:M(i +1, J-WI)+ VI

/ *:
/ * M(I,J)= {M(I +1,J),M(I +1,J-无线)+ VI},J> =无线
/ *
/ ***************************** ******************************************* /
/>

定义MAX(A,B)(((一)(二))(A)(B))
定义分(A,B)( ((A)(B))(A)(B))
模板
无效背包(类型* V,INT * W,C,N,类型* * M)
{
/ /递归的初始条件
诠释JMAX = MIN(W [N] - 1,C);
(J = 0;? <= JMAX J + +){
米[N] [J] = 0;
}

(J = W [N]; <= C; J + +){
M [N] [J] = V [N];
}

/ /我从2到n-1,分别为j> = wi和0 <= j的无线M(I,J)
(INT I = N-1; i> 1的,我 - ){
JMAX = MIN(W [我] - 1,C);
(J = 0; <= JMAX; J + +){
M [] [J] = M [+1] [J]; ...... />}
为(J = W [I]; <= C; + +){
M [] [J] = MAX(M [+1] [J], M [+1] [JW [我] + V [I]);
}
}

M [1] [C] = M [2] [ C];
(C> = W [1]){
M [1] [C] = MAX(M [1] [C],M [2] [CW [1]] + V [1]);
}

}

模板
无效回溯(类型**米,INT * W, INT整数N,C,* X)
{
为(int i = 1; <N,我+ +){
(M [] [C] == M [+1] [C])×[I] = 0;
其他{
×[我] = 1;
C - = W [I];
}
}
×[N] =(M [N] [C])? 1:0;
}

:(INT ARGC,字符*的argv [])
{
廉政n = 5;
>诠释W [6] = {-1,2,2,6,5,4};
INT V [6] = {-1,6,3,5,4,6};
>诠释C = 10;

** ppm =百新诠释[N +1];
(INT I = 0; I <N +1,我+ +){
PPM [] =新的int [C +1];
}

诠释x [6];

背包的(V ,W,C,N,PPM);
:回溯的(ppm时,W,C,N,X);

返回0;}

贪心算法求解0-1背包问题
的基本思路?贪婪的方法:
逐渐接近目标,给予尽可能快地获得一个更好的解决方案在地上 - 从最初的解决方案的问题。当达到一个算法的步骤不能继续向前走,算法停止。
的算法问题:
1)不能保证获得最终的解决方案;
2)不能被用来谋求最大或最小解决方案;
3)寻求满足一定的约束条件,可行的解决方案。

算法:
开始从最初的解决方案的问题;

获得可行的解决方案提出了总体目标,而移动元素的解决方案;
所有解决方案组件组合成一个可行的解决方案的问题;
2实例分析

1)。背包问题,背包,背包容量是M = 150。 7物品时,物品可分为任意的尺寸。的
要求尽可能多,以使总价值的物品装入背包,但不能超过总容量。

的项ABCDEFG
权重值35 30 60 50 40 10 25
10 40 30 50 35 40 30


目标函数: σpi,
约束装入的物品的总重量计不超过背包容量:Σwi<= M(M = 150)
(1)根据贪心策略,每个选定的值最大的装载物品的背包,得到的结果是最优的?
(2)每次选择负载最小的项目可以得到最佳的解决方案吗?
(3)每次你选择一个单位容量的最有价值的项目,解决这个问题的战略。

(环境:C + +)
包括与ltiostream.h的>
#定义最大100 / /最大数量的项目
无效排序(N,飘起了[MAX],持股量B [MAX])/ /密度值的排序
{
整数J,H,和K;
持股量T1,T2, T3,C [MAX];
(k = 1时,K <= N,K + +)
?[K] = [K] / B [K]
(H = 1,H <N,H + +)
(J = 1,J <= NH J + +)
(C [J] C [j +1]中)
{T1 = A [J],A [J] = A [j +1]中的[J +1] = T1;
T2 = B [J]; B [J] = B [j +1]中,B [J +1] = T2;
T3 = C [J] C [J] = C [j +1]中,C [J +1] = T3;
>}
}
无效背包(INT N,持股量limitw,浮动V [MAX],浮瓦特[MAX],诠释x [MAX])
{浮动C1 / / C1背包剩余的负载重量

排序(N,V,W)/ /排序价值密度的
C1 = limitw
为(i = 1我=我+ +)
{
(W [I]> C1)打破;
×[我] = 1; / / x [我]的文章我了解
C1 = C1-W [I];
}
}
无效的主要()
{N,I,X [MAX];
>浮法V [MAX],W [MAX],totalv = 0,totalw = 0,limitw;
法院<<“请输入n和limitw:”
CIN >> N >> limitw
为(i = 1; <=我+ +)
×[我] = 0; / /初始化为0的产品选择表
法院<<“请进入价值的项目依次:“<< endl;
(i = 1; <=我+ +)
CIN >> V [I];
法院<; <endl;
法院<<“请输入的重量反过来的项目:”<< endl;
(i = 1; <= N; + +)
CIN >> W [I];
法院<< endl;
背包(N,limitw,V,W,X);
法院<<“选择的是:”
>(i = 1; <= N; + +)
{
法院<< X [I];
(X [I] == 1) totalw = totalw + W [I];
}
法院<< endl;
法院<<“总重量的背包:<<totalw << endl; / /背包装载总重量
法院<<“总价值的背包为:的”“totalv << endl; / /背包的总价值
}
三回溯算法0 -1背包问题
1.0-l背包问题是选定的子集的问题。
一般,0-1背包问题是NP-难的。
0-1背包解决方案的可用空间的一个子集树。
喜欢回溯0-1背包问题装载问题的回溯是非常一流。搜索解空间树的搜索,只要其左子节点是一个可行的节点,进入其左子树。
右子树可能只包含右子树搜索的最佳解决方案,否则,切右子树。设r其余
商品的价值的总和;阴极保护电流值;当前最优值bestp。
子树可以削权当CP + R≤bestp。在右子树的上限解的计算是更好的方法,其余项目根据其单位重量排序的值
然后在打开加载项直到合适的,然后到的文章充满背包。得到的值是上限
右子树的解决方案。
2。解决方案的想法:
为了便于计算中上的第一个项目,每单位重量的约束根据降序自己的价值,并此后,只要测试的顺序
可以观察到的各种物品。在实现绑定在当前节点的上限计算。搜索解决方案空间树,只要其左子节点是一个可行的节点,搜索到的左子树,右子树可能含有进入右子树搜索的最佳解决方案之前,否则,切右子树。

回溯跳跃系统的搜索算法。它包含了所有问题的解空间树的解决方案,根据深度优先的策略,开始从根本上搜索解空间树。总是首先确定节点的解决方案不包含树算法搜索的解空间中的任何节点。当然不包含跳层和它的祖先节点的节点是根的子树的搜索回溯系统,否则,进入子树,继续深入搜索优先策略。回溯所有使用的解决方案,提出了一个问题,我们必须回头去根,根的子树搜索,直到结束。回溯,用乞求的问题,任何一个解决方案,只要搜索一个解决问题的办法可以结束了。被称为回溯的深度优先搜索算法的问题的解决方案,它适用于了解一些大量的组合。
2。算法框架: BR />问题的解空间:应用回溯法解决问题,首先应该明确的定义问题的解决空间问题的解空间中至少包含一个(最佳)的解决方案。
B。回溯基本思想是:以确定的组织结构对空间的理解,回溯从开始节点(根),深度优先搜索整个解空间的起始节点成为一个活结点,也是当前扩展节点。在当前的扩展结中,搜索到的深度方向移动到一个新的节点。新的节点是一个新的活结点,和目前的扩张节点,如果当前的交界处延长不能进一步移动,在深度方向上,那么当前的扩展结点死锁。换言之,此节点不再是一个活结点,在这一点上,应该向后移动(回溯)指向一个活结活结扩展节点停止回溯这样的递归地搜索解空间中工作,直到你找到解决方案或解决方案所需的空间活结点。
3。用回溯法解决问题通常是由以下三个步骤:
解决方案的空间一个给定的问题,定义问题;
B。确定易于搜索的解空间结构;
C。深度优先的方式搜索解空间,并在搜索过程中通过修剪功能,以避免无效搜索;
#包括

使用命名空间std;

类雷德克纳普
{
朋友整数背包(P [],诠释W [],诠释三,廉政n);


无效的print()
{

(M = 1,M <= N,M + +)
{
法院<< bestx [M] <<“”;
}
法院<< endl;
};

私人如下:
整数约束(int i)的;
无效回溯(int i)的;

诠释三;/ /背包容量
廉政n,/ /项目数
诠释* W / /菜单项权重数组
诠释* p ;/ /项目值的数组
诠释CW ;/ /重量
INT CP ;/ /电流值

诠释bestp ;/ /当前的最优值
诠释* bestx ;/ /当前最优解
诠释*

}

诠释雷德克纳普::绑定(一)
{
/ /计算x ;/ /当前的解决方案上限
诠释裂= C-CW ;/ /剩余容量
INT B = CP;
/ /项目单位重量价值递减的序列加载项
(I <= N &&瓦特[I] <=裂)
{
裂= W [I]
B + = [];
+ +;
>}
/ /填充背包
如果(i <= N)
B + = [I] / W [I] *裂;
回报B; ...... />}

无效雷德克纳普::回溯(一)
{
如果(I>)
{
( bestp <CP)
{
(J = 1; <= N; + +)
bestx [J] = [J]。
bestp = CP
}
回报;
}
(CW + W [I] <= C)/ /搜索左子树
{
×[我] = 1;
CW + = W [I];
CP + = [];
回溯 - 触控板(i +1);
CW - = W [I];
CP-= P [I];
}
如果“(绑定(i +1)> bestp)/ /搜索右子树
{ BR /> x [I] = 0;
回溯 - 触控板(i +1);
}

}

>类对象
{
朋友整数背包(P [],诠释W [],诠释三,廉政n);

内部操作符<=(对象A )常量
{
回报率(D> =广告);
}

私人如下:
整数ID;
浮动D; />};

整数背包(P [],诠释W [],诠释三,廉政n)
{
/ /为雷德克纳普: :Backtrack面板 - 触控板初始化
诠释W = 0;
诠释P = 0;
INT I = 1;
对象* Q =新的对象[N]; / a>(I = 1; <=我+ +)
{
Q [I-1]。ID =我;
Q [I-1],D = 1.0 * P [I] / W [I];
P + = [];
W + = W [I];
}
(W <= C)
回报P ;/ /加载的所有项目
/ /按产品单位重量排序
浮F;
为(i = 0; I <N;我+ +)
(整数J = I,J <N; J + +)
{
(Q [I],D <Q [J]。四)
{ BR /> F = Q [I],D;
Q [我],D = Q [J]。天;
Q [J]。D = F;
}
a>

}

雷德克纳普K;
KP =新的int [N +1];
KW =新的int [N +1];
KX = INT [N +1];
K.bestx =新的int [N +1];
KX [0] = 0;
K.bestx [0 ] = 0;
为(i = 1; <= N; + +)
{
KP [I] = P [Q [I-1]。ID]; BR /> KW [我] = W [Q [I-1]。ID];
}
K.cp = 0;
K.cw = 0;
KC = C;
KN = N;
K.bestp = 0;
/ /回溯搜索
K.Backtrack();
K.print的()
删除[] Q;
删除[]千瓦;
删除[] KP;
回报K.bestp;

}

无效的主要()
{
* P;
* W;
整数C = 0;
廉政n = 0;
> INT I = 0;
字符K表;
法院<<“0-1背包问题 - 回溯”<< endl;
法院<<“通过zbqplayer <而(K)
{
法院<<“请输入一个背包容量(C):”<< endl;
CIN >> C;
法院<<“请输入的项目数(n):“<< endl;
CIN >> N
P =新的int [N +1];
W =新的int [ +1];
P [0] = 0;
W [0] = 0;

法院<<“请输入一个值(P)的项目: << endl;
(i = 1; <= N; + +)
CIN >> P [I];

法院<<“请输入项目的重量(W):“<< endl;
为(i = 1; <=我+ +)
CIN >> W [I]
BR />法院<<“最佳的解决方案(bestx):”<< endl;
法院<<“最佳的值(bestp):”<< endl;
的cout <<背负的(P, W,C,N)<< endl;
法院<<“[S]重新启动”<< endl;
法院<<“[q]退出”<< endl; BR /> CIN >> K
}
四个分支定界法求解0-1背包问题
问题描述:已知的N项和一个背包可以容纳M个权重,权重我的体重,只认沽或不投入,解决如何把在背包中的物品的总收益的项目,可以使每一个项目。
2。设计的思考和分析:选择的项目,或不构成解决方案树,左子树不加载,正确的说,节点负载,以获得最佳的树检索问题的解决办法谢界杀不符合要求的节点。

包括

结构良好
{
诠释权重;
诠释的利益;
INT标志;/ /是否加载标签
};

整型数= 0 ;/ /项目数
诠释upbound = 0;
诠释curp = 0,curw = 0 ;/ /当前实际价值和重量
诠释maxweight = 0;
好包= NULL;

的虚空Init_good(){
>袋=新好[数字];

(INT I = 0;我数,我+ +)
{
法院<<“请输入第一块“。 ;“<< i +1 <<”重的项目:“
CIN >>袋[i]的重量;
法院<<”请输入片“<< i +1 << “项目的好处:”
CIN >>袋[I]。利益;
袋[I]。标志= 0 ;/ /初始标志没有被装入背包
法院<< endl;
}

}

诠释getbound(整数民, * bound_u)/ /返回节点c计和压力表ü
{
(瓦特= curw,P = curp,民数&&(W +包[NUM]重量)<= maxweight NUM + +)
{
W = W +包[NUM]。重量;
P = W +包[NUM]。效益;
}

* bound_u = P +包[NUM]。 ,效益;
回报(P +包[NUM]。利益((maxweight-W)/袋[NUM]。重量));
}

无效LCbag()
{
bound_u = 0,bound_c = 0 ;/ /当前节点c方向和u和约束

(i = 0;号码;我+ +层)/ /层穿越解决方案,树,以决定是否加载不同的项目
{(I +1,及bound_u))>
((bound_c =的getbound upbound)/ /遍历左子树
upbound = bound_u ;/ /改变有U计不会改变的标志
如果“(getbound(I&bound_u)> bound_c)/ /遍历右子树
/ /如果加载,以确定是否大于左子树的根的右子树C标尺C标尺负载
{的
upbound = bound_u ;/ /变化有u和约束
curp = curp +包[I]。利益的
curw = curw +包[i]的重量;/ /重量和效益从现有再加上新的项目
袋[I]。标志= 1; / /标记加载
}

}

}

显示()
{
/> COUT <<“项目可放置在背包里数:”;
(INT I = 0;我数,我+ +)
(袋[I]。标志> 0 )
法院<< i +1 <<“”;
法院<< endl;
删除[]袋;
}

扩展阅读:扫一扫一秒出答案 ... 分支定界法例题及详解 ... 扫一扫题目出答案 ... 免费作业拍照出答案 ... 0-1背包问题分支限界法 ... 分支定界法的剪支情况 ... 动态规划求0-1背包 ... 0 1背包问题回溯法 ... 队列分支限界法01背包图解 ...

本站交流只代表网友个人观点,与本站立场无关
欢迎反馈与建议,请联系电邮
2024© 车视网