2013数据结构课程设计之便利店选址(暴力枚举和随机函数两种做法)
[问题描述]
某小区决定在小区内部建一家便利店，现小区内部共有八栋楼，它们的地理坐标分别为：（10，20） （30，34） （19，25） （38，49.1） （9，38.1） （2，34） （5，8） （29，48）。同时，其中的住户人数分别为：30, 45, 28, 8, 36, 16, 78, 56。为了方便更多的住户购物，要求实现总体最优，请问便利店应该建立在哪里？
【提示】
1）便利店无论选址何处，八栋楼的居民均可直接到达，即八栋楼与便利店均相邻，且距离为直线距离；
2）八栋楼的居民人数为权重，应该方便大多数人，实现总体最优。
先是暴力找点的方法。
解题思路：自己开始想的就是暴力枚举，先找大范围，再找小范围。做这个题目就想到了的warmup2的1002题，但当时就是A不了。思路很简单，一步步地精确范围。先把整个地方划分成10*10的方格，再在里面找哪个最小，然后继续10*10每次都这样划分，精度确定跳出循环即可。详解见代码。
代码：
#include<iostream>
#include<cstring>
#include<string>
#include<cmath>
#include<cstdio>
#include<algorithm>
#include<ctime>
using namespace std;
int n; //n栋楼
double minx,miny,maxx,maxy,px,py; //找到四个边界
double ansx,ansy,sum; //最后的结果
struct mq
{
double x;
double y;
int peo;
};
mq node[1005];
double cal(double px,double py) //计算值
{
int i;
double sum=0;
for(i=0;i<n;i++)
sum+=sqrt((px-node[i].x)*(px-node[i].x)+(py-node[i].y)*(py-node[i].y))*node[i].peo;
return sum;
}
void solve()
{
double i,j;
double fenx=maxx-minx; //把x分成若干份
double feny=maxy-miny; //把y分成若干份
while(fenx>0.01&&feny>0.01) //暴力找点
{
fenx/=10.0,feny/=10.0;
for(i=minx;i<=maxx;i+=fenx)
for(j=miny;j<=maxy;j+=feny)
{
double tmp=cal(i,j);
if(tmp<sum)
{
sum=tmp;
ansx=i;
ansy=j;
}
}
minx=ansx-fenx;
miny=ansy-feny;
maxx=ansx+fenx;
maxy=ansy+feny;
}
}
int main()
{
int i;
//freopen("input.txt","r",stdin);
//freopen("output.txt","w",stdout);
while(~scanf("%d",&n)) //n栋楼
{
scanf("%lf%lf%d",&node[0].x,&node[0].y,&node[0].peo);
minx=maxx=node[0].x; miny=maxy=node[0].y;
//找到四个边界
for(i=1;i<n;i++)
{
scanf("%lf%lf%d",&node[i].x,&node[i].y,&node[i].peo);
if(node[i].x<minx) minx=node[i].x;
if(node[i].x>maxx) maxx=node[i].x;
if(node[i].y<miny) miny=node[i].y;
if(node[i].y>maxy) maxy=node[i].y;
}
sum=100000000;
solve();
cout<<"便利店选址坐标为:"<<endl;
cout<<"x: "<<ansx<<" "<<"y: "<<ansy<<endl;
cout<<"最优解为: "<<sum<<endl;
}
return 0;
}
/*
8
10 20 30
30 34 45
19 25 28
38 49.1 8
9 38.1 36
2 34 16
5 8 78
29 48 56
便利店选址坐标为:
x: 16.5404 y: 27.4362
最优解为: 5146.85
*/
然后就是LCM说的随机算法，也没想到那一块，被他一说是有些道理，但是点的范围一扩大的话，就没那么准确了。
随机算法，我是直接采用的随机数然后%lenth，区间宽度，相当于0~lenth-1。随机函数对范围要求比较严格吧，范围如果太大，随机也没用了。具体实现见代码：
代码：
#include<iostream>
#include<cstring>
#include<string>
#include<cmath>
#include<cstdio>
#include<algorithm>
#include<ctime>
using namespace std;
int n; //n栋楼
struct mq
{
double x;
double y;
int peo;
};
mq node[1005];
double cal(double px,double py) //计算值
{
int i;
double sum=0;
for(i=0;i<n;i++)
sum+=sqrt((px-node[i].x)*(px-node[i].x)+(py-node[i].y)*(py-node[i].y))*node[i].peo;
return sum;
}
int main()
{
int i;
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
double minx,miny,maxx,maxy,px,py; //找到四个边界
double ansx,ansy; //最后的结果
while(~scanf("%d",&n)) //n栋楼
{
scanf("%lf%lf%d",&node[0].x,&node[0].y,&node[0].peo);
minx=maxx=node[0].x; miny=maxy=node[0].y;
//找到四个边界
for(i=1;i<n;i++)
{
scanf("%lf%lf%d",&node[i].x,&node[i].y,&node[i].peo);
if(node[i].x<minx) minx=node[i].x;
if(node[i].x>maxx) maxx=node[i].x;
if(node[i].y<miny) miny=node[i].y;
if(node[i].y>maxy) maxy=node[i].y;
}
minx*=100,maxx*=100,miny*=100,maxy*=100; //边界扩大一百倍
//找到边界了就可以随机了
int lenx,leny;
lenx=ceil(maxx-minx),leny=ceil(maxy-miny);
double sum=1000000000;
srand((unsigned)time(NULL)); //播种
for(i=1;i<=500000;i++) //随机50W次
{
px=rand()%lenx+minx;
py=rand()%leny+miny;
px/=100.0;
py/=100.0;
double tmp=cal(px,py);
if(tmp<sum)
{
sum=tmp;
ansx=px;
ansy=py;
}
}
cout<<"便利店选址坐标为:"<<endl;
cout<<"x: "<<ansx<<" "<<"y: "<<ansy<<endl;
cout<<"最优解为: "<<sum<<endl;
}
return 0;
}
/*
8
10 20 30
30 34 45
19 25 28
38 49.1 8
9 38.1 36
2 34 16
5 8 78
29 48 56
便利店选址坐标为:
x: 16.56 y: 27.44
最优解为: 5146.85
*/