格点统计

题目大意

第一象限(由于本人思维不算太敏捷,不懂第一象限的只能给你们链接自己看)中,位于反比例函数xy=k的下方(含边界)的格点的个数。

思路详解

这道题,我们先将第一象限画出来,牛逼的可以脑中构图。


如上图,我画了一个k为4的第一象限图(无视边上的点)于是我们很成功的发现下表

(1,1)(2,1)(3,1)(4,1)
(2,1)(2,2)  
(3,1)   
(4,1)   

求出来有什么用呢?不用急,一步步来。

  1. 我们先按上表打一个为O(n)的程序,一层层做。这样我们的时间复杂度就从O(n^2)--->O(n)。
  2. 再仔细观察,我们可以发现,如纵列[(1,1),(2,1),(3,1),(4,1)]与横列[(1,1),(2,1),(3,1),(4,1)]把个数相加,重复了左上角的(1,1)。所以我们可以一层层这么枚举求出规律。首先定义一个len。len=k div i-i+1。这么做便可以求出被前面计算过的数。判断一下,如果这个数大于1,说明除(i,i)点(如(1,1)(2,2)这样的在表中位于左上角的点)外还有数,便将它*2,便能一起求出表这一列的数。但注意,要不能将(i,i)点也乘2,这样会"wrong answer"。这样时间复杂度便可优化到O(trunc(sqrt(k)))。

r
        i,j:longint;
        k,ans,len:int64;
begin
        assign(input,'count.in'); reset(input);
        assign(output,'count.out'); rewrite(output);
        readln(k);
        for i:=1 to trunc(sqrt(k)) do
        begin
                len:=k div i-(i-1);
                if (len>1) then len:=(len-1)*2+1;
                ans:=(ans+len) mod 998244353;
        end;
        writeln(ans mod 998244353);
end.

转载于:https://www.cnblogs.com/Sport-river/p/10390138.html

<think>好的,用户的问题是如何在QGIS中统计内的轨迹密度。我需要逐步思考解决方法,并确保回答结构清晰且正确。 首先,用户提到已经将国家数据存入PostGIS,但问题转向轨迹密度分析。可能需要先明确数据的准备情况:用户是否有轨迹的矢量数据,比如图层,以及是否已经创建了网图层。假设用户有轨迹数据,但需要创建网。 第一步,创建网。QGIS有创建网的工具,位于矢量菜单下的研究工具中。可以选择矩形或六边形网,设置范围和间距,输出为矢量图层。用户需要根据轨迹的分布范围调整参数,确保网覆盖所有轨迹。 第二步,将轨迹空间连接。使用“按位置连接”工具,将的属性连接到网面,统计每个网内的数。这里需要确认连接类型和统计字段的设置,确保正确计数。 第三步,计算密度。添加字段,用数除以网面积。需要注意单位的一致性,比如将面积从平方米转换为平方公里,或者根据需求调整。可能需要在字段计算器中使用表达式,例如数 / (面积 / 1000000) 得到每平方公里的密度。 另外,如果用户的数据存储在PostGIS中,可能需要使用空间查询。例如,使用ST_Intersects连接网和轨迹,用COUNT聚合函数统计每个网数,再用ST_Area计算面积。然后将结果导出到QGIS进行可视化。 可视化方面,建议用户对密度字段进行分级符号或热力图渲染,以便直观展示密度差异。同时,可能需要提醒用户检查数据投影,确保面积计算准确,以及网和轨迹的坐标系统一致。 用户可能需要进一步的优化,比如处理大量数据时的性能问题,或者动态计算密度的方法。需要简要提及这些可能性,但主要解决当前问题。</think>在QGIS中统计内轨迹密度的操作步骤如下,结合PostGIS数据库和QGIS空间分析工具实现: ### 一、数据准备阶段 1. **创建渔网(网)** - 使用QGIS菜单:`矢量 > 研究工具 > 创建网` - 选择网类型(矩形/六边形),设置范围参数时建议勾选`使用图层范围`并选择轨迹图层 - 网间距参数示例:若需要1km×1km网,在投影坐标系下输入1000(单位:米) 2. **建立空间索引** ```sql -- 对轨迹表创建空间索引 CREATE INDEX ON trajectory_points USING GIST(geom); -- 对网表创建空间索引 CREATE INDEX ON grid USING GIST(geom); ``` ### 二、密度计算阶段 #### 方法1:QGIS原生工具 1. **执行空间连接** ```python # 使用Processing Toolbox工具: # Vector > Analysis > Count points in polygon # 输入:网图层,轨迹图层 # 输出字段命名为point_count ``` 2. **计算单位面积密度** ```sql -- 在网图层字段计算器中: density = point_count / ($area / 1000000) -- 转换为平方公里密度 ``` #### 方法2:PostGIS空间查询 ```sql SELECT grid.id, COUNT(traj.geom) AS point_count, COUNT(traj.geom)/ST_Area(grid.geom) AS density FROM grid LEFT JOIN trajectory_points AS traj ON ST_Intersects(grid.geom, traj.geom) GROUP BY grid.id; ``` ### 三、可视化增强 1. **分级渲染设置** - 右键网图层选择`属性 > 符号化` - 选择`分级符号`模式,值字段选择density - 调整色带和分类方法(自然断/等间隔) 2. **叠加验证** ```python # 将计算结果导出为GeoJSON # 使用Qgis的Layer > Add Layer > Add/Edit Virtual Layer功能 # 叠加原始轨迹验证分布一致性 ``` ### 四、注意事项 1. 坐标系必须使用投影坐标系(如UTM),否则面积计算会出现偏差[^1] 2. 处理千万级数据时建议在PostGIS中完成计算 3. 可通过创建物化视图提升查询性能: ```sql CREATE MATERIALIZED VIEW density_map AS [上述查询语句]; ```
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符  | 博主筛选后可见
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值