GraphTheory

poj-1751Highways(最小生成树)

2019-07-29  本文已影响1人  雨落八千里

题意:

给了n个点的坐标,给了m条已经相连的边。输出建成最小生成树的其他边。

思路:

简单,普通的建立,之所以要写这篇博客,是因为它动不动超内存。为了以后不踩坑,于是写篇博客警醒自己。

克鲁斯卡尔

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int M=600000;//之前一直以为是这个大了才超的内存,但是参考了一下大神的博客,发现不是这个的问题
struct node
{
    int a,b;
    double w;
} edge[M];
struct tree
{
    int a,b;
} point[800];//之前是用两数组存,但是超内存,发现是这个两个数组的问题
int n,m;
int root[800];
int find(int x)
{
    if(x==root[x])
    {
        return x;
    }
    else
    {
        return root[x]=find(root[x]);
    }
}
int ouin(int a,int b)
{
    int x=find(a);
    int y=find(b);
    if(x==y)
    {
        return 0;
    }
    else
    {
        root[x]=y;
        return 1;
    }

}
bool cmp(node a,node b)
{
    return a.w<b.w;
}
int main( )
{
    int a,b;
    scanf("%d",&n);
    for(int i=1; i<=n; i++)
    {
        root[i]=i;
        scanf("%d%d",&point[i].a,&point[i].b);
    }
    int cnt=0;
    for(int i=1; i<=n; i++)
    {
        for(int j=i+1; j<=n; j++)//我以前都是j=1开始,但是这次超时,逼成这样,以前知道这样可以,但是偷懒不想
        {
            edge[cnt].a=i;
            edge[cnt].b=j;
            edge[cnt++].w=(point[i].a-point[j].a)*(point[i].a-point[j].a)+(point[i].b-point[j].b)*(point[i].b-point[j].b);
        }
    }
    scanf("%d",&m);
    for(int i=1; i<=m; i++)
    {
        scanf("%d%d",&a,&b);
        edge[cnt].a=a;
        edge[cnt].b=b;
        edge[cnt++].w=0;
    }
    sort(edge,edge+cnt,cmp);
    int res=0;
    int ans=0;
    for(int i=0; i<cnt; i++)
    {
        if(ouin(edge[i].a,edge[i].b))
        {
            if(edge[i].w!=0)
            {
               cout<<edge[i].a<<" "<<edge[i].b<<endl;
            }
            res++;
        }
        if(res==n-1)
        {
            break;
        }
    }
    return 0;
}

Prim

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define INF 0x3f3f3f3f
using namespace std;
int mp[800][800];
struct node
{
    int a,b;
}point[800];
int vis[800];
int pre[800];//记录每个节点的父节点
int low[800];
void prim(int cost[ ][800],int n)
{
    memset(vis,0,sizeof(vis));
    int i,j;
    int pos=1;//父节点从1开始
    for(i=1;i<=n;i++)
    {
        low[i]=mp[i][1];
        pre[i]=pos;
    }
    vis[1]=1;
    for(int i=2;i<=n;i++)
    {
        int mi=INF;
        for(int j=1;j<=n;j++)
        {
            if(!vis[j]&&low[j]<mi)
            {
                mi=low[j];
                pos=j;
            }
        }
        if(mp[pre[pos]][pos])
        {
            cout<<pre[pos]<<" "<<pos<<endl;
        }
        vis[pos]=1;
        for(int j=1;j<=n;j++)
        {
            if(!vis[j]&&low[j]>cost[pos][j])
            {
                low[j]=cost[pos][j];//如果通过当前的节点pos到达节点j更近,意味着到达节点j必须经过节点pos,所以节点j的父节点就是pos
                pre[j]=pos;
            }
        }
    }
}
int add(int i,int j)
{
    return (point[i].a-point[j].a)*(point[i].a-point[j].a)+(point[i].b-point[j].b)*(point[i].b-point[j].b);

}
int main( )
{
    int n,m;
    int x,y;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&point[i].a,&point[i].b);
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            mp[i][j]=add(i,j);
        }
    }
    cin>>m;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        mp[x][y]=mp[y][x]=0;
    }
    prim(mp,n);
    return 0;
}
上一篇下一篇

猜你喜欢

热点阅读