我需要浏览一个包含地图中点数的数组,并检查它们之间的距离。我需要计算每个节点在200米和50米范围内的节点数。它适用于较小数量的值。但是,当我尝试通过它运行更多值时(大约4000用于可伸缩性测试),发生错误,表示我已达到300秒的最大执行时间。如果可能的话,它需要能够在300秒内至少处理这么多。
我已阅读并发现有一种方法可以禁用/更改此限制,但是我想知道是否有更简单的方法来执行以下代码以便运行所需的时间它会减少。
for($i=0;$i<=count($data)-1;$i++)
{
$amount200a=0;
$amount200p=0;
$amount50a=0;
$amount50p=0;
$distance;
for($_i=0;$_i<=count($data)-1;$_i++)
{
$distance=0;
if($data[$i][0]===$data[$_i][0])
{
}
else
{
//echo "Comparing ".$data[$i][0]." and ".$data[$_i][0]." ";
$lat_a = $data[$i][1] * PI()/180;
$lat_b = $data[$_i][1] * PI()/180;
$long_a = $data[$i][2] * PI()/180;
$long_b = $data[$_i][2] * PI()/180;
$distance =
acos(
sin($lat_a ) * sin($lat_b) +
cos($lat_a) * cos($lat_b) * cos($long_b - $long_a)
) * 6371;
$distance*=1000;
if ($distance<=50)
{
$amount50a++;
$amount200a++;
}
else if ($distance<=200)
{
$amount200a++;
}
}
}
$amount200p=100*number_format($amount200a/count($data),2,'.','');
$amount50p=100*number_format($amount50a/count($data),2,'.','');
/*
$dist[$i][0]=$data[$i][0];
$dist[$i][1]=$amount200a;
$dist[$i][2]=$amount200p;
$dist[$i][3]=$amount50a;
$dist[$i][4]=$amount50p;
//*/
$dist.=$data[$i][0]."&&".$amount200a."&&".$amount200p."&&".$amount50a."&&".$amount50p."%%";
}
索引0包含每个节点的唯一ID,1包含每个节点的纬度和 index 2包含每个节点的经度。
错误发生在第一个循环内的第二个for循环中。此循环是将所选映射节点与其他节点进行比较的循环。我也在使用Haversine Formula。
答案 0 :(得分:0)
首先,你正在以大O表示法执行:O(数据^ 2),这将会很慢,实际上,要么有两种可能的解决方案。找到一种经过验证的算法,可以在更好的时间内解决相同的问题。或者,如果你不能,开始移动内在的for循环,并通过数学证明你是否可以将内部for循环转换为大多数简单的计算,这是经常你可以做的事情。
经过一些改写后,我看到了一些可能性: 如果$ data不是SPLFixedArray(具有FAR更好的访问时间),那么就进行。,因为您访问该数据的次数很多(4000 ^ 2)* 2。 更好,写清洁代码。虽然optizmier会尽力而为,但如果你不尝试缩小代码(这只会让它更具可读性),那么它可能无法做到尽可能好。并将中间结果移出循环,也就像数组的大小一样。
答案 1 :(得分:0)
目前,您正在针对所有其他点检查所有积分,实际上您只需要针对所有剩余点检查当前点。从A到B的距离与从B到A的距离相同,那么为什么要计算两次?
我可能会创建一个相邻的数组来计算在彼此范围内有多少个节点,并在我计算出两个节点在彼此的范围内之后增加该数组中的条目对。
你应该想出一个非常快速的距离近似值,可以在计算真实距离之前忽略尽可能多的节点(永远不会超快)。
一般来说,除了算法优化之外,优化的基本规则是:
不要做任何你不必做的处理:比如不要将$ distance乘以1000.只需将你测试的值从20和50改为0.02和0.05。
不要经常调用任何函数:只需要在任何处理开始之前调用count($ data)一次。
不要多次计算常数值:PI()/180
,例如。
将所有可能的处理移到循环之外。即尽可能预先计算。
另一个小问题,它将使您的代码更容易阅读:
for( $i = 0; $i <= count( $data ) - 1; $i++ )
与:
for( $i = 0; $i < count( $data ); $i++ )
答案 2 :(得分:0)
试试这个:
$max = count($data);
$CONST_PI = PI() / 180;
for($i=0;$i<$max;$i++)
{
$amount200a=0;
$amount50a=0;
$long_a = $data[$i][2] * $CONST_PI;
$lat_a = $data[$i][1] * $CONST_PI;
for($_i=0;$_i<=$max;$_i++)
//or use for($_i=($i+1);$_i<=$max;$_i++) if you did not need to calculate already calculated in other direction
{
$distance=0;
if($data[$i][0]===$data[$_i][0]) continue;
$lat_b = $data[$_i][1] * $CONST_PI;
$long_b = $data[$_i][2] * $CONST_PI;
$distance =
acos(
sin($lat_a ) * sin($lat_b) +
cos($lat_a) * cos($lat_b) * cos($long_b - $long_a)
) * 6371;
if ($distance<=0.2)
{
$amount200a++;
if ($distance<=0.05)
{
$amount50a++;
}
}
} // for %_i
$amount200p=100*number_format($amount200a/$max,2,'.','');
$amount50p=100*number_format($amount50a/$max,2,'.','');
$dist.=$data[$i][0]."&&".$amount200a."&&".$amount200p."&&".$amount50a."&&".$amount50p."%%";
} // for $i
我认为阅读会更好,如果你改变$ _i的注释行,它会更快:)