HDU 2489 Minimal Ratio Tree(DFS+prim)

1 篇文章 0 订阅

题意:在n个结点的完全图中,每条边和每个点都有一个权值。

定义:

从中选出m个结点组成一个连通图,顺序输出构成最小的ratio时所选取的m个点。

思路:边权之和最小,想到prim,同时n最大是15,直接暴力。先dfs出m个顶点,然后求出对应的ratio,找出最小的radio,输出结果就行了。

计算了一下15以内组合数的最大值是6000+,此时m为7或8,算上prim的O(n^2),时间是够的。一发AC!


#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<cmath>
#include<set>
#include<algorithm>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int N=18;
int n,m,g[N][N],node[N];
int arr[N];
double ans;
bool vis[N];
int Prim(){
    //lowcost[i]记录结果集到结点i的最短距离,mst[i]记录结果集中到i最短距离的端点
    int minn=inf,id=0,lowcost[N],mst[N];
    int res=0;
    for(int i=0;i<n;++i){
        if(vis[i]){
            id=i;
            break;
        }
    }
    for(int i=0;i<n;++i){
        if(vis[i]){
            lowcost[i]=g[id][i];
            mst[i]=id;
        }
    }
    for(int i=1;i<m;++i){
        minn=inf;
        for(int i=0;i<n;++i){
            if(vis[i]&&lowcost[i]&&lowcost[i]<minn){
                minn=lowcost[i];
                id=i;
            }
        }
        res+=minn;
        lowcost[id]=0;
        for(int i=0;i<n;++i){
            if(vis[i]&&lowcost[i]&&g[id][i]<lowcost[i]){
                lowcost[i]=g[id][i];
                mst[i]=id;
            }
        }
    }
    return res;
}
void dfs(int k,int num){//k是当前下标,num是已标记点数
    //printf("dfs:k=%d num=%d\n",k,num);
    if(num==m){
        int we=Prim();
        int wn=0;
        for(int i=0;i<n;++i){
            if(vis[i]){
                wn+=node[i];
            }
        }
        double res=(double)we/wn;
        if(res<ans){
            ans=res;
            for(int i=0,j=0;i<n;++i){
                if(vis[i]){
                    arr[j++]=i;//数组arr记录选取的结点
                }
            }
        }
    }else{
        if(k==n)return;//搜索终点
        if((n-k)<(m-num))return;//剪枝,减少递归层数
        vis[k]=1;
        dfs(k+1,num+1);
        vis[k]=0;
        dfs(k+1,num);
    }
}
int main(){
    while(~scanf("%d%d",&n,&m),m||n){
        for(int i=0;i<n;++i){
            scanf("%d",&node[i]);
        }
        for(int i=0;i<n;++i){
            for(int j=0;j<n;++j){
                scanf("%d",&g[i][j]);
            }
        }
        memset(vis,0,sizeof(vis));
        ans=99999999.0;
        dfs(0,0);
        for(int i=0;i<m-1;++i){
            printf("%d ",arr[i]+1);//题目下标从1开始
        }
        printf("%d\n",arr[m-1]+1);
    }
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值