【人工智能实验】蚁群算法解决TSP问题_蚁群算法解决tsp问题是人工智能的-程序员宅基地

技术标签: 算法  C++  c++  人工智能  

定义了Point类忘了用了结果代码变成了粪山
其中涉及贪婪算法,轮盘赌法,以及看不懂自己写的什么了

#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <deque>
#include <iostream>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>
#include <time.h>
#include<fstream>
#include<sstream>
using namespace std;

//created by jzyghsh

#define MAXPATH 9999999

struct Point {
     //city class
    string name;
    double x, y;
    int num;
};


void init();
double distance(Point a, Point b);
double greedyPath(vector<Point> v);
bool isIn(string str, vector<string> list);
void nextCity();
int findNumOfCity(string name);
string RWS(vector<string> visitList);
void update();
int n;//city number
int m;//ant number
int alpha = 1, beta = 2;
double rou = 0.5;
vector<Point> cities;//city group (unchaged)
vector<vector<double>> dist;//distance list (unchaged)
vector<vector<double>> message;//message level on every road
vector<vector<string>> visitedCities;//abandon list 
vector<double> p;
vector<double> C;
double minPath = FLT_MAX;
vector<string> bestPath;
int nextCityNum;

int main() {
    
#ifndef ONLINE_JUDGE
    freopen("E:\\txts\\wi29.txt", "r", stdin);
#endif
    int times = 0;
    int i;
    double difference = 0.0;
    vector<int> count;
    init();
    while (times<=100) {
    
        nextCity();
        update();
        times++;
    }
    nextCity();
    cout << "共生成" << times * m << "只蚂蚁" << endl;
    
    cout << "最终路径为:";
    for (i = 0; i < bestPath.size(); i++) {
    
        cout << bestPath[i] << "->";
    }
    cout << "最终路径长度为:" << minPath << endl;
    return 0;
}

void nextCity() {
    
    int i, j, k, begin, first;
    string nextVisit;
    int notVisit,nowVisit;
    double up, pi, totalUp = 0, ci = 0;
    C.clear();
    for (i = 0; i < m; i++) {
    
        notVisit = n - visitedCities[i].size();
        ci = 0;
        C.push_back(ci);
        vector<string> visitList;
        first = 0;
        while (notVisit > 0) {
    
            visitList.clear();//available visit list
            totalUp = 0;
            p.clear();
            for (j = 0; j < n; j++) {
    
                if (visitedCities[i].back() == cities[j].name) {
    
                    nowVisit = j;
                    if (first == 0) {
    
                        begin = j;
                        first = 1;
                    }
                    break;
                }
            }
            for (j = 0; j < n; j++) {
    
                if (!isIn(cities[j].name, visitedCities[i])) {
    
                    //cout << cities[j].name << "not in" << endl;
                    up = pow(message[nowVisit][j], alpha) * pow((1 / dist[nowVisit][j]), beta);//i有问题
                    totalUp += up;
                    p.push_back(up);
                    visitList.push_back(cities[j].name);
                    //cout << visitList.back() << " ";
                }
            }
            //cout << endl;
            //created by jzyghsh
            //cout << p.size() << endl;
            for (j = 0; j < p.size(); j++) {
    
                p[j] /= totalUp;
            }
            nextVisit = RWS(visitList);
            visitedCities[i].push_back(nextVisit);
            //cout << "add" << visitList[nextVisit] << endl;
            for (k = 0; k < n; k++) {
    
                if (nextVisit == cities[k].name) {
    
                    nextCityNum = k;
                }
            }
            C[i] += dist[nowVisit][nextCityNum];
            //cout << cities[nowVisit].name << " " << cities[nextCityNum].name << endl;
            notVisit--;
        }
        //cout << endl;
        C[i] += dist[nextCityNum][begin];
        //cout << cities[nextCityNum].name << " " << cities[begin].name << endl;
    }
    for (i = 0; i < m; i++) {
    
        cout << "蚂蚁" << i << "爬过的路径为:";
        visitedCities[i].push_back(visitedCities[i][0]);
        for (j = 0; j < visitedCities[i].size(); j++) {
    
            cout << visitedCities[i][j] << "->";
        }
        cout << "路径长度为:" << C[i] << endl;
        if (C[i] < minPath) {
    
            minPath = C[i];
            bestPath = visitedCities[i];
        }
    }
}

string RWS(vector<string> visitList) {
    //select a city
    double q, area = 0, temp;
    string t;
    int i, j;
    q = rand()/double(RAND_MAX);//0~1 double number
    //cout << "q=" << q;
    for (i = 0; i < visitList.size()-1; i++) {
    
        for (j = i; j < visitList.size() - 1; j++) {
    
            if (p[i] > p[i + 1]) {
    
                temp = p[i];
                p[i] = p[i + 1];
                p[i + 1] = temp;
                t = visitList[i];
                visitList[i] = visitList[i + 1];
                visitList[i + 1] = t;
            }
        }
    }
    for (i = 0; i < p.size(); i++) {
    
        area += p[i];
        if (q < area) {
    
            //cout << "choose" << i << endl;
            nextCityNum = i;
            return visitList[i];
        }
    }
    return visitList[0];
}

void update() {
    
    int i, j;
    int cityA = 0, cityB = 0;
    for (i = 0; i < n; i++) {
    // message fade
        for (j = 0; j < n; j++) {
    
            message[i][j] *= (1 - rou);
        }
    }
    for (i = 0; i < m; i++) {
    //message adding
        for (j = 0; j < visitedCities[i].size()-1; j++) {
    
            cityA = findNumOfCity(visitedCities[i][j]);
            cityB = findNumOfCity(visitedCities[i][j + 1]);
            message[cityA][cityB] += 1 / C[i];
            //cout << cities[cityA].name << "to" << cities[cityB].name << endl;
        }
        //cout << endl;
    }
    for (i = 0; i < m; i++) {
    //randomly put ants again
        int ran;
        visitedCities[i].clear();
        ran = rand() % n;
        visitedCities[i].push_back(cities[ran].name);
        cout << "蚂蚁" << i << "初始在城市" << visitedCities[i].back() << endl;
    }
}

int findNumOfCity(string name) {
    
    int i;
    for (i = 0; i < cities.size(); i++) {
    
        if (name == cities[i].name) {
    
            return i;
        }
    }
    return NULL;
}
void init() {
    
    int i = 0, j = 0;
    double Cnn,t0;
    cin >> n >> m;
    for (i = 0; i < n; i++) {
    
        Point p;
        cin >> p.name >> p.x >> p.y;
        p.num = i;
        cities.push_back(p);
    }
    for (i = 0; i < n; i++) {
    //init distance list
        vector<double> v;
        dist.push_back(v);
        for (j = 0; j < n; j++) {
    
            dist[i].push_back(distance(cities[i], cities[j]));
        }
    }
    Cnn = greedyPath(cities);
    t0 = m / Cnn;
    cout << "用贪婪算法得到的初始路径长度为:" << Cnn << endl;
    for (i = 0; i < n; i++) {
    //init message list
        vector<double> v;
        message.push_back(v);
        for (j = 0; j < n; j++) {
    
            if (i != j) {
    
                message[i].push_back(t0);
            }
            else {
    
                message[i].push_back(0);
            }
        }
    }
    for (i = 0; i < m; i++) {
    //randomly put ants
        int ran;
        vector<string> str;
        ran = rand() % n;
        str.push_back(cities[ran].name);
        visitedCities.push_back(str);
        str.clear();
        cout << "蚂蚁" << i << "初始在城市" << visitedCities[i].back() << endl;
    }
    
}

double distance(Point a,Point b) {
    
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}

bool isIn(string str,vector<string> list) {
    
    int i = 0, length = list.size();
    for (i = 0; i < length; i++) {
    
        if (str == list[i]) {
    
            return true;
        }
    }
    return false;
}

double greedyPath(vector<Point> v) {
    //greedy method get Cnn
    double res = 0, nowPath, minPath;
    int count = 0, begin, end;
    int i = 0, j = 0;
    vector<string> checked;
    string choseCity;
    int choseCityNum;
    srand((unsigned)time(NULL));
    begin = rand() % n;
    end = begin;
    checked.push_back(v[begin].name);
    choseCityNum = begin;
    while (count < v.size()-1) {
    
        minPath = MAXPATH;
        for (i = 0; i < n; i++) {
    
            if (i != begin) {
    
                nowPath = distance(v[begin], v[i]);
                if (nowPath < minPath && !isIn(v[i].name,checked)) {
    
                    minPath = nowPath;
                    choseCity = v[i].name;
                    choseCityNum = i;
                }
            }
        }
        checked.push_back(choseCity);
        begin = choseCityNum;
        res += minPath;
        count++;
    }
    checked.push_back(checked[0]);
    cout << "贪婪算法路径为:";
    for (i = 0; i < checked.size(); i++) {
    
        cout << checked[i] << "->";
    }
    cout << endl;
    res += dist[begin][end];
    //cout << begin << " " << end;
    return res;
}

在这里插入图片描述
在这里插入图片描述

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_48609829/article/details/121187770

智能推荐

从零开始搭建Hadoop_创建一个hadoop项目-程序员宅基地

文章浏览阅读331次。第一部分:准备工作1 安装虚拟机2 安装centos73 安装JDK以上三步是准备工作,至此已经完成一台已安装JDK的主机第二部分:准备3台虚拟机以下所有工作最好都在root权限下操作1 克隆上面已经有一台虚拟机了,现在对master进行克隆,克隆出另外2台子机;1.1 进行克隆21.2 下一步1.3 下一步1.4 下一步1.5 根据子机需要,命名和安装路径1.6 ..._创建一个hadoop项目

心脏滴血漏洞HeartBleed CVE-2014-0160深入代码层面的分析_heartbleed代码分析-程序员宅基地

文章浏览阅读1.7k次。心脏滴血漏洞HeartBleed CVE-2014-0160 是由heartbeat功能引入的,本文从深入码层面的分析该漏洞产生的原因_heartbleed代码分析

java读取ofd文档内容_ofd电子文档内容分析工具(分析文档、签章和证书)-程序员宅基地

文章浏览阅读1.4k次。前言ofd是国家文档标准,其对标的文档格式是pdf。ofd文档是容器格式文件,ofd其实就是压缩包。将ofd文件后缀改为.zip,解压后可看到文件包含的内容。ofd文件分析工具下载:点我下载。ofd文件解压后,可以看到如下内容: 对于xml文件,可以用文本工具查看。但是对于印章文件(Seal.esl)、签名文件(SignedValue.dat)就无法查看其内容了。本人开发一款ofd内容查看器,..._signedvalue.dat

基于FPGA的数据采集系统(一)_基于fpga的信息采集-程序员宅基地

文章浏览阅读1.8w次,点赞29次,收藏313次。整体系统设计本设计主要是对ADC和DAC的使用,主要实现功能流程为:首先通过串口向FPGA发送控制信号,控制DAC芯片tlv5618进行DA装换,转换的数据存在ROM中,转换开始时读取ROM中数据进行读取转换。其次用按键控制adc128s052进行模数转换100次,模数转换数据存储到FIFO中,再从FIFO中读取数据通过串口输出显示在pc上。其整体系统框图如下:图1:FPGA数据采集系统框图从图中可以看出,该系统主要包括9个模块:串口接收模块、按键消抖模块、按键控制模块、ROM模块、D.._基于fpga的信息采集

微服务 spring cloud zuul com.netflix.zuul.exception.ZuulException GENERAL-程序员宅基地

文章浏览阅读2.5w次。1.背景错误信息:-- [http-nio-9904-exec-5] o.s.c.n.z.filters.post.SendErrorFilter : Error during filteringcom.netflix.zuul.exception.ZuulException: Forwarding error at org.springframework.cloud..._com.netflix.zuul.exception.zuulexception

邻接矩阵-建立图-程序员宅基地

文章浏览阅读358次。1.介绍图的相关概念  图是由顶点的有穷非空集和一个描述顶点之间关系-边(或者弧)的集合组成。通常,图中的数据元素被称为顶点,顶点间的关系用边表示,图通常用字母G表示,图的顶点通常用字母V表示,所以图可以定义为:  G=(V,E)其中,V(G)是图中顶点的有穷非空集合,E(G)是V(G)中顶点的边的有穷集合1.1 无向图:图中任意两个顶点构成的边是没有方向的1.2 有向图:图中..._给定一个邻接矩阵未必能够造出一个图

随便推点

MDT2012部署系列之11 WDS安装与配置-程序员宅基地

文章浏览阅读321次。(十二)、WDS服务器安装通过前面的测试我们会发现,每次安装的时候需要加域光盘映像,这是一个比较麻烦的事情,试想一个上万个的公司,你天天带着一个光盘与光驱去给别人装系统,这将是一个多么痛苦的事情啊,有什么方法可以解决这个问题了?答案是肯定的,下面我们就来简单说一下。WDS服务器,它是Windows自带的一个免费的基于系统本身角色的一个功能,它主要提供一种简单、安全的通过网络快速、远程将Window..._doc server2012上通过wds+mdt无人值守部署win11系统.doc

python--xlrd/xlwt/xlutils_xlutils模块可以读xlsx吗-程序员宅基地

文章浏览阅读219次。python–xlrd/xlwt/xlutilsxlrd只能读取,不能改,支持 xlsx和xls 格式xlwt只能改,不能读xlwt只能保存为.xls格式xlutils能将xlrd.Book转为xlwt.Workbook,从而得以在现有xls的基础上修改数据,并创建一个新的xls,实现修改xlrd打开文件import xlrdexcel=xlrd.open_workbook('E:/test.xlsx') 返回值为xlrd.book.Book对象,不能修改获取sheett_xlutils模块可以读xlsx吗

关于新版本selenium定位元素报错:‘WebDriver‘ object has no attribute ‘find_element_by_id‘等问题_unresolved attribute reference 'find_element_by_id-程序员宅基地

文章浏览阅读8.2w次,点赞267次,收藏656次。运行Selenium出现'WebDriver' object has no attribute 'find_element_by_id'或AttributeError: 'WebDriver' object has no attribute 'find_element_by_xpath'等定位元素代码错误,是因为selenium更新到了新的版本,以前的一些语法经过改动。..............._unresolved attribute reference 'find_element_by_id' for class 'webdriver

DOM对象转换成jQuery对象转换与子页面获取父页面DOM对象-程序员宅基地

文章浏览阅读198次。一:模态窗口//父页面JSwindow.showModalDialog(ifrmehref, window, 'dialogWidth:550px;dialogHeight:150px;help:no;resizable:no;status:no');//子页面获取父页面DOM对象//window.showModalDialog的DOM对象var v=parentWin..._jquery获取父window下的dom对象

什么是算法?-程序员宅基地

文章浏览阅读1.7w次,点赞15次,收藏129次。算法(algorithm)是解决一系列问题的清晰指令,也就是,能对一定规范的输入,在有限的时间内获得所要求的输出。 简单来说,算法就是解决一个问题的具体方法和步骤。算法是程序的灵 魂。二、算法的特征1.可行性 算法中执行的任何计算步骤都可以分解为基本可执行的操作步,即每个计算步都可以在有限时间里完成(也称之为有效性) 算法的每一步都要有确切的意义,不能有二义性。例如“增加x的值”,并没有说增加多少,计算机就无法执行明确的运算。 _算法

【网络安全】网络安全的标准和规范_网络安全标准规范-程序员宅基地

文章浏览阅读1.5k次,点赞18次,收藏26次。网络安全的标准和规范是网络安全领域的重要组成部分。它们为网络安全提供了技术依据,规定了网络安全的技术要求和操作方式,帮助我们构建安全的网络环境。下面,我们将详细介绍一些主要的网络安全标准和规范,以及它们在实际操作中的应用。_网络安全标准规范