粒子群算法之01背包问题（C语言实现）

参数说明

	Vi(t+1)=ωVi(t)+c1r1(pbesti−Xi(t))+c2r2(gbest−Xi(t))

根据粒子群算法速度迭代公式：
w：惯性因子
c1，c2：学习率
r1，r2：[0,1]之间的随机数
该式右端三项分别代表：自身惯性贡献、自身认知贡献以及群体经验贡献。

n：物品个数
m：鸟群大小
W：背包容量

算法注意事项

比起遗传算法，粒子群算法的超参数多一些，对于速度的大小，我选取的是[0,0.02]，然后就是三个贡献的系数选择要合理，如果c1,c2过大可能会导致陷入局部解。最大速度步长也需要控制一下，否则有可能导致一直走出定义域。鸟群的大小也会影响是否陷入局部最优解，鸟群越大越可能搜索到较优解。

代码

#include<bits/stdc++.h>
using namespace std;
#define maxn 105
#define maxm 1006
#define ll long long int
#define INF 0x3f3f3f3f
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
#define mem(a) memset(a,0,sizeof(a))
#define sqr(x) (x*x)
#define inf (ll)2e18+1
#define PI acos(-1)
#define mod 10007
#define db double
#define auto(i,x) for(int i=head[x];i;i=ed[i].nxt)
ll read(){
    ll x=0,f=1ll;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
     while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
     return f*x;
}
int n,m,W;
db xiaorand(db low,db up,db jin){
    db down=min(0.0,low);
    low-=down,up-=down;
    int l=(int)low*jin,r=(int)up*jin;
    int len=r-l+1;
    int num=rand()%len + l;
    return 1.0*num/jin + down;
}
struct node{int w,val;}a[maxn];
db w,c1,c2,r1,r2;
struct xiaoniao{
    db x[maxn],v[maxn],pbest[maxn];
    int value;
    bool operator < (const xiaoniao &rid)const{
        return value < rid.value;
    }
}gg[maxm],gbest;
int f(xiaoniao niao){
    int sumv=0,sumw=0;
    inc(i,1,n)if(niao.x[i]>0.5){sumv+=a[i].val,sumw+=a[i].w;}
    return sumw <= W ? sumv : -INF;
}
void init_gg(){
    gbest.value=-INF;
    inc(i,1,m)inc(j,1,n){
        gg[i].x[j]=xiaorand(0,1,1000);
        gg[i].v[j]=xiaorand(0,0.02,1000);
        gg[i].pbest[j]=gg[i].x[j];
        gg[i].value=f(gg[i]);
        if(gbest<gg[i])gbest=gg[i];
    }
}
void start_search(){
    inc(i,1,m){
        inc(j,1,n){
            r1=xiaorand(0,1,1000);r2=xiaorand(0,1,1000);
            gg[i].v[j]=w*gg[i].v[j]+c1*r1*(gg[i].pbest[j]-gg[i].x[j])+c2*r2*(gbest.x[j]-gg[i].x[j]);
            if(gg[i].v[j]<-0.05)gg[i].v[j]=-0.05;
            else if(gg[i].v[j]>0.05)gg[i].v[j]=0.05;
            gg[i].x[j]+=gg[i].v[j];
            if(gg[i].x[j]<0.0)gg[i].x[j]=0.0;
            else if(gg[i].x[j]>1.0)gg[i].x[j]=1.0;
        }
        int val=f(gg[i]);
        if(val>gg[i].value){
            gg[i].value=val;
            inc(j,1,n)gg[i].pbest[j]=gg[i].x[j];
        }
        if(gbest<gg[i])gbest=gg[i];
    }
}
int main()
{
    srand(time(NULL));
    n=read();m=read();W=read();
    inc(i,1,n)a[i].w=read(),a[i].val=read();
    w=1,c1=0.05,c2=0.05;
    init_gg();
    inc(i,1,1000){
        start_search();
        printf("%4d : %d\n",i,gbest.value);
    }
    return 0;
}
/*
12 50 25
2 5
2 5
4 8
10 2
2 3
5 7
3 6
8 6
4 7
6 6
3 3
2 1
*/