POJ1258

计算机科学与技术系上机实验报告
课程名称：算法设计与分析班级：计科081 实验日期：2010-11-13
姓名：学号：指导教师：
实验序号：3 实验成绩：
一、实验名称
贪心算法求农场网络
二、实验目的及要求
1、学会使用在线测评的算法题目评分系统；
2、通过应用问题，加深对贪心算法的理解；
三、实验环境
C、C++编写调试工具，北京大学ICPC在线测评系统POJ
四、实验内容
2、阅读题目，建立程序的流程图；
3、使用贪心算法，实现本题；
4、进行简单测试，完成之后提交到POJ系统。

五、算法描述及实验步骤
贪心算法原理：
用贪心算法，使用当前最优的解决方案设计程序，最终得到的算法也是最优解，或者相对最优解
农场问题描述：
描述
农夫约翰被选举为他所在城市的市长，他的一个竞选诺言就是将这块地方所有的农场安装上互联网。

当然，他需要你的帮助。

约翰给他自己的农场已经安装了一个高速的互联网连接，他准备让其它农场共享这个连接。

为了花费最小，他想铺设最少的光纤到其它农场。

你得到了连接每两个农场需要花费的光纤长度的列表，你必须找到连接所有农场的最小花费。

两个农场的距离值不会超过100，000。

输入
输入包含几组例子。

每组例子的第一行是农场数量N (3 <= N <= 100)，接着的N行是一个NxN的邻接矩阵，记录了每两个农场之间的距离的整数值。

输出
对于每个案例，输出一个行，这个一行有一个整数，就算连接全部农场需要的最小的光纤长度。

例子输入
4
0 4 9 21
4 0 8 17
9 8 0 16
21 17 16 0
例子输出
28
此办法应该归入到kruskal算法。

函数体里面的v是最主要的程序体，主要用来进行最小边的查找和插入。

1、输入数据。

int a[101][101];//定义数组矩阵存储农场之间的距离。

int b[101];//定义数组b存储已被连接在子树里面的节点序号
cout<<"请输入邻接矩阵"<<endl;
for(i=1;i<=n;i++)
{for(j=1;j<=n;j++)
cin>>a[i][j];
}//输入矩阵
2、请给出算法框架！（流程图、纯文字语言、有注释的关键C代码均可）
v:
x=100002;p=0;q=0;
for(i=1;i<=n;i++)
{
for(j=1;j<i;j++)
if(x>a[i][j])/*因为俩个农场间连接的距离相同所以只用计算矩阵的上半部分即能得到正确解*/
{x=a[i][j];
c=i,l=j;}
}
for (w=1;w<=k;w++)
{if (b[w]==c)
p=c;
if (b[w]==l)
q=l;}//判断最短边的端点是否已在生成树中
if (p==c&&q==l)
a[c][l]=100001;/*置最短边的长度为100001 下次查询的时候就不会选到此边*/
else
{ s=s+x;
a[c][l]=100001; /*置最短边的长度为100001 下次查询的时候就不会选到此边*/
if(p==c&&q!=l)
{k=k+1;
b[k]=l;}
if(p!=c&&q==l)
{k=k+1;
b[k]=c;}
if(p!=c&&q!=l)
{k=k+1;
b[k]=c;
k=k+1;
b[k]=l;}
}//如果不在则插入数组
//cout<<"k"<<k<<endl;
if(k<n)
goto v;
3、输出
cout<<s<<endl;
六、调试过程及实验结果
详细记录程序在调试过程中出现的问题及解决方法。

进入程序执行，选择输入的数组的行列运行。

出现如下面输出的正确结果：
七、总结
对上机实践结果进行分析，问题回答，上机的心得体会及改进意见.
题目意思，输入距离的矩阵，存放在多维数组里。

用贪心算法求连接各个农场所需的最短距离。

由于对算法不太了解，和编程基础、习惯不好导致程序冗余度大，并且在vc++6.0上可以运行但在系统上不能正确运行。

以后再编程的过程中还要努力和改进。

return type of 'main' should be 'int' instead of 'void'
这是在北大在线测评系统上出现的错误，还需改进。

八、附录
源程序（核心代码）清单或使用说明书，可另附纸:
#include"iostream"
using namespace std;
#include"stdio.h"
int a[101][101];//定义矩阵
int b[101];
void main()
{int n,s=0,i,j,k=0,l,c,w,p,q,x;
cout<<"请输入农场数量n(3<=n<=100)"<<endl;
cin>>n;// 定义农场数量即矩阵规模
cout<<"请输入邻接矩阵"<<endl;
for(i=1;i<=n;i++)
{for(j=1;j<=n;j++)
cin>>a[i][j];
}//输入矩阵
v:
x=100002;p=0;q=0;
for(i=1;i<=n;i++)
{
for(j=1;j<i;j++)
if(x>a[i][j])
{x=a[i][j];
c=i,l=j;}
}
for (w=1;w<=k;w++)
{if (b[w]==c)
p=c;
if (b[w]==l)
q=l;}//判断最短边的端点是否已在生成树中if (p==c&&q==l)
a[c][l]=100001;
else
{ s=s+x;
a[c][l]=100001;
if(p==c&&q!=l)
{k=k+1;
b[k]=l;}
if(p!=c&&q==l)
{k=k+1;
b[k]=c;}
if(p!=c&&q!=l)
{k=k+1;
b[k]=c;
k=k+1;
b[k]=l;}
}//如果不在则插入数组if(k<n)
goto v;
cout<<s<<endl;}。