1295: [SCOI2009]最长距离

本文介绍了一种基于暴力SPFA算法解决特定最短路径问题的方法。通过遍历地图上的每个点作为起点,利用SPFA算法计算到其他点的最短距离,并在此过程中考虑障碍物的影响。最终目的是找到所有距离小于给定阈值T的点中,与起点的最大欧几里得距离。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

暴力SPFA。

若该点有一个障碍物,则经过该点后距离+1,对于所有距当前起点的距离小于T的点,最大距离(欧几里得距离)。

即得答案。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
using namespace std;
const int dx[4]={0,1,0,-1};
const int dy[4]={1,0,-1,0};
int mp[35][35],d[35][35],T,n,m;
bool inq[35][35];
double ans;
void spfa(int s,int t){
	memset(d,0x3f,sizeof(d));
	queue<int>qx,qy;qx.push(s);qy.push(t);
	d[s][t]=mp[s][t];
	while(!qx.empty()){
		int u=qx.front(),v=qy.front();qx.pop();qy.pop();
		inq[u][v]=false;
		for(int i=0;i<4;i++){
			int tx=u+dx[i],ty=v+dy[i];
			if(tx<1||tx>n||ty<1||ty>m)continue;
			if(d[tx][ty]>d[u][v]+mp[tx][ty]){
				d[tx][ty]=d[u][v]+mp[tx][ty];
				if(!inq[tx][ty]){qx.push(tx);qy.push(ty);inq[tx][ty]=true;}
			}
		}
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
	if(d[i][j]<=T)
	ans=max(ans,sqrt((i-s)*(i-s)+(j-t)*(j-t)));
}
char ch[50];
int main(){
	scanf("%d%d%d",&n,&m,&T);
	for(int i=1;i<=n;i++){
		scanf("%s",ch+1);
		for(int j=1;j<=m;j++)
		mp[i][j]=ch[j]-'0';
	}
	ans=0.0;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
	spfa(i,j);
	printf("%.6lf",ans);
	return 0;
}


### 关于 SCOI2009 WINDY 数的解法 #### 定义与问题描述 WINDY数是指对于任意两个相邻位置上的数字,它们之间的差至少为\(2\)。给定正整数区间\([L, R]\),计算该范围内有多少个WINDY数。 #### 动态规划方法解析 为了高效解决这个问题,可以采用动态规划的方法来处理。定义状态`dp[i][j]`表示长度为`i`且高位是`j`的WINDY数的数量[^3]。 - **初始化** 对于单个数字的情况(即只有一位),显然每一位都可以单独构成一个合法的WINDY数,因此有: ```cpp dp[1][d] = 1; // d ∈ {0, 1,...,9} ``` - **状态转移方程** 当考虑多位数时,如果当前位选择了某个特定数值,则下一位的选择会受到限制——它必须满足与前一位相差不小于2的要求。具体来说就是当上一高位为`pre`时,当前位置可选范围取决于`pre`的具体取值: - 如果`pre >= 2`, 则可以选择`{0... pre-2}` - 否则只能从剩余的有效集合中选取 这样就可以通过遍历所有可能的状态来进行状态间的转换并累加结果。 - **边界条件处理** 特殊情况下需要注意的是,在实际应用过程中还需要考虑到给出区间的上下限约束。可以通过逐位比较的方式判断是否越界,并据此调整有效状态空间大小。 ```cpp // 计算不超过num的大windy数数量 int calc(int num){ int f[15], g[15]; memset(f, 0, sizeof(f)); string s = to_string(num); n = s.size(); for (char c : s) { a[++len] = c - '0'; } // 初始化f数组 for (int i=0;i<=9;++i)f[1][i]=1; // DP过程省略... return sum; } long long solve(long long L,long long R){ return calc(R)-calc(L-1); } ``` 此代码片段展示了如何利用预处理好的`dp`表快速查询指定范围内的WINDY数总量。其中`solve()`函数用于返回闭区间\[L,R\]内符合条件的总数;而辅助函数`calc()`负责根据传入参数构建相应的状态序列并终得出答案。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值