A. HDU1142 Walk Through the Forest
单源最短路径+记忆化搜索。
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#define INF 10000000
using namespace std;
int G[1005][1005],dis[1005],vis[1005];
int dp[1005];
int n;
void dijk(int s)
{
int i,j;
for(i=1; i<=n; i++)
{
vis[i]=0;
dis[i]=G[s][i];
}
vis[s]=1;
for(i=1; i<=n; i++)
{
int MIN=INF,f=0;
for(j=1; j<=n; j++) if(!vis[j])
if(dis[j]<MIN) MIN=dis[j],f=j;
vis[f]=1;
for(j=1; j<=n; j++) if(!vis[j])
dis[j]=min(dis[j],MIN+G[f][j]);
}
return ;
}
int dfs(int x)
{
if(dp[x]!=-1) return dp[x];
dp[x]=0;
for(int i=1;i<=n;i++)
{
if(G[x][i]<INF&&dis[x]>dis[i])
dp[x]+=dfs(i);
}
return dp[x];
}
int main()
{
int m,s,d,x,y,v,i,j;
while(~scanf("%d",&n)&&n)
{
scanf("%d",&m);
for(i=0;i<=n;i++)
for(j=0;j<=n;j++)
G[i][j]=i==j?0:INF;
for(i=0;i<m;i++)
{
scanf("%d%d%d",&x,&y,&v);
G[y][x]=G[x][y]=v;
}
dijk(2);
memset(dp,-1,sizeof(dp));
dp[2]=1;
dfs(1);
printf("%d\n",dp[1]);
}
return 0;
}
B.一个人的旅行 HDU2066
单源最短路径,对于他家相连的城市直接把距离设为0即可。注意有重边。
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#define INF 10000000
using namespace std;
int G[1005][1005],dis[1005],vis[1005];
int n=1005;
void dijk(int s)
{
int i,j;
for(i=0; i<n; i++)
{
vis[i]=0;
dis[i]=G[s][i];
}
vis[0]=1;
for(i=0; i<n; i++)
{
int MIN=INF,f=0;
for(j=0; j<n; j++) if(!vis[j])
if(dis[j]<MIN) MIN=dis[j],f=j;
vis[f]=1;
for(j=0; j<n; j++) if(!vis[j])
dis[j]=min(dis[j],MIN+G[f][j]);
}
return ;
}
int main()
{
int m,s,d,x,y,v,i,j;
while(~scanf("%d%d%d",&m,&s,&d))
{
int ans=INF;
for(i=0; i<n; i++)
for(j=0; j<n; j++)
G[i][j]= i==j?0:INF;
for(i=0; i<m; i++)
{
scanf("%d%d%d",&x,&y,&v);
G[y][x]=G[x][y]=min(v,G[x][y]);
}
for(i=0; i<s; i++)
{
scanf("%d",&x);
G[0][x]=0;
}
dijk(0);
for(i=0;i<d;i++)
{
scanf("%d",&x);
ans=min(ans,dis[x]);
}
printf("%d\n",ans);
}
return 0;
}
C - Minimum Transport Cost HDU1385
Floyd多源最短路径,而且要按照路径输出。
其实Floyd最短路是基于DP思想的。首先要理解它的状态表示方式和状态转移方程。
而且路径输出要求在距离相同时输出字典序最小的路径,所以要处理一下。
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#define INF 10000000
using namespace std;
int G[1005][1005];
int pre[1005][1005];
int dis[1005][1005];
int tax[1005];
int n;
void floyd()
{
int i,j,k;
for (i=0; i<n; i++)
for (j=0; j<n; j++)
dis[i][j]=G[i][j],pre[i][j]=j;
for (k=0; k<n; k++)
for (i=0; i<n; i++)
for (j=0; j<n; j++)
{
if(dis[i][k]>=INF||dis[k][j]>=INF) continue;
int t=dis[i][k]+dis[k][j]+tax[k];
if(t<dis[i][j])
{
dis[i][j]=t;
pre[i][j]=pre[i][k];
}
else if(t==dis[i][j])
{
pre[i][j]=min(pre[i][j],pre[i][k]);
}
}
}
int main()
{
int d,i,j;
while(scanf("%d",&n)!=EOF&&n)
{
for(i=0; i<n; i++)
for(j=0; j<n; j++)
{
scanf("%d",&d);
if(d==-1) G[i][j]=INF;
else G[i][j]=d;
}
for(i=0;i<n;i++)
scanf("%d",&tax[i]);
floyd();
int a,b;
while(scanf("%d%d",&a,&b)&&a+b>0)
{
int t=a-1;
printf("From %d to %d :\n",a,b);
printf("Path: %d",a);
while(t+1!=b)
{
t=pre[t][b-1];
printf("-->%d",t+1);
}
printf("\nTotal cost : %d\n\n",dis[a-1][b-1]);
}
}
return 0;
}
D - Minimal Ratio Tree HDU2489
给你n个点 (n<15),从中取m个,组成一棵树,并且使得 (所有边权值之和 / 所有点权值之和)最小。
n是非常小的,所以可以用深搜枚举所有的可能的组合。
在所有组合中取最小值。
#include <algorithm>
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#define INF 0x7fffffff
using namespace std;
int n,m;
bool choose[20];
bool vis[20];
int dis[20];
int w[20],G[20][20];
int q[20];
double ans;
int prim(int s)
{
int i,j;
int sum=0;
for(i=0; i<n; i++)
{
dis[i]=G[i][s];
if(choose[i]) vis[i]=0;
else vis[i]=1;
}
vis[s]=1;
for(i=1; i<m; i++)
{
int MIN=INF,f=0;
for(j=0; j<n; j++)
if(!vis[j]&&dis[j]<MIN) f=j,MIN=dis[j];
vis[f]=1;
sum+=MIN;
for(j=0; j<n; j++) if(!vis[j])
dis[j]=min(dis[j],G[f][j]);
}
return sum;
}
void dfs(int last,int t)
{
if(t==m)
{
double sum_node=0;
for(int i=0; i<n; i++) if(choose[i])
sum_node+=w[i];
int sum_edge=prim(last);
double rat=sum_edge*1.0/sum_node;
if(ans>rat)
{
ans=rat;
for(int k=0,i=0; i<n; i++) if(choose[i])
q[k++]=i;
}
}
for(int i=last; i<n; i++) if(!choose[i])
{
choose[i]=1;
dfs(i,t+1);
choose[i]=0;
}
}
int main()
{
int x,i,j;
while(scanf("%d%d",&n,&m)!=EOF&&m+n)
{
ans=1e10;
memset(choose,0,sizeof(choose));
for(i=0; i<n; i++)
scanf("%d",&w[i]);
for(i=0; i<n; i++)
for(j=0; j<n; j++)
{
scanf("%d",&x);
G[i][j]=x;
}
dfs(0,0);
for(i=0; i<m-1; i++)
printf("%d ",q[i]+1);
printf("%d\n",q[i]+1);
}
return 0;
}
E - Slim Span
题目给出n个点,m条边,要求组成一棵树并且使得 最大边与最小边尽量接近。
很容易想到用kruskal,在最短边确定的情况下,依据贪心的方法,所得的最小生成树的最大边也是尽量小的。
排序后取最小的一条边作为最短边,进行一次kruskal。
依次枚举每条边作为最短边,进行kurskal,最后取最大边与最小边最接近的即可。
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#define INF 0x7fffffff
using namespace std;
int a[1005];
struct node
{
int x,y,w;
}e[100004];
int cmp(node a ,node b)
{
return a.w<b.w;
}
int Find(int x)
{
if(x!=a[x]) return a[x]=Find(a[x]);
else return x;
}
int main()
{
int n,m,i,j;
while(scanf("%d%d",&n,&m)!=EOF&&m+n)
{
for(i=0;i<m;i++)
scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].w);
sort(e,e+m,cmp);
int ans=INF;
for(j=0;j<m;j++)
{
int cnt=0,L=INF,R=0;
for(i=0;i<=n;i++)
a[i]=i;
for(i=j;i<m;i++)
{
int x=Find(e[i].x),y=Find(e[i].y);
if(x!=y)
{
a[x]=y;
cnt++;
L=min(L,e[i].w);
R=max(R,e[i].w);
}
if(cnt==n-1) break;
}
if(cnt==n-1) ans=min(ans,R-L);
}
if(ans!=INF) printf("%d\n",ans);
else printf("-1\n");
}
return 0;
}
F - Highways poj1751
最小生成树,但是要输出在生成树过程中添加的边。
用prim算法,只要在其中加一个数组记录当前点与生成树中距离最近的点。 比如 dis_p[2]=1, 表示 树外一点 2 , 距离最近的树内的点为 1.
所以每次更新dis的时候同时更新dis_p即可。
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#define INF 0x7fffffff
using namespace std;
int x[1005];
int y[1005];
bool vis[1005];
int dis[1005];
int dis_p[1005];
int G[1005][1005];
struct node
{
int x,y;
int dis;
}ans[1005];
int getdis(int a,int b)
{
return (x[a]-x[b])*(x[a]-x[b])+(y[a]-y[b])*(y[a]-y[b]);
}
int main()
{
int n,m,i,j;
while(~scanf("%d",&n))
{
for(i=0;i<n;i++)
scanf("%d%d",&x[i],&y[i]);
for(i=0;i<n;i++)
for(j=i;j<n;j++)
{
if(j==i) G[i][j]=0;
else G[i][j]=G[j][i]=getdis(i,j);
}
scanf("%d",&m);
for(i=0;i<m;i++)
{
int a,b;
scanf("%d%d",&a,&b);
G[a-1][b-1]=G[b-1][a-1]=0;
}
memset(vis,0,sizeof(vis));
for(i=0;i<n;i++) dis[i]=G[i][0],dis_p[i]=0;
vis[0]=1;
for(i=1;i<n;i++)
{
int MIN=INF,f=0;
for(j=0;j<n;j++)
if(!vis[j]&&dis[j]<MIN) MIN=dis[j],f=j;
vis[f]=1;
if(MIN!=0) printf("%d %d\n",dis_p[f]+1,f+1);
for(j=0;j<n;j++) if(!vis[j])
if(dis[j]>G[f][j]) dis[j]=G[f][j],dis_p[j]=f;
}
}
return 0;
}