算法设计与分析--回溯算法

 

01:素数环-子集树和排列树分别实现(重点)

总时间限制: 

3000ms

内存限制: 

65535kB

描述

素数环:从1到n这n个数摆成一个环,要求相邻的两个数的和是一个素数。

输出可能的方案数

输入

一个整数n(n<=20)

输出

可能的方案个数

样例输入

2

样例输出

2

子集数

#include <iostream>
#include <cstdlib>
#include <cmath>
using namespace std;

bool b[21]={0};//标志i是否出现在素数环中
int total=0,a[21]={0};//a记录素数环中的每一个数
void search(int t);//回溯过程。形参表示素数环中的数的编号
void print();//输出方案
bool pd(int,int);//判断两数之和是不是素数

void search(int t)//寻找所有解
{
    int i;
    for(i=1;i<=20;i++)
    {
        if(pd(a[t-1],i)&&(!b[i])) //与前一个数是否构成素数及该数是否可用

        {
            a[t]=i;//i进入素数环
            b[i]=1;//标志1
            if(t==20)//到最后得到一个解
            {
                if(pd(a[20],a[1]))//判断最后一个和第一个
                    print();
            }
            else
                search(t+1);
            b[i]=0;
        }
    }
}

void print()
{
    total++;
    cout<<"<"<<total<<">";
    for(int j=1;j<=20;j++)
        cout<<a[j]<<" ";
    cout<<endl;
}

bool pd(int x,int y)
{
    int k=2,i=x+y;
    while(k<=sqrt(i)&&i%k!=0)
        k++;
    if(k>sqrt(i))
        return 1;
    else
        return 0;
}

int main()
{
    for(int i=1;i<=20;i++)
        a[i]=i;
    search(1);
    cout<<total<<endl; //输出总方案数
    return 0;
}

排列树

#include <iostream>
#include <cstdlib>
#include <cmath>
using namespace std;

bool b[21]={0};//标志i是否出现在素数环中
int total=0,a[21]={0};//a记录素数环中的每一个数
void search(int t);//回溯过程。形参表示素数环中的数的编号
void print();//输出方案
bool pd(int,int);//判断两数之和是不是素数

void search(int t)//寻找所有解
{
    int i;
    for(i=t;i<=20;i++)
    {
        swap(a[t],a[i]);
        if(pd(a[t-1],a[t]))
        {
            if(t==20)
            {
                if(pd(a[20],a[1]))
                    print();
            }
            else
                search(t+1);
        }
        swap(a[t],a[i]);
    }
}

void print()
{
    total++;
    cout<<"<"<<total<<">";
    for(int j=1;j<=20;j++)
        cout<<a[j]<<" ";
    cout<<endl;
}

bool pd(int x,int y)
{
    int k=2,i=x+y;
    while(k<=sqrt(i)&&i%k!=0)
        k++;
    if(k>sqrt(i))
        return 1;
    else
        return 0;
}

int main()
{
    for(int i=1;i<=20;i++)
        a[i]=i;
    search(1);
    cout<<total<<endl; //输出总方案数
    return 0;
}

02:装载问题--最大装载量(重点)

分析: 由于集装箱问题是从n个集装箱里选择一部分集装箱,假设解向量为X(x1, x2, …, xn),其中xi∈{0, 1}, xi =1表示集装箱i装上轮船, xi =0表示集装箱i不装上轮船。

输入:每组测试数据:第1行有2个整数c和n。c是轮船的载重量(0<c<30000),n是集装箱的个数(n≤20)。第2行有n个整数w1, w2, …, wn,分别表示n个集装箱的重量。

输出:对每个测试例,输出两行:第1行是装载到轮船的最大载重量,第2行是集装箱的编号。

#include <iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;

class goods{
    int weight;//集装箱重量
public:
    goods(int w=0):weight(w)
    {}
    int get_w()
    {
        return weight;
    }
    void set(int w)
    {
        weight=w;
    }
};

void load(goods *g, int *x, int t, int n,int cw, int &bestcw ,int *best,int r,int c)//轮船载重

//goods *g集装箱列表,*x满足当前最大荷载的装载方案,t子集树的层号,n集装箱的总数

//cw当前的轮船的荷载,bestcw当前的最大荷载,*best待求解的最优装载方案;r剩余的未考察的集装箱重量和
{
    if(t==n)//已经遍历的到叶子结点
    {
        if(cw>bestcw)//当前的轮船的荷载>当前的最大荷载
        {
            for(int i=0;i<n;i++)
            {
                best[i]=x[i];
            }
             bestcw=cw;
        }
    }
    else
        {
            r=r-g[t].get_w();//剩余未处理的物品的重量和,与是否选取当前物品无关
            if(cw+g[t].get_w()<=c)//约束条件进行剪枝
            {
                x[t]=1;
                cw=cw+g[t].get_w();//当前装入的物品的重量和
                load(g,x,t+1,n,cw,bestcw,best,r,c);
                cw=cw-g[t].get_w(); //回溯的需要
            }
            if(cw+r>bestcw)//限界规则
            {
                x[t]=0;
                load(g,x,t+1,n,cw,bestcw,best,r,c);
            }
            r=r+g[t].get_w(); //回溯的需要
        }
}

int main()
{
    int n,c,bestcw=0;
    int *x,*best, r=0;
    cout<< "请输入轮船的装载重量c和物品的件数n:";
    cin>>c>>n;
    goods *g;
    g=new goods[n];
    x=new int [n];
    best=new int[n];
    cout<<"请输入每件物品的重量:";
    for(int i=0;i<n;i++)
    {
        int w;
        cin>>w;
        g[i].set(w);
        r=r+w;
    }
    load(g,x,0,n,0,bestcw,best,r,c);
    cout<<bestcw<<endl;
    for(int i=0;i<n;i++)
        cout<<best[i]<<"  ";
    cout<<endl;
    return 0;
}

03:0-1背包问题(重点)

输入: 第一个数据是背包的容量为c(1≤c≤1500),第二个数据是物品的数量为n(1≤n≤50)。接下来n行是物品i的重量是wi,其价值为vi。所有的数据全部为整数,且保证输入数据中物品的总重量大于背包的容量。当c=0时,表示输入数据结束。

输出: 对每组测试数据,输出装入背包中物品的最大价值。

分析:0-1背包问题(子集数)

令cw(i)表示目前搜索到第i层已经装入背包的物品总重量
令cv(i)表示目前到第i层结点已经装入背包的物品价值
对于左子树, xi =1 ,其约束函数为:construction(i)=cw(i-1)+wi。若constraint(i)>W,则停止搜索左子树,否则继续搜索。
对于右子树,r(i)表示剩余物品的总价值,限界函数Bound(i)=cv(i)+r(i)
r(i)越小, Bound(i)越小,剪掉的分支就越多(越靠近树根剪枝,剪掉的分支越多。从而能加快搜索速度)。
为了构造更小的r(i) ,将物品以单位重量价值比di=vi/wi递减的顺序进行排列(贪心策略)

#include <bits/stdc++.h>
using namespace std;

int c,n;//背包容量、物品数量
int cw,cv;//背包当前重量,当前价值
int bestv;//当前最优价值

struct bag
{
    int w; //物品的重量
    int v; //物品的价值
    double d; //物品的单位重量价值比
}a[101];

bool cmp(bag a,bag b)//按照价值比从大到小排序
{
    if(a.d>=b.d)
        return true;
    else
        return false;
}

int Bound(int i)//限界函数
{
    int cleft=c-cw;//背包剩余的容量
    int b=cv;//上界
    //尽量装满背包
    while(i<n&&a[i].w<=cleft)
    {
        cleft-=a[i].w;
        b+=a[i].v;
        i++;
    }
    /*剩余的部分空间也装满。0-1是可能装不满的。
    但此处主要计算最大值,所以,需要从装满的角度考虑该问题
    (尽管这不是一个可行解,但可以证明其价值是最优值的上界)
    */
    if(i<n)
        b+=1.0*cleft*a[i].v/a[i].w;
    return b;
}

int ff(int i)//形参i是回溯的深度,从0开始.商品编号从0开始编号
{
    if(i==n)
    {
         bestv=cv;//当前最优价值=当前价值
         //return ;
    }
    //进入左子树搜索, 表示选择第i件物品
    if(cw+a[i].w<=c)//满足约束条件
    {
        cw+=a[i].w;//选择第i件物品
        cv+=a[i].v;
        ff(i+1);
        cw-=a[i].w;//回溯的需要
        cv-=a[i].v;//回溯需要
    }
    //进入右子树搜索,表示不选择第i件物品,相关的cw, cv不改变
    if(Bound(i+1)>bestv)//大于当前最优价值,满足条件不剪枝
    {
        ff(i+1);
    }
    return 0;
}

int main()
{
    cin>>c>>n;
    for(int i=0;i<n;i++)
    {
        cin>>a[i].w>>a[i].v;
        a[i].d=1.0*a[i].v/a[i].w;
    }
    sort(a,a+n,cmp);
    ff(0);
    cout<<bestv<<endl;
}

04:n皇后问题(重点)

 

 

#include<bits/stdc++.h>
using namespace std;

int n;//棋盘大小
int x[101];//解向量
int sum=0;//当前已经找到的可行方案数

inline bool Place(int t)//判断是否在同一行,同一列,同一对角线
{
    int i;
    for(i=1;i<t;i++)
    {
        if ((abs(t-i) == abs(x[i]-x[t])) || (x[i] == x[t]))
            //判断是否在同一对角线、同一行
            return false;
    }
    return true;
}

void Back(int t)//t是回溯的深度,从1开始
{
    int i;
    if(t>n)//到达叶子结点
    {
        sum++;//方案数+1
        for(i=1;i<=n;i++)
            cout<<x[i]<<" ";
        cout<<endl;
    }
    else
    {
        for(i=1;i<=n;i++)
        {
            x[t]=i;
            if(Place(t))//判断是否符合条件
                Back(t+1);
        }
    }
}

int main()
{
    cin>>n;
    Back(1);
    cout<<sum<<endl;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值