-
本版专家分：3607结帖率 100% -
有两个数值段列表, 例如:
A列表: 1 ~ 5, 10 ~ 16, 19 ~ 23, 34 ~ 50 共4个数值段
B列表: 3 ~ 7, 9 ~ 13, 21 ~ 33 共3个数值段
求两个列表的相交后的数值段, 结果应该为:
结果: 3 ~ 5, 10 ~ 13, 21 ~ 23 的这样一个列表
有什么高效算法来计算这个结果?
-
本版专家分：8810 -
用2个map保存数据， 轮询其中一个，再查找另外一个map中有没有相等的元素
-
本版专家分：5761
-
typedef struct
{
int a,b;
}INTPARE;
void GetIntersect(CArray<INTPARE,INTPARE> &TableA,CArray<INTPARE,INTPARE> &TableB,CArray<INTPARE,INTPARE> &TableRes)
{
TableRes.RemoveAll();
INTPARE PareA,PareB;
int a,b,c,d;
for (int i=0;i<TableA.GetCount();i++)
{
PareA=TableA[i];
for (int j=0;j<TableB.GetCount();j++)
{
PareB=TableB[j];
a=PareA.a-PareB.a;
b=PareA.a-PareB.b;
c=PareA.b-PareB.a;
d=PareA.b-PareB.b;
if(a!=0) a/=abs(a);
if(b!=0) b/=abs(b);
if(c!=0) c/=abs(c);
if(d!=0) d/=abs(d);
if(a*b<=0 || b*c<=0 || c*d<=0)//有交叉
{
INTPARE pare;
pare.a=max(PareA.a,PareB.a);
pare.b=min(PareA.b,PareB.b);
TableRes.Add(pare);
}
}
}
};
-
本版专家分：3607
-
非常感谢提供的代码, 但这种方式效率不算高吧. 里面还有-,*和/运算, -可以忽略不计, 但* 和 /开销挺大的. 而且很明显, 只要发现B当前范围均比A当前范围要大的时候, 就应该停止B循环, 进入下一个A范围的处理了.
我之前的条件漏说了个, 是范围段都是从小到大排好顺序的. 加上这个条件, 还得要加入下次B的循环应该是从上次循环最后被修改的范围开始查找, 不需要每次都从头查找.
虽然说是这样说, 但写代码的时候, 却发现所说的方案很难实现, 变化非常多, 人眼看上去就容易得出结果, 但机器要去做的时候, 找不到循环点..
-
本版专家分：5761
-
去掉*/很简单：
a+=abs(a);
b+=abs(b);
c+=abs(c);
d+=abs(d);
if(a==0 || b==0 || c==0 || d==0)//有交叉
{
INTPARE pare;
pare.a=max(PareA.a,PareB.a);
pare.b=min(PareA.b,PareB.b);
TableRes.Add(pare);
}
-
本版专家分：5761
-
上面改的不对，再换个方法
-
本版专家分：5761
-
下面的完全满足你的要求了：
void GetIntersect(CArray<INTPARE,INTPARE> &TableA,CArray<INTPARE,INTPARE> &TableB,CArray<INTPARE,INTPARE> &TableRes)
{
TableRes.RemoveAll();
INTPARE PareA,PareB;
int a,b,c,d,nStartPoint=0;
for (int i=0;i<TableA.GetCount();i++)
{
PareA=TableA[i];
for (int j=nStartPoint;j<TableB.GetCount();j++)
{
PareB=TableB[j];
if(PareA.b<PareB.a)
{
break;
}
a=PareA.a-PareB.a;
b=PareA.a-PareB.b;
c=PareA.b-PareB.a;
d=PareA.b-PareB.b;
if(a==0 || b==0 || c==0 || d==0 || (a^b)>>31 || (b^c)>>31 || (c^d)>>31)//有交叉
{
INTPARE pare;
pare.a=max(PareA.a,PareB.a);
pare.b=min(PareA.b,PareB.b);
TableRes.Add(pare);
nStartPoint=j;
}
}
}
}
-
本版专家分：3607
-
非常感谢, 这个效率就高很多了, B是只会循环一次了. 只要在A循环里面加多一个判断条件, 假如很明显A剩下的范围已经超出了所有B的范围, 就应该跳出A的循环了, 因为已经没有相交部分.
看了你的例子后, 我也想到了一个新的办法, 目前用相交的的方式, 好像可以采用并运算的方式来简化过程.
像上述的例子:
A列表: 1 ~ 5, 10 ~ 16, 19 ~ 23, 34 ~ 50 共4个数值段
B列表: 3 ~ 7, 9 ~ 13, 21 ~ 33 共3个数值段
把B列表取补集, 变成 0 ~ 2, 8 ~ 8, 14 ~ 20, 34 ~ MAX, 这些就是B列表的补集.(范围是0 ~ MAX)
利用这个补集, 好像可以更使判断过程更简单.
我写好了也发上来.
-
本版专家分：6384
-
楼主继续
我旁观
-
本版专家分：3607
-
采用补集并运算的过程大概流程描述就是:
整个过程就是归纳为:
先用小端查找
查找到时
当前段为起点
判大端是否也在本段内
在本段内
小端跟小端合并, 大端跟大端合并, 即插入一个新段(假如刚好为边界值是不用分段)
不在本段内
本段的大端变成当前查找的小端作为大端.
用大端查找位置, 所经过的段假如大端比当前查找的大端要小, 直接删除
查找不到, 删除完毕即可.
查找得到, 根据是否刚好为边界来进行决定分段还是不分段.
查找不时
当前段为起点
用大端查找后段, 所经过的段假如大端比当前查找的大端要小, 直接删除
查找不到, 删除完毕即可
查找得到, 根据是否刚好为边界来进行决定分段还是不分段.
即查找到跟查找不到时, 是小了一个在本段内进行比较的过程.
由于有相应的业务条件, A 交 B 的结果要放回 A 当中, 所以内部也有相应的删除双向链表操作. 并且假如交数 B 当中不存在 X~MAX 段的话, 那么可以确保 X+1 是属于补集范围内的, 所以最后会有一个单值段的并运算.
代码如下:
BOOL VXIntersect( DWORD32* pVXRanges, ULONG ulCnt ) // 交运算
{
DWORD32 dwMin = 0;
DWORD32 dwMax;
DWORD32 dwOldMin;
BOOL bFind;
if (0 != ulCnt % 2 || 0 == ulCnt) // 不是2的倍数肯定是一个错误的范围数组, 不能够进行交运算
return FALSE;
ULONG i = 0;
m_pVXCursor = m_pVXHead;
do
{
dwMax = pVXRanges[i] - 1;
bFind = FALSE;
while (NULL != m_pVXCursor)
{
if (m_pVXCursor->dwVXMin > dwMin) // 查找不到对应段
break;
if (dwMin >= m_pVXCursor->dwVXMin && m_pVXCursor->dwVXMax >= dwMin) // 查找到对应段
{
bFind = TRUE;
break;
}
m_pVXCursor = m_pVXCursor->pVXNext;
}
if (TRUE == bFind && dwMax <= m_pVXCursor->dwVXMax) // 小端查找到, 并且大端也在本范围段内时
{
if (dwMin == m_pVXCursor->dwVXMin && dwMax == m_pVXCursor->dwVXMax) // 刚好复盖整个范围, 当前段删除即可
VXDelete();
else if (dwMin == m_pVXCursor->dwVXMin) // 小段刚好全部合并, 大段写进去即可
{
m_pVXCursor->dwVXMin = dwMax;
m_pVXCursor = m_pVXCursor->pVXNext;
}
else if (dwMax == m_pVXCursor->dwVXMax) // 大段刚好全部合并, 留下小段即可
{
m_pVXCursor->dwVXMax = dwMin;
m_pVXCursor = m_pVXCursor->pVXNext;
}
else // 不能够闭合, 要创建新的段
{
dwOldMin = m_pVXCursor->dwVXMin;
m_pVXCursor->dwVXMin = dwMax;
VXInsertBefore( dwOldMin, dwMin );
}
}
else
{
if (TRUE == bFind)
{
if (m_pVXCursor->dwVXMin == dwMin) // 刚好合并时, 删除当前段即可
VXDelete();
else // 不能够闭合, 把当前段的最大值修改即可
{
if (0 == dwMin)
{
printf("error! impossible be zeor of min here!");
system("pause");
}
m_pVXCursor->dwVXMax = dwMin - 1;
}
}
while (NULL != m_pVXCursor)
{
if (m_pVXCursor->dwVXMax < dwMax) // 只要小于大端, 就可以直接删除
{
VXDelete();
continue;
}
if (m_pVXCursor->dwVXMin > dwMax) // 当前块的小端比当前比较的大端还大, 代表不需要再继续处理了
break;
if (m_pVXCursor->dwVXMax >= dwMax)
{
if (m_pVXCursor->dwVXMax == dwMax) // 刚好闭合, 删除这个块
VXDelete();
else // 不能闭合, 修改最小值
m_pVXCursor->dwVXMin = dwMax + 1;
break;
}
m_pVXCursor = m_pVXCursor->pVXNext;
}
}
i ++;
dwMin = pVXRanges[i] + 1;
i ++;
}while (i < ulCnt);
if (0 != dwMin)
{
if (NULL == m_pVXCursor)
m_pVXCursor = &m_VXEnd;
if (dwMin == m_pVXCursor->dwVXMin) // 实现插入单个区间, 刚好为边界就直接加即可
m_pVXCursor->dwVXMin = dwMin + 1;
else // 非边界时要创建一个新的区域
{
dwOldMin = m_pVXCursor->dwVXMin;
m_pVXCursor->dwVXMin = dwMin + 1;
VXInsertBefore( dwOldMin, dwMin - 1 );
}
}
m_pVXCursor = m_pVXHead;
}
分支挺多, 流程是很简单的, 大部分分支是为了不乱去new和delete. 并且实现确保不会存在 2 ~ 4, 5 ~ 8 这样的可合并分段, 2 ~ 4, 5 ~ 8应该是直接为 2 ~ 8 分段的.
结贴了.
-
本版专家分：3607
-
之前的代码没有测试, 这个代码是经过测试的正确代码
带有默认的条件:
1. 范围总是按从小到大的顺序的.
2. 输入范围不会存在像 12 ~ 13, 13 ~ 14 或 12 ~ 13, 14 ~ 15的区域, 虽然不允许, 但还是可以运算, 只是运算出来的范围可能会存在 13 ~ 13, 14 ~ 14, 但实际是可以合并成 13 ~ 14的.
3. 当输入的范围不存在 X ~ MAX 段, 默认会拿最后一个范围的最大值+2 ~ MAX 作补充. 即 一个范围段为: 1 ~ 3, 4 ~ 6, 在运算时, 会默认为 1 ~ 3, 4 ~ 6, 8 ~ MAX 作为输入.
BOOL VXIntersect( DWORD32* pVXRanges, ULONG ulCnt ) // 交运算
{
DWORD32 dwMin = 0;
DWORD32 dwMax;
DWORD32 dwOldMin;
BOOL bFind;
if (0 != ulCnt % 2 || 0 == ulCnt) // 不是2的倍数肯定是一个错误的范围数组, 不能够进行交运算
return FALSE;
ULONG i = 0;
m_pVXCursor = m_pVXHead;
do
{
dwMax = pVXRanges[i] - 1;
bFind = FALSE;
while (NULL != m_pVXCursor)
{
if (m_pVXCursor->dwVXMin > dwMin) // 查找不到对应段
break;
if (dwMin >= m_pVXCursor->dwVXMin && m_pVXCursor->dwVXMax >= dwMin) // 查找到对应段
{
bFind = TRUE;
break;
}
m_pVXCursor = m_pVXCursor->pVXNext;
}
if (TRUE == bFind && dwMax <= m_pVXCursor->dwVXMax) // 小端查找到, 并且大端也在本范围段内时
{
if (dwMin == m_pVXCursor->dwVXMin && dwMax == m_pVXCursor->dwVXMax) // 刚好复盖整个范围, 当前段删除即可
VXDelete();
else if (dwMin == m_pVXCursor->dwVXMin) // 小段刚好全部合并, 大段写进去即可
{
m_pVXCursor->dwVXMin = dwMax + 1;
//m_pVXCursor = m_pVXCursor->pVXNext;
}
else if (dwMax == m_pVXCursor->dwVXMax) // 大段刚好全部合并, 留下小段即可
{
m_pVXCursor->dwVXMax = dwMin - 1;
m_pVXCursor = m_pVXCursor->pVXNext;
}
else // 不能够闭合, 要创建新的段
{
dwOldMin = m_pVXCursor->dwVXMin;
m_pVXCursor->dwVXMin = dwMax + 1;
VXInsertBefore( dwOldMin, dwMin - 1 );
}
}
else
{
if (TRUE == bFind)
{
if (m_pVXCursor->dwVXMin == dwMin) // 刚好合并时, 删除当前段即可
VXDelete();
else // 不能够闭合, 把当前段的最大值修改即可
{
if (0 == dwMin)
{
printf("error! impossible be zeor of min here!");
system("pause");
}
m_pVXCursor->dwVXMax = dwMin - 1;
m_pVXCursor = m_pVXCursor->pVXNext;
}
}
while (NULL != m_pVXCursor)
{
if (m_pVXCursor->dwVXMax <= dwMax) // 只要小于大端, 就可以直接删除
{
VXDelete();
continue;
}
if (m_pVXCursor->dwVXMin > dwMax) // 当前块的小端比当前比较的大端还大, 代表不需要再继续处理了
break;
if (m_pVXCursor->dwVXMax >= dwMax)
{
if (m_pVXCursor->dwVXMax == dwMax) // 刚好闭合, 删除这个块
VXDelete();
else // 不能闭合, 修改最小值
{
m_pVXCursor->dwVXMin = dwMax + 1;
//m_pVXCursor = m_pVXCursor->pVXNext;
}
break;
}
m_pVXCursor = m_pVXCursor->pVXNext;
}
}
i ++;
dwMin = pVXRanges[i] + 1;
i ++;
}while (i < ulCnt);
if (0 != dwMin)
{
if (NULL == m_pVXCursor)
m_pVXCursor = &m_VXEnd;
while (NULL != m_pVXCursor)
{
if (m_pVXCursor->dwVXMin > dwMin) // 查找不到对应段
break;
if (dwMin >= m_pVXCursor->dwVXMin && m_pVXCursor->dwVXMax >= dwMin) // 查找到对应段
{
if (dwMin == m_pVXCursor->dwVXMin && dwMin == m_pVXCursor->dwVXMax) // 刚好复盖整个范围, 当前段删除即可
VXDelete();
else if (dwMin == m_pVXCursor->dwVXMin) // 实现插入单个区间, 刚好为边界就直接加即可
m_pVXCursor->dwVXMin = dwMin + 1;
else if (dwMin == m_pVXCursor->dwVXMax) // 大段刚好全部合并, 留下小段即可
m_pVXCursor->dwVXMax = dwMin - 1;
else // 非边界时要创建一个新的区域
{
dwOldMin = m_pVXCursor->dwVXMin;
m_pVXCursor->dwVXMin = dwMin + 1;
VXInsertBefore( dwOldMin, dwMin - 1 );
}
break;
}
m_pVXCursor = m_pVXCursor->pVXNext;
}
}
m_pVXCursor = m_pVXHead;
return TRUE;
}
- java遗传算法_求函数值，java界面版
java遗传算法_求函数值，java界面版，求函数最值的最优化问题
- Delaunay三角剖分算法
Delaunay三角剖分算法 1. 三角剖分与Delaunay剖分的定义 如何把一个散点集合剖分成不均匀的三角形网格，这就是散点集的三角剖分问题，散点集的三角剖分，对数值分析以及图形学来说，都是极为重要的一项预处理技术。 ...
- 三维裂隙渗流模拟中裂隙对相交判断过程的优化
引入定向包围盒技术对传统相交分析"两步法"进行算法优化,提出运算成本估算公式和基于假想裂隙面的参数选取方法。经多组相交实验表明,在大规模裂隙模拟中采用定向包围盒的整体运算时间与范围盒和轴向包围盒相比可节约...
- 《计算机图形学原理及算法教程》（Visual C++版）电子教案（课件）+配书源代码
<br>详细了解计算机图形学原理及算法教程（Visual C++版） <br>目录： <br>第1章 基本图形的生成 1 1.1 直线 1 1.1.1 DDA（数值微分）算法 2 1.1.2 生成直线的中点画线法 2 1.1.3 Bresenham...
- 牛顿迭代求根算法的分析与实现 论文 完整版
对这些问题，都可以利用数值方法来求解，在计算机中实现的数值方法也称为数值算法。牛顿迭代法是数值分析中一个重要的计算方法和思想。迭代法的主要功能：计算方程时可以比较快速。 在工程实践中，有许多问题往往...
- 关于二维O（N）玻色子模型的S矩阵自举的注释
这是应用数学中的一个标准问题，与半定规划有关，可以通过快速可用的数值算法求解。 数值解提供的信息足以在不使用可积性的情况下重现已知的解析解，即Yang-Baxter方程。 这种情况似乎很普遍，因此我们希望通过最大...
- GIS面试题,希望大家喜欢
对于任意两个矩形，空间上有包含、相离、相交三种空间关系，现要求设计一个算法判断两个矩形的空间关系，如果两个矩形相交，计算相交部分。 4. 某国土局需要建设一个综合管理信息平台以满足国土资源管理的业务需要...
- 3D游戏引擎设计实时计算机图形学的应用方法第2版
实时计算机图形学的应用方法（第2版）》详细阐述了与3D游戏引擎设计相关的高效解决方案及相应的数据结构和算法，主要包括图形系统、渲染器、场景图、控制器动画、空间排序、细节级别、碰撞检测、物理学、标准对象、...
- 浙江大学ACM模板 计算几何，图论，数据结构，经典题的模板
1.16.9 pku_1151_矩形相交的面积 76 1.16.10 pku_1118_共线最多的点的个数 78 1.16.11 pku2826_线段围成的区域可储水量 80 1.16.12 Pick公式 84 1.16.13 N点中三个点组成三角形面积最大 86 1.16.14 直线关于圆的反射...
- 第十一节 图像处理之霍夫检测直线
y2 = int(y0 - 1000 * a) #计算直线终点纵坐标 注：这里的数值1000给出了画出的线段长度范围大小，数值越小，画出的线段越短，数值越大，画出的线段越长 cv.line(image, (x1, y1), (x2, y2), (0, 0, 255), 2) #点...
- MAPGIS地质制图工具
剪断相交线：选择线剪断相交线或者直接拉一条线剪断相交线，自相交及节点断线。 等分线：定距离或定段数等分一条线。还可以随机按偏移量自动等分线。 图元筛选：根据参数或者属性筛选图元，还可选择相同的图元。 ...
- WGAN
掌握基于深度学习的计算机视觉原理与算法 掌握基于深度学习的计算机视觉的项目实践方法 提高深度学习计算机视觉方向的文献阅读能力
- 68_直线与三角形的交点C++实践1of2
掌握游戏编程和计算机图形学的基本数学知识。 掌握线性代数、几何变换、运动学、3D物理和相关数值运算的基本方法。
- 69_直线与三角形的交点C++实践2of2
掌握游戏编程和计算机图形学的基本数学知识。 掌握线性代数、几何变换、运动学、3D物理和相关数值运算的基本方法。
- [excel2010应用大全].ExcelHome.扫描版.pdf
14.5.1 数字条码校验码的模拟算法 14.5.2 计扣个人所得税 14.5.3 专利发明奖励统计 14.5.4 数字校验应用 14.5.5 返回指定有效数字 第15章 日期和时间计算 15.1 认识日期数据 15.1.1 了解日期序数 15.1.2 了解时间...
- 第二周直播
第一阶段：机器学习原理及商品推荐系统实现（5 weeks） 第二阶段：深度学习原理及实战项目强化训练（6 weeks） 第三阶段：四个工业级实战项目（可选）及成果展示（31 days）
- 土豆SDK(Java版)-非官方
由于土豆SDK一直建设中，最近几天抽空写了一套java的SDK。包含了现有的所有请求协议。本套SDK中仅提供了oAuth的方式(引用oAuth.net的java版示例)，并没有在框架中实现，涉及到登录验证的，就麻烦自己去完善框架了。
- MFC上位机与STM32下位机通讯精讲
本课程主要介绍C++类库MFC上位机与STM32单片机的RS232、RS422、RS485、USB、LWIP以太网、CAN等接口进行稳定通信。课程主要从MFC和STM32基础开始，以编写上位机以及下位机为主，非常注重实践。 本课程主要目标是让大家学会MFC上位机的编程、STM32下位机编程。单单的会上位机或者下位机有时候体会不到通讯的乐趣，如果你单单会下位机，你学了本课程，你可以把你的数据以任何方式在windows上呈现出来，如果你只会上位机，你学了本课程你可以通过windows的串口、以太网、USB控制任何硬件设备。本课程的目标就是补偿大家的短处，让大家成为综合性人才，同时让大家体会到通讯中的更多乐趣。
- QT漂亮界面
QT界面开发。QT漂亮界面。 欢迎大家一起学习。。。。。
- navicat简体中文版 绿色版 （64位）
解压后安装navicat，打开navicat执行PatchNavicat即破解成功。可以正常使用啦。