NOIP2013 DAY1 T3

本文介绍了一种结合Kruskal算法和LCA算法的解决方案,用于求解最大生成树及两点间的最小权值路径问题。通过具体实现代码展示了如何使用这两种算法,并提供了完整的C++代码示例。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目链接

Kruskal+LCA求解
Kruskal求最大生成树
LCA查找两点间最小权值

可以参见另一题
星际导航

上代码
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MIN(a,b) a<b?a:b
using namespace std;

int n,m,fa[500000],khead[500000],head[500000],edge_num,t_num,deep[100009],anc[100009][21];
int weight[100009][21],kruedge_num,q;
struct E{
    int next,to,value,from;
}edge[500000],kruedge[500000];

inline void Swap(int &a,int &b){
    int t=a;
    a=b;b=t;
}

bool cmp(const E &x,const E &y){
    return x.value>y.value;
}
int Find(int x){
    if(fa[x]!=x) return fa[x]=Find(fa[x]);
    return x;
}
inline int Read(){
    int x=0,flag=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') flag=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*flag;
}

void addedge(int x,int y,int z){
    edge[++edge_num].next=head[x];
    edge[edge_num].from=x;
    edge[edge_num].to=y;
    edge[edge_num].value=z;
    head[x]=edge_num;
}

void addedgekru(int x,int y,int z){
    kruedge[++kruedge_num].next=khead[x];
    kruedge[kruedge_num].from=x;
    kruedge[kruedge_num].to=y;
    kruedge[kruedge_num].value=z;
    khead[x]=kruedge_num;
}

void Kruskal(){
    sort(edge+1,edge+m+1,cmp);
    int i,k=0;
    for(i=1;i<=n;i++){
        fa[i]=i;
    }
    for(i=1;i<=m;i++){
        int a=edge[i].from,b=edge[i].to,v=edge[i].value;
        int aa=Find(a),bb=Find(b);
        if(aa!=bb){
            fa[aa]=bb;
            addedgekru(a,b,v);
            addedgekru(b,a,v);
            k++;
            if(k==n-1) return;
        }
    }
}

void DFS(int x){
    int i;
    for(i=khead[x];i;i=kruedge[i].next){
        if(kruedge[i].to==anc[x][0]) continue;
        anc[kruedge[i].to][0]=x;
        deep[kruedge[i].to]=deep[x]+1;
        weight[kruedge[i].to][0]=kruedge[i].value;
        DFS(kruedge[i].to);
    }
}

void PRE(){
    int i,j;
    for(i=1;i<=n;i++){
        if(!anc[i][0]){
            deep[i]=1;
            DFS(i);
        }
    }
    for(j=1;(1<<j)<=n;j++){
        for(i=1;i<=n;i++){
            if(anc[i][j-1]){
                anc[i][j]=anc[anc[i][j-1]][j-1];
                weight[i][j]=MIN(weight[i][j-1],weight[anc[i][j-1]][j-1]);
            }
        }
    }
}

int LCA(int x,int y){
    int ans=2e9;
    if(deep[x]<deep[y]) Swap(x,y);
    int i;
    int maxlog=floor(log(n)/log(2));
    for(i=maxlog;i>=0;i--){
        if(deep[x]-(1<<i)>=deep[y]){
            ans=MIN(ans,weight[x][i]);
            x=anc[x][i];
        }
    }
    if(x==y) return ans;

    for(i=maxlog;i>=0;i--){//
        if(anc[x][i] && anc[x][i]!=anc[y][i]){
            ans=MIN(ans,weight[x][i]);ans=MIN(ans,weight[y][i]);
            x=anc[x][i];y=anc[y][i];
        }
    }
    ans=MIN(ans,weight[x][0]);ans=MIN(ans,weight[y][0]);
    return ans;
}

int main(){
    freopen("truck.in","r",stdin);
    freopen("truck.out","w",stdout);
    n=Read();m=Read();
    int i;
    for(i=1;i<=m;i++){
        int a=Read(),b=Read(),c=Read();
        addedge(a,b,c);
        //addedge(b,a,c);
    }
    Kruskal();
    PRE();
    q=Read();
    for(i=1;i<=q;i++){
        int a=Read(),b=Read();
        int aa=Find(a),bb=Find(b);
        if(aa!=bb) printf("-1\n");
        else printf("%d\n",LCA(a,b));
    }
    return 0;
}
内容概要:本文探讨了在MATLAB/SimuLink环境中进行三相STATCOM(静态同步补偿器)无功补偿的技术方法及其仿真过程。首先介绍了STATCOM作为无功功率补偿装置的工作原理,即通过调节交流电压的幅值和相位来实现对无功功率的有效管理。接着详细描述了在MATLAB/SimuLink平台下构建三相STATCOM仿真模型的具体步骤,包括创建新模型、添加电源和负载、搭建主电路、加入控制模块以及完成整个电路的连接。然后阐述了如何通过对STATCOM输出电压和电流的精确调控达到无功补偿的目的,并展示了具体的仿真结果分析方法,如读取仿真数据、提取关键参数、绘制无功功率变化曲线等。最后指出,这种技术可以显著提升电力系统的稳定性与电能质量,展望了STATCOM在未来的发展潜力。 适合人群:电气工程专业学生、从事电力系统相关工作的技术人员、希望深入了解无功补偿技术的研究人员。 使用场景及目标:适用于想要掌握MATLAB/SimuLink软件操作技能的人群,特别是那些专注于电力电子领域的从业者;旨在帮助他们学会建立复杂的电力系统仿真模型,以便更好地理解STATCOM的工作机制,进而优化实际项目中的无功补偿方案。 其他说明:文中提供的实例代码可以帮助读者直观地了解如何从零开始构建一个完整的三相STATCOM仿真环境,并通过图形化的方式展示无功补偿的效果,便于进一步的学习与研究。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值