【先进的算法】Lasvegas算法3SAT问题(C++实现代码)

转载请注明出处:http://blog.csdn.net/zhoubin1992/article/details/46469557


1.SAT问题描写叙述

命题逻辑中合取范式 (CNF) 的可满足性问题 (SAT)是当代理论计算机科学的核心问题, 是一典型的NP 全然问题.在定义可满足性问题SAT之前。先引进一些逻辑符号。

一个 SAT 问题是指: 对于给定的 CNF 是否存在一组关于命题变元的真值指派使得A 为真. 显然, 假设A 为真, 则 CNF 的每一个子句中必有一个命题变元为 1 (真) 。


2.Las Vegas算法

Las Vegas 算法是利用随机值做出随机选择的一种概率算法,而且不会产生不对的答案。在计算过程中所做出的随机选择,可能使算法比其它算法更快地得到所要求的解。

拉斯维加斯算法不会得到不对的解。一旦用拉斯维加斯算法找到一个解,这个解就一定是正确解。但有时用拉斯维加斯算法找不到解。与蒙特卡罗算法相似,拉斯维加斯算法找到正确解的概率随着它所用的计算时间的添加而提高。

对于所求解问题的任一实例,用同一拉斯维加斯算法反复对该实例求解足够多次,可使求解失败的概率随意小。

Las Vegas 算法用来搜索包括目标结点的解空间。它用一些随机选择来移动,而不须要在每一个结点都计算一个新的结点。假设成功结点的比例在解空间中相当高,则找到目标结点的概率可能非常高。

当下一个结点的计算比較困难或者系统化地搜索没有什么必要时。採用Las Vegas 算法,会提高计算的效率。当然,下一个结点的随机选择有可能导致找不到成功的结点,可是我们能够反复多次运行,来提高目标结点的效率。

拉斯维加斯算法的一个显著特征是它所作的随机性决策有可能导致算法找不到所需的解,可是通过反复多次运行来克服,在求解NP难问题时,用它往往会收到奇效。


3.C++实现代码

// lasvegas3SAT.cpp : 定义控制台应用程序的入口点。
//
/*********************************
-----------------------------------
Lasvegas算法解决3SAT问题(C++实现代码)
-----------------------------------
Author:牧之丶  Date:2014年
Email:[email protected]
**********************************/
#include "stdafx.h"
#include <stdlib.h>
#include <string.h>
#include <time.h>
#include <iostream>

const int n=250;
int M[n][3];
int sign[3*n+1];
int x[101],y[101];
int ok[2];
bool Place( int k)
{
    //memset(y,1,101);
    int t;
    /*for( int j = 1; j <= k - 1; j++)
    {
        y[j]=~x[j];
    }*/
    for(int i = 0; i < n; i++)
    {
        t=0;
        for(int j = 0; j < 3; j++)
        {
            if(M[i][j]>k)
                t = t+1;
            else
                t = t+!(sign[i*3+j]^x[M[i][j]] );
        }
        if(t<1)
            return false;
    }
    return true;
}
bool SAT_True(int x[])
{
    int k = 1;
    int count = 0;
    int i;
    while( k <= 100 )
    {
        count = 0;
        for( i = 0; i <= 1 ; i++ )
        {
            x[k] = i;
            if( Place(k))
            {
                ok[count] = i;
                count ++;
            }
        }
        if( count == 0 ) return false;
        i = ok[rand() % count];
        x[k] = i;
        k++;
    }
    return true;
}

int _tmain(int argc, _TCHAR* argv[])
{       //反复运行20次
    //int ncase = 20;
    double run_time = 0.0; //运行时间
    double run_num = 0.0; //运行次数
    time_t start,end;
    srand(time(0));
    //while(ncase--)
    //{
        start = clock();
        for(int i=0;i<n;i++)
            for(int j=0;j<3;j++)
                M[i][j] = rand()%100+1;
        for(int i=1;i<=3*n;i++)
            sign[i] = rand()%2;

        memset(x,0,101*sizeof(int));
        int k=1;
         while(!SAT_True(x) )
        {
            k++;
            if(k > 100000)
            {
                printf("failed!\n");
                break;
            }
        }
        end = clock();
        run_num = k;
        run_time += (end - start)/CLOCKS_PER_SEC;
        if(k <= 100000)
            std::cout << "运行了" << run_num << "次" << std::endl;
    //}

    printf("the running time is : %f\n", run_time);
    system("pause");
    return 0;
}


4.实验结果及分析

为了測试Las Vegas 的计算效果, 我们用随机产生的3-SAT 模型(每一个子句的长度 l= 3, 且子句里的变元两两不同) 做实例。每种取值运行20次。考虑有可能找不到解的情况,当搜索次数超过十万次,觉得此例子不可满足。

得到的结果为:



參考文献

[1] 张德富.算法设计与分析(高级教程)[M].国防工业出版社,2007.

时间: 2024-10-18 13:47:45

【先进的算法】Lasvegas算法3SAT问题(C++实现代码)的相关文章

模拟退火算法解决3SAT问题(C++实现代码)

转载请注明出处:http://blog.csdn.net/zhoubin1992/article/details/46453761 1 SAT问题描述 命题逻辑中合取范式 (CNF)的可满足性问题 (SAT)是当代理论计算机科学的核心问题,是一典型的NP完全问题.在定义可满足性问题SAT之前,先引进一些逻辑符号. 2  模拟退火算法 模拟退火算法来源于固体退火原理,将固体加温至充分高,再让其徐徐冷却,加温时,固体内部粒子随温升变为无序状,内能增大,而徐徐冷却时粒子渐趋有序,在每个温度都达到平衡态

【啊哈!算法】算法7:Dijkstra最短路算法

上周我们介绍了神奇的只有五行的Floyd最短路算法,它可以方便的求得任意两点的最短路径,这称为“多源最短路”.本周来来介绍指定一个点(源点)到其余各个顶点的最短路径,也叫做“单源最短路径”.例如求下图中的1号顶点到2.3.4.5.6号顶点的最短路径. <ignore_js_op> 与Floyd-Warshall算法一样这里仍然使用二维数组e来存储顶点之间边的关系,初始值如下. <ignore_js_op> 我们还需要用一个一维数组dis来存储1号顶点到其余各个顶点的初始路程,如下.

【啊哈!算法】算法9:开启树之旅

这是什么?是一个图?不对,确切的说这是一棵树.这哪里像树呢?不要着急我们来变换一下. 是不是很像一棵倒挂的树,也就是说它是根朝上,而叶子朝下的.不像?哈哈,看完下面这幅图你就会觉得像啦. 你可能会问:树和图有什么区别?这个称之为树的东西貌似和无向图差不多嘛.不要着急,继续往下看.树其实就是不包含回路的连通无向图.你可能还是无法理解这其中的差异,举个例子,如下.          上面这个例子中左边的是一棵树,而右边的是一个图.因为左边的没有回路,而右边的存在1->2->5->3->

【啊哈!算法】算法10:二叉树

二叉树是一种特殊的树.二叉树的特点是每个结点最多有两个儿子,左边的叫做左儿子,右边的叫做右儿子,或者说每个结点最多有两棵子树.更加严格的递归定义是:二叉树要么为空,要么由根结点.左子树和右子树组成,而左子树和右子树分别是一棵二叉树. 下面这棵树就是一棵二叉树. 二叉树的使用范围最广,一棵多叉树也可以转化为二叉树,因此我们将着重讲解二叉树. 二叉树中还有连两种特殊的二叉树叫做满二叉树和完全二叉树.如果二叉树中每个内部结点都有两个儿子,这样的二叉树叫做满二叉树.或者说满二叉树所有的叶结点都有同样的深

九章算法 基础算法 强化算法 系统设计 大数据 安卓 leetcode 高清视频

leetcode 直播视频讲座录像 九章算法视频录像,PPT 算法班,算法强化班,Java入门与基础算法班,big data项目实战班,Andriod项目实战班 九章算法下载 九章算法面试 九章算法leetcode 九章算法答案 九章算法mitbbs 九章算法班 九章算法ppt 九章算法录像 九章算法培训 九章算法微博 leetcode 视频 九章算法偷录 算法培训 算法班课程大纲: 1 从strStr谈面试技巧与Coding Style(免费试听) 2 二分搜索与旋转排序数组 Binary S

编程算法 - 把数组排成最小的数 代码(C)

把数组排成最小的数 代码(C) 本文地址: http://blog.csdn.net/caroline_wendy 题目: 输入一个正整数数组, 把数组里所有数字拼接起来排成一个数, 打印能拼接出的所有数字中最小的一个. 大数转换为字符串, 重载快速排序的比较方法, 进行排序, 最后拼接. 代码: /* * main.cpp * * Created on: 2014.6.12 * Author: Spike */ /*eclipse cdt, gcc 4.8.1*/ #include <stdi

编程算法 - 最小的k个数 红黑树 代码(C++)

最小的k个数 红黑树 代码(C++) 本文地址: http://blog.csdn.net/caroline_wendy 题目: 输入n个整数, 找出其中的最小k个数. 使用红黑树(multiset), 每次替换最大的值, 依次迭代. 时间复杂度: O(nlogk). 代码: /* * main.cpp * * Created on: 2014年6月29日 * Author: wang */ #include <iostream> #include <vector> #includ

[算法]有趣算法合辑[31-40]

题目31:八进制转换为十进制 1.程序分析: 2.程序源代码: main() { char *p,s[6];int n; p=s; gets(p); n=0; while(*(p)!='\0') {n=n*8+*p-'0'; p++;} printf("%d",n); } 题目32:求0-7所能组成的奇数个数. 1.程序分析: 2.程序源代码: main() { long sum=4,s=4; int j; for(j=2;j<=8;j++)/*j is place of num

【坐在马桶上看算法】算法9:开启“树”之旅

我们先来看一个例子. 这是什么?是一个图?不对,确切的说这是一棵树.这哪里像树呢?不要着急我们来变换一下. 是不是很像一棵倒挂的树,也就是说它是根朝上,而叶子朝下的.不像?哈哈,看完下面这幅图你就会觉得像啦. 你可能会问:树和图有什么区别?这个称之为树的东西貌似和无向图差不多嘛.不要着急,继续往下看.树其实就是不包含回路的连通无向图.你可能还是无法理解这其中的差异,举个例子,如下.          上面这个例子中左边的是一棵树,而右边的是一个图.因为左边的没有回路,而右边的存在1->2->5