保留道路

这里写图片描述
这里写图片描述
这里写图片描述
这里写图片描述
这里写图片描述
将边排序,以g为第一关键字,s为第二关键字
每次枚举一个maxg,将g

if(a[i].s+a[i].g>ans) continue;

如果当前ans<当前加入的边的g+s,那么这条边对于ans无任何贡献(因为在边集中的边的s的最大值(加入此边后)一定>=当前S),那么这些边生成的最小生成树所得出ans一定比当前ans大,无意义
不加60,加了AC

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
struct node{
    int x,y;
    long long g,s;
    bool vis;
}a[510000];
node now_tree[520000],tree[510000];
int fat[510000];
int n,m,max_g,max_s,cnt;
long long ans=1e17;
inline int read()
{
    int data=0,w=1; char ch=0;
    while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
    if(ch=='-') w=-1,ch=getchar();
    while(ch>='0' && ch<='9') data=data*10+ch-'0',ch=getchar();
    return data*w;
}
bool comp(const node&s1,const node&s2)
{
    if(s1.g!=s2.g) return s1.g<s2.g;
    return s1.s<s2.s;
}
int find(int x)
{
    if(fat[x]==x) return x;
    return fat[x]=find(fat[x]);
}
void init()
{
    for(int i=1;i<=n;i++)
     fat[i]=i;
    for(int i=1;i<=cnt;i++)
     now_tree[i]=tree[i],now_tree[i].vis=0;
}
void get_ans(long long maxg)
{
    long long maxs=0;
    long long tot=0;
    init();
    for(int i=1;i<=cnt;i++)
     {
        int f1=find(tree[i].x),f2=find(tree[i].y);
        if(f1!=f2)
        {
            tot++;
            fat[f1]=f2;
            now_tree[i].vis=1;
            maxs=max(maxs,tree[i].s);
        }
        if(tot==n-1)
        {
            int sum=0;
            for(int j=1;j<=cnt;j++)
             if(now_tree[j].vis) tree[++sum]=now_tree[j];
            cnt=sum;
            ans=min(ans,maxs+maxg);
            break;
        }
     }
     return;
}
int main()
{
    freopen("road.in","r",stdin);
    freopen("road.out","w",stdout);
    scanf("%d%d%d%d",&n,&m,&max_g,&max_s);
    for(int i=1,g,s;i<=m;i++)
    {
        a[i].x=read(),a[i].y=read();
        g=read(),s=read();
        a[i].g=(long long)g*max_g,a[i].s=(long long)s*max_s;
    }
    sort(a+1,a+m+1,comp);
    int flag=0;
    for(int i=1;i<=m;i++)
     {
        if(a[i].s+a[i].g>ans) continue;
        int w=cnt+1;
        for(int j=1;j<=cnt;j++)
         if(tree[j].s>a[i].s)
          {
            w=j;
            break;
          } 

        if(w==cnt+1) tree[++cnt]=a[i];
        else
        {
          cnt++;
          for(int j=cnt;j>=w+1;j--)
           tree[j]=tree[j-1];
          tree[w]=a[i];
        }
        get_ans(a[i].g);
     } 
     if(ans!=1e16) 
     printf("%lld",ans);
     else
     printf("-1");
}
针对如何将广东省道路矢量据从Shp格式转换为GeoJSON格式,并保留道路分类属性及属性信息的问题,您可参考《广东省道路矢量层及Shp转GeoJSON代码发布》提供的资源。该资源中包含的Python代码是一个非常实用的工具,能够帮助您实现上述需求。具体步骤如下: 参考资源链接:[广东省道路矢量层及Shp转GeoJSON代码发布](https://wenku.csdn.net/doc/7cqyirghsw) 首先,确保您的环境中安装了必要的Python库,如Fiona或GDAL/OGR,这些库能够帮助您读取和写入Shp以及GeoJSON文件。接下来,您可以使用以下代码结构进行转换操作: ```python import fiona from shapely.geometry import shape # 打开Shp文件 with fiona.open('道路路网矢量据.shp', 'r') as source: # 创建一个新的GeoJSON文件 schema = source.schema.copy() schema['geometry'] = 'Point' with open('道路路网矢量据.geojson', 'w', encoding='utf-8') as dest: # 遍历Shp文件中的每条记录 for item in source: # 将Shp的几何形状转换为GeoJSON的几何形状 geom = shape(item['geometry']) # 创建GeoJSON的feature feature = { 'type': 'Feature', 'id': item['id'], 'properties': item['properties'], 'geometry': geom.__geo_interface__ } # 将feature写入GeoJSON文件 dest.write(feature) # 检查属性信息是否完整 with open('道路路网矢量据.geojson', 'r') as *** *** ***['features']: if '道路分类属性' in feature['properties']: print('属性信息完整保存。') ``` 在上述代码中,我们使用了Fiona库来读取和写入Shp及GeoJSON文件,同时利用shapely库转换几何形状。这样可以确保在转换过程中,不仅空间据得到保留道路的多级分类属性及其它属性信息也同样得到完整转换。 完成上述步骤后,您将获得一个包含广东省道路路网矢量据的GeoJSON文件,该文件既方便在Web端进行展示和分析,也可以作为进一步GIS分析的据源。此外,如果您希望深入了解GIS据转换的更多细节和高级应用,建议详细阅读《广东省道路矢量层及Shp转GeoJSON代码发布》资源,以便更加全面地掌握GIS据处理的技巧。 参考资源链接:[广东省道路矢量层及Shp转GeoJSON代码发布](https://wenku.csdn.net/doc/7cqyirghsw)
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值