关于最短路
2015-10-08 23:46
176 查看
把最近做的几道最短路一起总结一下吧
poj1062
中文题,不解释题意,对于渣渣来说,是先看别人题解后才做的,以酋长的允诺作为目标点,每一个物品与替换物之间会有一条边,由于交换有限制,所以建完图后,以每一个点为终点或是说最高点,先遍历一次,把不能交换或间接交换的点去掉,然后再进行最短路寻找。实际就是把所有可能都搞了一遍找最短
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
#include <queue>
#include <map>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <stack>
using namespace std;
const int inf=0x7fffffff;
int m,n;
int mapp[105][105],lev[105],x[105],dis[105];//分别是每两个物品之间的优惠,对应物品拥有者的等级,替换物品数,以i为目标点的距离
bool vis[105];
void dijkstra()
{
for(int i=1;i<=n;i++)dis[i]=mapp[0][i];//初始化每个点本身的值
for(int i=1;i<=n;i++)
{
int s=inf,node=0;
for(int j=1;j<=n;j++)
{
if(s>dis[j]&&!vis[j])
{
s=dis[j];
node=j;
}
}//找到值最小的那个点,以它为起点松弛
if(node==0)break;//找不到就说明已遍历完毕
vis[node]=1;
for(int j=1;j<=n;j++)
{
if(!vis[j]&&mapp[node][j]>0)
dis[j]=min(dis[j],mapp[node][j]+dis[node]);
}
}
}
int main()
{
cin>>m>>n;
memset(x,0,sizeof(x));
memset(mapp,0,sizeof(mapp));
memset(lev,0,sizeof(lev));
memset(dis,inf,sizeof(dis));
for(int i=1;i<=n;i++)
{
cin>>mapp[0][i]>>lev[i]>>x[i];
for(int j=1;j<=x[i];j++)
{
int v,w;
cin>>v>>w;
mapp[v][i]=w;
}
}
int maxn=inf;
for(int i=1;i<=n;i++)
{
int level=lev[i];
for(int j=1;j<=n;j++)
{
if(lev[j]>level||level-lev[j]>m)
vis[j]=true;
else vis[j]=false;
}
dijkstra();
if(dis[1]<maxn)maxn=dis[1];
}
/*for(int i=1;i<=n;i++)
cout<<dis[i]<<' ';
cout<<endl;
*/
cout<<maxn<<endl;
return 0;
}
还有另一种解法,dp
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
#include <queue>
#include <map>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <stack>
using namespace std;
int dp[105][105][105],lev[105],mapp[105][105];
int n,m;
int dfs(int x,int y,int z)
{
if(dp[x][y][z])return dp[x][y][z];
dp[x][y][z]=mapp[0][x];
for(int i=1;i<=n;i++)
{
if(lev[i]-y<=m&&z-lev[i]<=m&&mapp[i][x])
dp[x][y][z]=min(dp[x][y][z],mapp[i][x]+dfs(i,min(lev[i],y),max(lev[i],z)));
}
return dp[x][y][z];
}
int main()
{
memset(dp,0,sizeof(dp));
cin>>m>>n;
for(int i=1;i<=n;i++)
{
int x;
cin>>mapp[0][i]>>lev[i]>>x;
while(x--)
{
int v,w;
cin>>v>>w;
mapp[v][i]=w;
}
}
cout<<dfs(1,lev[1],lev[1])<<endl;
return 0;
}
poj2253
接下来几道题都是floyed
poj1125
poj2240
也可以跟poj1860那样去解
poj1062
中文题,不解释题意,对于渣渣来说,是先看别人题解后才做的,以酋长的允诺作为目标点,每一个物品与替换物之间会有一条边,由于交换有限制,所以建完图后,以每一个点为终点或是说最高点,先遍历一次,把不能交换或间接交换的点去掉,然后再进行最短路寻找。实际就是把所有可能都搞了一遍找最短
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
#include <queue>
#include <map>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <stack>
using namespace std;
const int inf=0x7fffffff;
int m,n;
int mapp[105][105],lev[105],x[105],dis[105];//分别是每两个物品之间的优惠,对应物品拥有者的等级,替换物品数,以i为目标点的距离
bool vis[105];
void dijkstra()
{
for(int i=1;i<=n;i++)dis[i]=mapp[0][i];//初始化每个点本身的值
for(int i=1;i<=n;i++)
{
int s=inf,node=0;
for(int j=1;j<=n;j++)
{
if(s>dis[j]&&!vis[j])
{
s=dis[j];
node=j;
}
}//找到值最小的那个点,以它为起点松弛
if(node==0)break;//找不到就说明已遍历完毕
vis[node]=1;
for(int j=1;j<=n;j++)
{
if(!vis[j]&&mapp[node][j]>0)
dis[j]=min(dis[j],mapp[node][j]+dis[node]);
}
}
}
int main()
{
cin>>m>>n;
memset(x,0,sizeof(x));
memset(mapp,0,sizeof(mapp));
memset(lev,0,sizeof(lev));
memset(dis,inf,sizeof(dis));
for(int i=1;i<=n;i++)
{
cin>>mapp[0][i]>>lev[i]>>x[i];
for(int j=1;j<=x[i];j++)
{
int v,w;
cin>>v>>w;
mapp[v][i]=w;
}
}
int maxn=inf;
for(int i=1;i<=n;i++)
{
int level=lev[i];
for(int j=1;j<=n;j++)
{
if(lev[j]>level||level-lev[j]>m)
vis[j]=true;
else vis[j]=false;
}
dijkstra();
if(dis[1]<maxn)maxn=dis[1];
}
/*for(int i=1;i<=n;i++)
cout<<dis[i]<<' ';
cout<<endl;
*/
cout<<maxn<<endl;
return 0;
}
还有另一种解法,dp
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
#include <queue>
#include <map>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <stack>
using namespace std;
int dp[105][105][105],lev[105],mapp[105][105];
int n,m;
int dfs(int x,int y,int z)
{
if(dp[x][y][z])return dp[x][y][z];
dp[x][y][z]=mapp[0][x];
for(int i=1;i<=n;i++)
{
if(lev[i]-y<=m&&z-lev[i]<=m&&mapp[i][x])
dp[x][y][z]=min(dp[x][y][z],mapp[i][x]+dfs(i,min(lev[i],y),max(lev[i],z)));
}
return dp[x][y][z];
}
int main()
{
memset(dp,0,sizeof(dp));
cin>>m>>n;
for(int i=1;i<=n;i++)
{
int x;
cin>>mapp[0][i]>>lev[i]>>x;
while(x--)
{
int v,w;
cin>>v>>w;
mapp[v][i]=w;
}
}
cout<<dfs(1,lev[1],lev[1])<<endl;
return 0;
}
poj2253
接下来几道题都是floyed
#include <iostream> #include <cstring> #include <string> #include <cmath> #include <algorithm> #include <queue> #include <map> #include <vector> #include <cstdio> #include <cstdlib> #include <stack> #include <iomanip> using namespace std; const int inf=0x7fffffff; struct pointt { int x; int y; }p[205]; double dis[205][205]; double dist(pointt a,pointt b) { return sqrt((double)((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y))); } int main() { int n,r=1; while(cin>>n&&n!=0) { for(int i=1;i<=n;i++) scanf("%d%d",&p[i].x,&p[i].y); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) { double d=dist(p[i],p[j]); dis[i][j]=dis[j][i]=d; } /*for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) cout<<dis[i][j]<<' '; cout<<endl; } */ for(int k=1;k<=n;k++) for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) { double d=max(dis[i][k],dis[k][j]); if(dis[i][j]-d>0.001)dis[i][j]=d; } /*for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) cout<<dis[i][j]<<' '; cout<<endl; } */ cout<<"Scenario #"<<r++<<'\n'<<"Frog Distance = "<<setiosflags(ios::fixed)<<setprecision(3)<<dis[1][2]<<endl; cout<<endl; } return 0; }
poj1125
#include <iostream> #include <cstring> #include <string> #include <cmath> #include <algorithm> #include <queue> #include <map> #include <vector> #include <cstdio> #include <cstdlib> #include <stack> #include <iomanip> using namespace std; int mapp[105][105]; int main() { int n; while(~scanf("%d",&n)) { if(n==0)break; memset(mapp,0x3f,sizeof(mapp)); for(int i=1;i<=n;i++) { int x,v,w; cin>>x; while(x--) { cin>>v>>w; mapp[i][v]=w; } } /* for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { cout<<mapp[i][j]<<' '; } cout<<endl; } */ for(int k=1;k<=n;k++) for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) mapp[i][j]=min(mapp[i][j],mapp[i][k]+mapp[k][j]); int maxlen=0x3f3f3f3f,sorce; for(int i=1;i<=n;i++) { int maxlen1=0; for(int j=1;j<=n;j++) { if(i!=j) maxlen1=max(maxlen1,mapp[i][j]); } if(maxlen1<maxlen) { maxlen=maxlen1; sorce=i; } } if(maxlen==0x3f3f3f3f)cout<<"disjoint\n"<<endl; else cout<<sorce<<' '<<maxlen<<endl; } return 0; }
poj2240
也可以跟poj1860那样去解
#include <iostream> #include <cstring> #include <string> #include <cmath> #include <algorithm> #include <queue> #include <map> #include <vector> #include <cstdio> #include <cstdlib> #include <stack> #include <iomanip> using namespace std; double dis[35][35]; const double eps=1e-12; map<string,int>mp; int main() { int n,m,r=1; double ans; while(cin>>n&&n) { memset(dis,0.0,sizeof(dis)); mp.clear(); for(int i=1;i<=n;i++) { string s; cin>>s; mp[s]=i; dis[i][i]=1.0; } cin>>m; while(m--) { string s1,s2; cin>>s1>>ans>>s2; dis[mp[s1]][mp[s2]]=ans; } for(int k=1;k<=n;k++) for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if (dis[i][j]-dis[i][k]*dis[k][j]<eps) dis[i][j]=dis[i][k]*dis[k][j]; bool flag=false; for(int i=1;i<=n;i++) { if(dis[i][i]-1.0>eps) { flag=true; break; } }printf("Case %d: ",r++); if(flag)printf("Yes\n"); else printf("No\n"); } return 0; }
相关文章推荐
- Java之——简单的网络爬虫实现
- base64的一个应用情景
- 个人数据备份的常见方案
- vim使用技巧
- HTTP协议请求过程
- Linux设置某软件开机自动启动的方法
- 《Web安全测试》读书笔记
- 用sphinx写文档
- Webpack初试
- [置顶] 【动态规划】0-1背包问题
- SNPassoc全基因组关联分析
- UVa 489 - Hangman Judge【字符串...】
- 马良建模学习笔记——软包制作十三
- jsp的转跳方法
- 《SQL注入攻击与防御》读书笔记
- Laravel数据库操作
- foj2204 Problem 2204 7 dp
- 用Node.js实现基于https的Restful风格webservice
- hdu 4421 Bit Magic 2-sat
- 设计模式之模板方法模式&&迪米特法则(代码Objective-C展示)