OpenMP并行加速笛卡尔乘积
1.字典字符集的笛卡尔乘积
问题描述:
对于给定的由字典字符集组合而成的表达式,求该表达式构成的所有元素。例如表达式[0-9][a-z]
,其中0-9表示10个数字,a-z表示26个小写字母,构成的所有元素就是0a,0b,…,0z,1a,1b,…9z。字典字符集的笛卡尔乘积示意如下:
问题分析:
对于任意的一个字典字符集构成的表达式[dic0][dic1]...[dicn]
,从左至右可以看作按照高位到低位的一个由字典元素组成的一个“数”,这样比较符合我们日常表示数值的高低位的习惯。比如如果字典都是[0-9],那么表达式[0-9][0-9]
表示的就是一个数值字符串00~99。笛卡尔乘积的空间是各个字典高度的乘积,给定其空间中的任意一个元素下标,就可以对应到每个字典中的元素下标。比如[0-9[0-9]
的笛卡尔乘积的空间是各个字典高度的乘积10*10=100,空间中第0个元素就是00,第99个元素就是99。
每一个字典元素都有一个位权重。表达式[0-9[0-9]
从左至右第一个字典的位权重w=10,第二个字典的位权重w=1。我们平常所说的个位、十位、百位,就是根据数值位的位权重来称呼它。位权重的意义在于,数值是其位权重的多少倍,那么它就取第几个元素。比如第99个元素(下标从0开始),那么数值99是十位的位权重w=10的9倍,所以元素为字符‘9’,对数值99取w=10的余数得9,同理9是个位的位权重w=1的9倍,所以元素为字符’9’,所以构成了字符串99。
实现示例:
对表达式[0-9][a-z[A-Z]
,其实现笛卡尔乘积的具体过程可以描述如下:
(1)对从左至右(高位到低位)将各个字典字符集的所在数位的计算单位计算出来,由当前字典右边的字典的高度相乘得到,如[0-9]的计数单位w=26*26=676,[a-z]的计数单位w=26*1=26,[A-Z]的计数单位w=1。
(2)给定笛卡尔乘积空间的元素下标i,根据i找到各个字典内的元素下标的过程如下,从高位开始查找,即从左开始查找。
(2.1)查找字典[0-9]中的元素下标:[0-9].index=i/[0-9].w;
(2.2)查找字典[a-z]中的元素下标[a-z].index:i=i%[0-9].w; [a-z].index=i/[a-z].w;
(2.3)查找字典[A-Z]中的元素下标[A-Z].index:i=i%[a-z].w;[A-Z].index=i/1=i。
(3)将i=0递增至笛卡尔乘积的空间大小减一,即10*26*26-1
,重复步骤2,即可完成表达式[0-9][a-z[A-Z]
的笛卡尔乘积。
例如给定第677个笛卡尔乘积的元素,那么[0-9].index=1,所以取[0-9]内的素‘1’,[a-z].index=671%676/26=0,所以取出元素‘a’,[a-z].index=1/1=1,所以取出元素‘B’,所以第677个元素就是“1aB”。
2.源码
以下代码在Linux平台编译运行,稍作修改可移植到Windows平台。其功能是完成多个字典字符集的笛卡尔乘积。并通过OpenMP并行加速。正确性已在实际项目中通过验证。
#include <pthread.h>
#include <omp.h>#include <iostream>
#include <map>
#include <string>
using namespace std;typedef unsigned char uint8;//字典字符集与段字符集
struct charset_mem{int high,width; //字符集的宽度和高度int mem_size; //字符集data所占用的内存,单位字节uint8 *data; //字符集的数据char name[128]; //字符集名称
};map<string,charset_mem*> dic_utf8_charset_map; //全局字典字符集缓存
map<string,charset_mem*> dic_ucs2_charset_map; //全局字典字符集缓存
map<string,charset_mem*> seg_charset_map; //全局段字符集缓存pthread_mutex_t charset_mutex; //功能:根据多个字典字符集生成相应的笛卡尔乘积
//参数:charsetID:笛卡尔乘积结果字符集名称,dicNum:字典字符集数目,dicName:字典字符集名称数组指针,encode:字典字符编码类型
//返回值:成功返回true,失败返回false
bool cartesianProduct(string charsetID, int dicNum, char(*dicName)[128],uint8 encode)
{pthread_mutex_lock(&charset_mutex); //对字符集的map关联容器修改需要加锁string charsetNewedID=charsetID;map<string,charset_mem*>::iterator iter;charset_mem* segNewedCharset=new charset_mem;memset(segNewedCharset,0,sizeof(charset_mem));strcpy(segNewedCharset->name,charsetNewedID.c_str());#define MAX_WORD_LEN 40#define MAX_DIC_NUM 32//笛卡尔乘积(cartesian product)准备工作map<string,charset_mem*>& dic_charset_map=(0==encode)?dic_utf8_charset_map:dic_ucs2_charset_map;int high=1,width=0;int s[MAX_DIC_NUM]={0}; //字典段进制位for(int i=dicNum-1; i>=0; --i){cout<<"dicName["<<i<<"]:"<<dicName[i]<<endl;iter=dic_charset_map.find(dicName[i]);if(iter==dic_charset_map.end()){cout<<dicName[i]<<"do not exist"<<endl;pthread_mutex_unlock(&charset_mutex);return false;}s[i]=high;high*=iter->second->high;width+=iter->second->width;}segNewedCharset->high=high;segNewedCharset->width=width;segNewedCharset->data=new uint8[high*width]; //笛卡尔乘积int thread_num=omp_get_max_threads();//获取处理器最大可并行的线程数#pragma omp parallel for num_threads(thread_num) for(int i=0;i<high;++i){uint8 wordTmp[MAX_WORD_LEN]={0};map<string,charset_mem*>::iterator iterTmp;int offset=0;int charpos=i; for(int j=0;j<dicNum;++j){iterTmp=dic_charset_map.find(dicName[j]); int indexDic=charpos/s[j];int offsetDic=indexDic*iterTmp->second->width;memcpy(wordTmp+offset,iterTmp->second->data+offsetDic,iterTmp->second->width);charpos=charpos%s[j];offset+=iterTmp->second->width;}memcpy(segNewedCharset->data+i*segNewedCharset->width,wordTmp,segNewedCharset->width);}//将结果字符集添加到,map映射表seg_charset_map.insert(pair<string,charset_mem*>(charsetNewedID,segNewedCharset));pthread_mutex_unlock(&charset_mutex);return true;
}
3.优化
写毕业论文的时候,经实验室的小伙伴提醒,发现其实不用事先求出各个字典所在数位的计数单位,也可以根据给定的笛卡尔乘积的元素下标唯一的找到各个字典中对应的元素。为了避免与论文查重时重复,只贴出图片。具体实现已经抽象为如下算法:
算法中注释中的热词就是上文提到字典,其实现的原理是从表达式的低位到高位计算每一个字典的元素下标,上面未优化的方法是从高位到低位顺序计算。从低位到高位来计算的话,无需事先求出各个字典位的计数单位。因为:当字典位的计数单位的为w=1时,可以通过笛卡尔乘积的元素下标i对其高度h取余,即得到最低字典位字典内的元素下标。当对下一个字典求其元素下标时,需要将下一个字典位的计数单位w’变为1,具体做法就是i除以当前字典的高度向下取整。依次类推,就可以求出各个字典内的元素下标了。具体描述见上面的算法。
以表达式[0-9][a-z[A-Z]
,求笛卡尔乘积中第677个(从0开始)元素的各个字典内的元素下标的过程描述如下:
(1)求字典[A-Z]的元素下标index=i%[A-Z].h=677%26=1,所以去元素’B’;
(2)求字典[a-z]的元素下标index:
(2.1)将[a-z]的计数单位变为1,做法是i=i/[A_Z].h=677/26=26;
(2.2)求[a-z].index=i%[a-z].h=26%26=0。所以取元素’a’。
(3)求字典[0-9]的元素下标index:
(3.1)将[0-9]的计数单位变为1,做法是i=i/[a_z].h=26/26=1;
(3.2)求[a-z].index=i%[0-9].h=1%10=1。所以取元素’1’。
所以第677个笛卡尔乘积的元素就是“1aB”,与上面的算法殊途同归。
4.再优化
仔细阅读上面的算法描述,你会发现算法的内层循环存在重复的字典元素拷贝,比如笛卡尔乘积元素下标0~25对应的字典[0-9]和[a-z]内的元素下标始终是0,那么就重复拷贝了[0-9]和[a-z]中元素25次。针对该问题,可以对上面的算法做进一步的优化。
以一次字典元素拷贝作为基本操作, 那么第二小节和第三小节的时间复杂度是O(hn),h为笛卡尔乘积空间大小,n为字典个数。
再优化算法描述如下:
再优化步骤描述如下:
(1)选取高度最高的字典SkS_k;
(2)循环h次,h为其它字典高度的乘积;
(2.1)将其它字典元素拼接在一起;
(2.2)循环最高字典高度HkH_k次,k为最高字典的下标,将元素填充到临时字符串s中后,将s加入笛卡尔乘积集合。
时间复杂度为O(h0h_0(n-1)+h0h1h_0h_1)=O(h0h_0(h1+n−1h_1+n-1))。其中h0h_0为非最高字典高度乘积,h1h_1为最高字典高度,n为字典个数,n≥2n\ge2。上文中未优化的时间复杂度是优化后的倍数t=h0h1nh0(h1+n−1)=n1+n−1h1t=\frac{h_0h_1n}{h_0(h_1+n-1)}=\frac{n}{1+\frac{n-1}{h_1}},可见,n和h1h_1越大,优化效果越明显。假设h1h_1很大,那么优化的倍数大约为字典的个数。
OpenMP并行加速笛卡尔乘积相关推荐
- JavaScript实现cartesianProduct笛卡尔乘积算法(附完整源码)
JavaScript实现cartesianProduct笛卡尔乘积算法(附完整源码) cartesianProduct.js完整源代码 cartesianProduct.test.js完整源代码 ca ...
- MySQL连接查询—笛卡尔乘积
连接查询 含义:连接查询又称为多表查询,当查询的字段来自多个表时,就会用到连接查询. 笛卡尔乘积 笛卡尔乘积现象:表1有 m 行,表2有 n 行,则查询结果有 m*n 行. SELECT `name` ...
- 使用Java8新特性(stream流、Lambda表达式)实现多个List 的笛卡尔乘积 返回需要的List<JavaBean>
需求分析: 有两个Long类型的集合 : List<Long> tagsIds; List<Long> attributesIds; 现在需要将这两个Long类型的集合进行组合 ...
- 一文完整MySQL连接查询,笛卡尔乘积,内连接外连接交叉连接
文章目录 笛卡尔乘积 连接查询分类 等值连接 非等值连接 自连接 外连接 交叉连接 连接查询又称为多表查询,当查询的字段来自于多个表时,使用连接查询. 笛卡尔乘积 笛卡尔乘积现象:表1有m行,表2有n ...
- 笛卡尔乘积算法的体现
一个商品有多项属性,多项属性的组合就产生不同的商品型号如: 衣服: 颜色:红.绿 尺寸:大.中.小 产生的系列就有:红|大.红|中.红|小.绿|大.绿|中.绿|小 如果商品的属性不至两 ...
- 克罗内克积(Kronecker product)与笛卡尔乘积(Cartesian product)
数学上,克罗内克积是两个任意大小的矩阵间的运算. 假设A为mxn的矩阵,B为pxq的矩阵. ,. 在进行计算的时候,克罗内克积可表示为,为mpxnq大小的矩阵.在MATLAB中可以使用kron函数. ...
- 【Python】求笛卡尔乘积
在数学中,两个集合X和Y的笛卡尓乘积(Cartesian product),又称直积,表示为X × Y,第一个对象是X的成员而第二个对象是Y的所有可能有序对的其中一个成员.假设集合A={a,b},集合 ...
- mysql 迪卡尔运算_「笛卡尔乘积」mysql笛卡儿积 - seo实验室
笛卡尔乘积 1.笛卡儿积数学理解 A={2,3} :B= {0,2,4} AxB = {(2,0),(3,0),(2,2),(3,2),(2,4),(3,4)},则AxB的结果集既是笛卡儿积: 但其不 ...
- python 笛卡尔_Python实现求笛卡尔乘积的方法
本文实例讲述了Python实现求笛卡尔乘积的方法.分享给大家供大家参考,具体如下: 在数学中,两个集合X和Y的笛卡尓乘积(Cartesian product),又称直积,表示为X × Y,第一个对象是 ...
最新文章
- FMDB支持的事务类型
- VTK:图片之ImageIslandRemoval2D
- 授权码模式-获取令牌
- 利用SVN工具下载OpenCore代码
- mysql执行一条语句会加锁吗_一条简单的更新语句,MySQL是如何加锁的?
- 大圆距离matlab代码,python – cartopy:大圆距离线的更高分辨率
- kafka与Spring的集成
- Android Ndef Message解析
- 阶段5 3.微服务项目【学成在线】_day05 消息中间件RabbitMQ_5.RabbitMQ研究-入门程序-生产者...
- ubuntu下载Linux 内核,ubuntu下载linux内核源码
- 加工机械双探头高频读写器CK-FR102AN用户开发手册
- 基于Packet Tracer 的校园网络设计方案(计算机网络与课程实验)(一)——子网划分方案规划
- ICPC North Central NA Contest 2017 H.Zebras and Ocelots
- 单例模式(DCL、holder等)
- 庄懂21课 塞博小人
- 画圆怎么编程python,python实现画圆功能
- unity 延迟等待执行
- requir与include的区别
- TokenInsight 对话首席——钱包安全与发展
- Flex和Flash的区别和联系
热门文章
- 思科警告: Nexus 交换机受严重缺陷影响
- Kubernetes 竞争条件漏洞的发现过程简述
- 我摸鱼写的Java片段意外称霸Stack Overflow十年、征服6000多GitHub开源项目: 有bug!...
- 无人机在倾斜摄影时丢片的解决方案
- 【初学】python执行系统命令四种方法比较
- 测试工程师的核心竞争力----打卡第九天
- 学校中有老师和学生两类人,而在职研究生既是老师又是学生,对学生的管理和对教师的管理在他们身上都有体现。...
- 调查显示存储虚拟化技术尚需3到5年才能成熟
- 关联容器——map、set
- Linux Shell 编程笔记