NKOI 2747 最小生成树1

本文深入探讨了最小生成树的概念,并详细介绍了Kruskal和Prim算法的实现过程。通过实例代码,展示了如何利用这两种算法解决最小生成树问题,包括输入输出规范和算法的时间复杂度分析。

最小生成树1

Time Limit:10000MS  Memory Limit:234567K
Total Submit:108 Accepted:54
Case Time Limit:1000MS

Description

给定一个N个节点M条边的带权无向图,节点从1到N开始编号,求最小生成树

Input

第一行两个数N,M,表示点数和边数
接下来M行每行三个数u,v,w,表示u,v之间有一条权为w的边
可能有重边,但是没有自环

Output

如果最小生成树存在,输出树上所有边的权值和
如果不存在,输出“Keng Die!”

Sample Input

4 5
1 2 1
2 3 2
3 4 3
4 1 4
1 3 5

Sample Output

6

Hint

1<=N<=200000
N-1<=M<=500000
0<=每条边的权值<=5000

Source



先用 kruscal:

#include <cstdio>
#include <cstdlib>
using namespace std;
const int maxn=200005,maxm=500005;
int fath[maxn],x[maxm],y[maxm],z[maxm],N,M;
inline void _getnum(int &xx)
{
    char tt=getchar();
    while(tt<'0'||tt>'9')tt=getchar();
    for(xx=0;tt>='0'&&tt<='9';tt=getchar())xx=xx*10+(tt-'0');
}
void _qst(int l,int r)
{
    int i,j,m,t;
    i=l;j=r;
    m=z[(i+j)>>1];
    while(i<=j)
    {
        while(z[i]<m)i++;
        while(z[j]>m)j--;
        if(i<=j)
        {
            t=x[i];x[i]=x[j];x[j]=t;
            t=y[i];y[i]=y[j];y[j]=t;
            t=z[i];z[i]=z[j];z[j]=t;
            i++;j--;
        }
    }
    if(i<r)_qst(i,r);
    if(l<j)_qst(l,j);
}
int _getfath(int i)
{
    if(fath[i]!=i)fath[i]=_getfath(fath[i]);
    return fath[i];
}
int main()
{
    int i,j,ans=0,X,Y;
    _getnum(N);_getnum(M);
    for(i=1;i<=M;i++)_getnum(x[i]),_getnum(y[i]),_getnum(z[i]);
    _qst(1,M);
    for(i=1;i<=N;i++)fath[i]=i;
    for(i=j=1;i<=M&&j<N;i++)
    {
        X=_getfath(x[i]);
        Y=_getfath(y[i]);
        if(X!=Y)fath[X]=Y,j++,ans+=z[i];
    }
    if(j==N)printf("%d",ans);
    else printf("Keng Die!");
    return 0;
}
然后是prim:

#include <cstdio>
#include <cstdlib>
using namespace std;
const int maxn=200005,maxm=1000005;
int N,M,dis[maxn],H[maxn],pH[maxn],tot;
int y[maxm],z[maxm],next[maxm],last[maxn];
inline void _getnum(int &xx)
{
    char tt=getchar();
    while(tt<'0'||tt>'9')tt=getchar();
    for(xx=0;tt>='0'&&tt<='9';tt=getchar())xx=xx*10+(tt-'0');
}
void _shift_up(int i)
{
    int j=i>>1,t=H[i];
    while(j>=1)
    {
        if(dis[H[j]]>dis[t])
        {
            H[i]=H[j];
            pH[H[i]]=i;
            i=j;j=i>>1;
        }
        else break;
    }
    H[i]=t;pH[H[i]]=i;
}
void _shift_down(int i)
{
    int j=i<<1,t=H[i];
    while(j<tot)
    {
        if(j+1<tot&&dis[H[j]]>dis[H[j+1]])j++;
        if(dis[t]>dis[H[j]])
        {
            H[i]=H[j];
            pH[H[i]]=i;
            i=j;j=i<<1;
        }
        else break;
    }
    H[i]=t;pH[H[i]]=i;
}
void _prim()
{
    int i,j,k,t,ans=0;
    for(i=1;i<=N;i++)
    {
        H[i]=pH[i]=i;
        dis[i]=999999999;
    }
    dis[1]=0;tot=N+1;
    for(i=1;i<=N;i++)
    {
        k=H[1];H[1]=H[--tot];pH[H[1]]=1;pH[k]=0;
        _shift_down(1);
        if(dis[k]==999999999){printf("Keng Die!");return;}
        ans+=dis[k];
        for(j=last[k];j>0;j=next[j])
            if(pH[y[j]]!=0&&dis[y[j]]>z[j])
            {
                dis[y[j]]=z[j];
                _shift_up(pH[y[j]]);
            }
    }
    printf("%d",ans);
}
int main()
{
    int i,j,k,t;
    _getnum(N);_getnum(M);
    for(i=j=1;i<=M;i++)
    {
        _getnum(k);_getnum(t);_getnum(z[j]);
        y[j]=t;next[j]=last[k];last[k]=j++;
        y[j]=k;z[j]=z[j-1];next[j]=last[t];last[t]=j++;
    }
    _prim();
    return 0;
}

内容概要:本文详细介绍了基于Matlab实现的“梯级水光互补系统最大化可消纳电量期望短期优化调度模型”,属于电力系统领域高水平科研成果的复现(EI级别)。该模型聚焦于梯级水电站与光伏发电系统的协同优化调度,通过构建短期优化调度框架,旨在提升可再生能源的电量消纳能力并最大化系统综合效益。研究采用先进的数学优化方法对水光资源进行联合调度,充分考虑了光伏出力的不确定性、水资源约束、系统运行边界条件及电力平衡要求,实现了在多重约束下的电量期望最大化目标。模型不仅具备严谨的理论基础,还具有良好的工程应用前景,适用于新能源高比例渗透背景下电力系统的优化调度研究与实践。; 适合人群:具备电力系统分析、可再生能源利用或优化建模背景的研究生、科研人员及工程技术人员,特别适合致力于复现高水平学术论文(EI/顶刊)研究成果的学习者与开发者。; 使用场景及目标:① 学习并掌握梯级水电与光伏系统协同调度的建模思路与关键技术;② 熟悉基于Matlab的混合整数线性规划(MILP)或其他非线性优化方法在能源系统中的实际应用;③ 提升在新能源消纳、短期调度优化等方向的科研建模能力与代码实现水平,支持二次开发与创新研究。; 阅读建议:建议结合Matlab代码与优化理论同步研读,重点理解目标函数的设计逻辑、各类物理与运行约束的数学表达以及求解器的调用流程,推荐使用YALMIP等建模工具辅助实现,以提高模型构建效率与可读性,便于深入理解与后续拓展。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值