这道题也是我拿来测板子的
题意:给定一个r*c(r<=100,c<=100)的矩阵,每个元素表示该点高度,给定初始速度v,当从一点移动到另一点后,速度变为原先的2^(该点高度-移动到的点的高度),花费时间是1/原先高度,求从1,1到r,c的最短时间
题解:可以推出,每个点速度是一样的,这样就直接最短路了
#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
#include <cmath>
using namespace std ;
const int MAXN = 105 ;
const int MAXM = MAXN * MAXN * 20 ;
double INF = 9999999999.0 ;
double EPS = 0.000000001 ;
struct type
{
int end , next ;
double length ;
} ;
type edge[MAXM] ;
int val[MAXN][MAXN] ;
int head[MAXN*MAXN] ;
int cnt[MAXN*MAXN] ;
int V , R , C , Count ;
bool vis[MAXN*MAXN] ;
double dis[MAXN*MAXN] ;
double speed[MAXN][MAXN] ;
int dx[4] = { 0 , 0 , 1 , -1 } ;
int dy[4] = { 1 , -1 , 0 , 0 } ;
void addedge( int start , int end , double length )
{
edge[Count].end = end ;
edge[Count].length = length ;
edge[Count].next = head[start] ;
head[start] = Count ++ ;
}
double SPFA( int start , int n )
{
int i ;
memset( cnt , 0 , sizeof( cnt ) ) ;
memset( vis , 0 , sizeof( vis ) ) ;
for( i = 1 ; i <= n ; i ++ ) dis[i] = INF ;
dis[start] = 0 ;
queue < int > q ;
q.push( start ) ; vis[start] = true ; ++ cnt[start] ;
while( ! q.empty() )
{
int u , v ;
u = q.front() ; q.pop() ; vis[u] = false ;
for( i = head[u] ; i != -1 ; i = edge[i].next ) {
v = edge[i].end ;
if( dis[v] > dis[u] + edge[i].length + EPS )
{
dis[v] = dis[u] + edge[i].length ;
if( !vis[v] ){
q.push( v ) ;
vis[v] = true ;
if( ++cnt[v] > n ) return -1 ;
}
}
}
}
return 0 ;
}
int main()
{
int i , j , k ;
memset( head , -1 , sizeof( head ) ) ;
memset( vis , 0 , sizeof( vis ) ) ;
memset( speed , -1 , sizeof( speed ) ) ;
Count = 0 ;
scanf( "%d%d%d" , & V , & R , & C ) ;
for( i = 0 ; i < R ; i ++ )
for( j = 0 ; j < C ; j ++ )
scanf( "%d" , & val[i][j] ) ;
speed[0][0] = V ;
for( i = 1 ; i < C ; i ++ )
speed[0][i] = speed[0][i-1] * pow( 2 , double ( val[0][i-1] - val[0][i] ) ) ;
for( i = 1 ; i < R ; i ++ )
for( j = 0 ; j < C ; j ++ )
speed[i][j] = speed[i-1][j] * pow( 2 , double ( val[i-1][j] - val[i][j] ) ) ;
for( i = 0 ; i < R ; i ++ )
for( j = 0 ; j < C ; j ++ )
for( k = 0 ; k < 4 ; k ++ )
{
int tmpx = i + dx[k] ;
int tmpy = j + dy[k] ;
if( tmpx >= 0 && tmpx < R && tmpy >= 0 && tmpy < C )
addedge( ( i * C + j ) , ( tmpx * C + tmpy ) , 1.0 / speed[i][j] ) ;
}
SPFA( 0 , R * C ) ;
printf( "%.2f\n" , dis[R * C - 1] ) ;
return 0 ;
}
代码搓了点
分享到:
相关推荐
poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题...
poj 3414解题报告poj 3414解题报告poj 3414解题报告poj 3414解题报告
poj 1012解题报告poj 1012解题报告poj 1012解题报告poj 1012解题报告
poj 2329解题报告poj 2329解题报告poj 2329解题报告poj 2329解题报告
poj 1440解题报告 poj 1440解题报告 poj 1440解题报告 poj 1440解题报告
poj 3083解题报告poj 3083解题报告poj 3083解题报告poj 3083解题报告
poj 1659解题报告poj 1659解题报告poj 1659解题报告poj 1659解题报告
poj 3720解题报告poj 3720解题报告poj 3720解题报告poj 3720解题报告
poj1691解题报告 题目来源:http://acm.pku.edu.cn/JudgeOnline/showproblem?problem_id=1691(POJ No.1691) 解法: 搜索
poj2828解题报告,希望能帮到志同道合的算法爱好者
北大poj解题报告,希望能帮到软件工程的同学,每天一道,持之以恒,熟能生巧,与您共勉!
北大ACM1316解题报告
Problem 1061 青蛙的约会 poj解题报告。有源码。可直接提交C++实现
这是我发的第二篇解题报告,写的不怎么的。呵呵!!!!!
ACM Poj Pku 解题报告答案 打包 下载 600多题 史上最全 不是网上乱传的200多题,更不是100多题就挂着10分才能下的题 下了这个 大家也不要浪费分数去下载其它版本的了,基本上都有 共享 一起进步 中国加油 ACMer...
ACM POJ 解题报告北大POJ 大量解题代码
本人整理的POJ解题报告,一共有250道题
北大ACM在线评测系统POJ的题目解题报告。涵盖各种类型的acm题目。值得参考借鉴。打包下载。
不错的资源 收集了网上的一些解题报告 方便同学们学习参考
2遍dp poj_3613解题报告 poj_3613解题报告