【动态规划】【图论】【C++算法】1575统计所有可行路径

作者推荐

【动态规划】【字符串】【行程码】1531. 压缩字符串

本文涉及知识点

动态规划汇总
图论

LeetCode1575统计所有可行路径

给你一个 互不相同 的整数数组，其中 locations[i] 表示第 i 个城市的位置。同时给你 start，finish 和 fuel 分别表示出发城市、目的地城市和你初始拥有的汽油总量每一步中，如果你在城市 i ，你可以选择任意一个城市 j ，满足 j != i 且 0 <= j < locations.length ，并移动到城市 j 。从城市 i 移动到 j 消耗的汽油量为 |locations[i] - locations[j]|，|x| 表示 x 的绝对值。
请注意， fuel 任何时刻都 不能 为负，且你 可以 经过任意城市超过一次（包括 start 和 finish ）。
请你返回从 start 到 finish 所有可能路径的数目。
由于答案可能很大， 请将它对 10^9 + 7 取余后返回。
示例 1：
输入：locations = [2,3,6,8,4], start = 1, finish = 3, fuel = 5
输出：4
解释：以下为所有可能路径，每一条都用了 5 单位的汽油：
1 -> 3
1 -> 2 -> 3
1 -> 4 -> 3
1 -> 4 -> 2 -> 3
示例 2：
输入：locations = [4,3,1], start = 1, finish = 0, fuel = 6
输出：5
解释：以下为所有可能的路径：
1 -> 0，使用汽油量为 fuel = 1
1 -> 2 -> 0，使用汽油量为 fuel = 5
1 -> 2 -> 1 -> 0，使用汽油量为 fuel = 5
1 -> 0 -> 1 -> 0，使用汽油量为 fuel = 3
1 -> 0 -> 1 -> 0 -> 1 -> 0，使用汽油量为 fuel = 5
示例 3：
输入：locations = [5,2,1], start = 0, finish = 2, fuel = 3
输出：0
解释：没有办法只用 3 单位的汽油从 0 到达 2 。因为最短路径需要 4 单位的汽油。
提示：
2 <= locations.length <= 100
1 <= locations[i] <= 109
所有 locations 中的整数 互不相同 。
0 <= start, finish < locations.length
1 <= fuel <= 200

动态规划

令n = locations.length

动态规划的状态表示

dp[l][f] 表示 使用f单位的汽油，终点是l的可行路径数。状态数共：nm，故空间复杂度:O(nm)。

动态规划的转移方程

通过前置条件转移后置条件。
枚举符合以下条件的l1：
l1!=l
f1 = |loc[l1]-[l]| + f <= fuel。
dp[l1][[f1] += dp[l][f]
时间复杂度: O(nnm)

动态规划的填表顺序

f从小到大，确保动态规划的无后效性。

动态规划的初始状态

dp[start][0] =1 ，其它全部为0。

动态规划的返回值

dp[finish]之和

代码

核心代码

template<int MOD = 1000000007>
class C1097Int
{
public:
	C1097Int(long long llData = 0) :m_iData(llData% MOD)
	{

	}
	C1097Int  operator+(const C1097Int& o)const
	{
		return C1097Int(((long long)m_iData + o.m_iData) % MOD);
	}
	C1097Int& operator+=(const C1097Int& o)
	{
		m_iData = ((long long)m_iData + o.m_iData) % MOD;
		return *this;
	}
	C1097Int& operator-=(const C1097Int& o)
	{
		m_iData = (m_iData + MOD - o.m_iData) % MOD;
		return *this;
	}
	C1097Int  operator-(const C1097Int& o)
	{
		return C1097Int((m_iData + MOD - o.m_iData) % MOD);
	}
	C1097Int  operator*(const C1097Int& o)const
	{
		return((long long)m_iData * o.m_iData) % MOD;
	}
	C1097Int& operator*=(const C1097Int& o)
	{
		m_iData = ((long long)m_iData * o.m_iData) % MOD;
		return *this;
	}
	bool operator<(const C1097Int& o)const
	{
		return m_iData < o.m_iData;
	}
	C1097Int pow(long long n)const
	{
		C1097Int iRet = 1, iCur = *this;
		while (n)
		{
			if (n & 1)
			{
				iRet *= iCur;
			}
			iCur *= iCur;
			n >>= 1;
		}
		return iRet;
	}
	C1097Int PowNegative1()const
	{
		return pow(MOD - 2);
	}
	int ToInt()const
	{
		return m_iData;
	}
private:
	int m_iData = 0;;
};

class Solution {
public:
	int countRoutes(vector<int>& locations, int start, int finish, int fuel) {
		const int n = locations.size();
		vector<vector<C1097Int<> >> dp(n, vector<C1097Int<>>(fuel + 1));
		dp[start][0] = 1;
		for (int f = 0; f < fuel; f++)
		{
			for (int l = 0; l < n; l++)
			{
				for (int l1 = 0; l1 < n; l1++)
				{
					if (l1 == l)
					{
						continue;
					}
					int f1 = f + abs(locations[l1] - locations[l]);
					if (f1 <= fuel)
					{
						dp[l1][f1] += dp[l][f];
					}
				}
			}
		}
		return std::accumulate(dp[finish].begin(), dp[finish].end(), C1097Int()).ToInt();
	}
};

测试用例

template<class T>
void Assert(const T& t1, const T& t2)
{
	assert(t1 == t2);
}

template<class T>
void Assert(const vector<T>& v1, const vector<T>& v2)
{
	if (v1.size() != v2.size())
	{
		assert(false);
		return;
	}
	for (int i = 0; i < v1.size(); i++)
	{
		Assert(v1[i], v2[i]);
	}

}

int main()
{	
	vector<int> locations;
	int start,  finish,  fuel;
	{
		Solution sln;
		locations = { 2, 3, 6, 8, 4 }, start = 1, finish = 3, fuel = 5;
		auto res = sln.countRoutes(locations, start, finish, fuel);
		Assert(4, res);
	}

	{
		Solution sln;
		locations = { 4, 3, 1 }, start = 1, finish = 0, fuel = 6;
		auto res = sln.countRoutes(locations, start, finish, fuel);
		Assert(5, res);
	}

	{
		Solution sln;
		locations = { 5, 2, 1 }, start = 0, finish = 2, fuel = 3;
		auto res = sln.countRoutes(locations, start, finish, fuel);
		Assert(0, res);
	}
}

小幅优化

排序后，向右（左）第一个油量不够的城市就停止。

class Solution {
public:
	int countRoutes(vector<int>& locations, int start, int finish, int fuel) {
		const int n = locations.size();
		const int startLoc = locations[start];
		const int finishLoc = locations[finish];
		sort(locations.begin(), locations.end());
		start = std::find(locations.begin(), locations.end(), startLoc)- locations.begin();
		finish = std::find(locations.begin(), locations.end(), finishLoc) - locations.begin();
		vector<vector<C1097Int<> >> dp(n, vector<C1097Int<>>(fuel + 1));
		dp[start][0] = 1;
		for (int f = 0; f < fuel; f++)
		{
			for (int l = 0; l < n; l++)
			{
				int use = 0;
				for (int l1 = l+1 ;(l1 < n )&&( f + (use = locations[l1] - locations[l]) <= fuel); l1++)
				{
					dp[l1][f+use] += dp[l][f];					
				}
				for (int l1 = l - 1; (l1 >= 0 ) && (f + (use = locations[l] - locations[l1]) <= fuel); l1--)
				{
					dp[l1][f + use] += dp[l][f];
				}
			}
		}
		return std::accumulate(dp[finish].begin(), dp[finish].end(), C1097Int()).ToInt();
	}
};

再次优化

填表顺序，还是从使用汽油少的到使用汽油多的。
前一个位置可能在当前位置的左边，也可能是右边。不失一般性，只讨论从左边过来。
假定当前使用汽油是f，位置是l。前一站使用的汽油是f1,位置是l1。则
行驶的路径等于消耗的汽油 → \rightarrow → f-f1= l - l1 → \rightarrow → f-l = f1-l1 性质一
符合性质一的状态分以下三类：
一，使用汽油比当前使用的汽油少。此状态就是本状态的前置状态。
二，使用汽油和当前状态使用的汽油一样，就是本状态。
三，使用汽油比当前汽油多，还没有更新。
结论： 符合性质一的状态路径和，就是本状态的路径和。
从右边来，类似：
f-f1=l1-l → \rightarrow → f+l==f1+l1 性质二
m1，记录性质一；m2 记录性质二。
时间复杂度：将为O(nfuel)。

class Solution {
public:
	int countRoutes(vector<int>& locations, int start, int finish, int fuel) {
		const int n = locations.size();
		vector<vector<C1097Int<> >> dp(n, vector<C1097Int<>>(fuel + 1));
		unordered_map<int, C1097Int<>> m1,m2;
		dp[start][0] = 1;
		m1[0 - locations[start]] +=1 ;
		m2[0 + locations[start]] += 1;
		for (int f = 1; f <= fuel; f++)
		{
			for (int l = 0; l < n; l++)
			{
				dp[l][f] = m1[f- locations[l]]+ m2[f + locations[l]];
				m1[f - locations[l]] += dp[l][f];
				m2[f + locations[l]] += dp[l][f];
			}
		}
		return std::accumulate(dp[finish].begin(), dp[finish].end(), C1097Int()).ToInt();
	}
};

2023年2月第一版

class C1097Int
{
public:
C1097Int(int iData = 0) :m_iData(iData)
{
}
C1097Int operator+(const C1097Int& o)const
{
return C1097Int((m_iData + o.m_iData) % s_iMod);
}
C1097Int& operator+=(const C1097Int& o)
{
m_iData = (m_iData + o.m_iData) % s_iMod;
return this;
}
C1097Int operator(const C1097Int& o)const
{
return((long long)m_iData o.m_iData) % s_iMod;
}
C1097Int& operator=(const C1097Int& o)
{
m_iData =((long long)m_iData *o.m_iData) % s_iMod;
return *this;
}
int ToInt()const
{
return m_iData;
}
private:
int m_iData = 0;;
static const int s_iMod = 1000000007;
};

int operator+(int iData, const C1097Int& int1097)
{
int iRet = int1097.operator+(C1097Int(iData)).ToInt();
return iRet;
}

int& operator+=(int& iData, const C1097Int& int1097)
{
iData = int1097.operator+(C1097Int(iData)).ToInt();
return iData;
}

class Solution {
public:
int countRoutes(vector& locations, int start, int finish, int fuel) {
m_c = locations.size();
vector<vector> vFuelPos(fuel + 1, vector(m_c));
vFuelPos[fuel][start] = 1;
for (int iCurFuel = fuel-1; iCurFuel >= 0; iCurFuel–)
{
for (int iPos = 0; iPos < m_c; iPos++)
{
for (int iPrePos = 0; iPrePos < m_c; iPrePos++)
{
if (iPrePos == iPos)
{
continue;
}
const int iNeedFuel = iCurFuel + abs(locations[iPos] - locations[iPrePos]);
if (iNeedFuel <= fuel)
{
vFuelPos[iCurFuel][iPos] += vFuelPos[iNeedFuel][iPrePos];
}
}
}
}
C1097Int iNum = 0;
for (int iCurFuel = fuel ; iCurFuel >= 0; iCurFuel–)
{
iNum += vFuelPos[iCurFuel][finish];
}
return iNum.ToInt();
}
int m_c;
};

2023年9月版

class Solution {
public:
int countRoutes(vector& locations, int start, int finish, int fuel) {
m_c = locations.size();
int startValue = locations[start], finishValue = locations[finish];
std::sort(locations.begin(), locations.end());
start = std::lower_bound(locations.begin(), locations.end(), startValue) - locations.begin();
finish = std::lower_bound(locations.begin(), locations.end(), finishValue) - locations.begin();
m_vLeft.assign(m_c, vector<C1097Int<>>(fuel + 1));
m_vRight = m_vLeft;
int iRemain;
if ((start+1 < m_c)&&(( iRemain = fuel - locations[start+1]+ locations[start ]) >=0 ))
{
m_vRight[start+1][iRemain] = 1;
}
if ((start > 0) && ((iRemain = fuel - locations[start] + locations[start - 1]) >= 0))
{
m_vLeft[start-1][iRemain] = 1;
}
for (int iFuel = fuel-1; iFuel > 0; iFuel–)
{
for (int city = 0; city < m_c; city++)
{
RightMoveRight(iFuel, city, locations);
RightMoveLeft(iFuel, city, locations);
LeftMoveRight(iFuel, city, locations);
LeftMoveLeft(iFuel, city, locations);
}
}
C1097Int<> biRet = 0;
for (int iFuel = fuel; iFuel >= 0; iFuel–)
{
biRet += m_vLeft[finish][iFuel];
biRet += m_vRight[finish][iFuel];
}
if (start == finish)
{
biRet += 1;
}
return biRet.ToInt();
}
void RightMoveRight(int iFuel, int city, const vector& loc)
{
if (city + 1 >= m_c)
{
return;
}
const int iNeedFuel = loc[city + 1] - loc[city];
if (iNeedFuel > iFuel)
{
return;
}
m_vRight[city + 1][iFuel - iNeedFuel] += m_vRight[city][iFuel] * 2;
}
void RightMoveLeft(int iFuel, int city, const vector& loc)
{
if (0 == city)
{
return;
}
const int iNeedFuel = loc[city] - loc[city - 1];
if (iNeedFuel > iFuel)
{
return;
}
m_vLeft[city - 1][iFuel - iNeedFuel] += m_vRight[city][iFuel];
}
void LeftMoveRight(int iFuel, int city, const vector& loc)
{
if (city + 1 >= m_c)
{
return;
}
const int iNeedFuel = loc[city + 1] - loc[city];
if (iNeedFuel > iFuel)
{
return;
}
m_vRight[city + 1][iFuel - iNeedFuel] += m_vLeft[city][iFuel];
}
void LeftMoveLeft(int iFuel, int city, const vector& loc)
{
if (0 == city)
{
return;
}
const int iNeedFuel = loc[city] - loc[city - 1];
if (iNeedFuel > iFuel)
{
return;
}
m_vLeft[city - 1][iFuel - iNeedFuel] += m_vLeft[city][iFuel] * 2;
}
vector<vector<C1097Int<>>> m_vLeft, m_vRight;
int m_c;
};

扩展阅读

视频课程

有效学习：明确的目
标 及时的反馈 拉伸区（难度合适），可以先学简单的课程，请移步CSDN学院，听白银讲师（也就是鄙人）的讲解。
https://edu.csdn.net/course/detail/38771

如何你想快速形成战斗了，为老板分忧，请学习C#入职培训、C++入职培训等课程
https://edu.csdn.net/lecturer/6176

相关下载

想高屋建瓴的学习算法，请下载《喜缺全书算法册》doc版
https://download.csdn.net/download/he_zhidan/88348653

我想对大家说的话
闻缺陷则喜是一个美好的愿望，早发现问题，早修改问题，给老板节约钱。
子墨子言之：事无终始，无务多业。也就是我们常说的专业的人做专业的事。
如果程序是一条龙，那算法就是他的是睛

测试环境

操作系统：win7 开发环境： VS2019 C++17
或者 操作系统：win10 开发环境： VS2022 C++17
如无特殊说明，本算法用**C++**实现。

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：/a/357035.html

如若内容造成侵权/违法违规/事实不符，请联系我们进行投诉反馈qq邮箱809451989@qq.com，一经查实，立即删除！

相关文章

如何使用YOLOv8训练自己的模型

如何使用YOLOv8训练自己的模型

本文介绍如何用YOLO8训练自己的模型&#xff0c;我们开门见山&#xff0c;直接步入正题。 前言&#xff1a;用yolo8在自己的数据集上训练模型首先需要配置好YOLO8的环境&#xff0c;如果不会配置YOLO8环境可以参考本人主页的另一篇文章 提醒&#xff1a;使用GPU训练会大幅度加…
阅读更多...
实习日志7

实习日志7

1.试试pdf发票识别 1.1.添加文件类型判断 //判断文件类型 if (getFileType(imgCodeCell.getValue()) "jpg"||getFileType(imgCodeCell.getValue()) "png"||getFileType(imgCodeCell.getValue()) "jpeg"||getFileType(imgCodeCell.getValue(…
阅读更多...
Python爬虫解析库安装

Python爬虫解析库安装

解析库的安装 抓取网页代码之后&#xff0c;下一步就是从网页中提取信息。提取信息的方式有多种多样&#xff0c;可以使用正则来提取&#xff0c;但是写起来相对比较烦琐。这里还有许多强大的解析库&#xff0c;如 lxml、Beautiful Soup、pyquery 等。此外&#xff0c;还提供了…
阅读更多...
R语言(数据导入，清洗，可视化，特征工程，建模)

R语言(数据导入，清洗，可视化，特征工程，建模)

记录一下痛失的超级轻松的数据分析实习&#xff08;线上&#xff09;&#xff0c;hr问我有没有相关经历&#xff0c;我说我会用jupyter book进行数据导入&#xff0c;清洗&#xff0c;可视化&#xff0c;特征工程&#xff0c;建模&#xff0c;python学和用的比较多&#xff0c;…
阅读更多...
Vue学习之使用开发工具创建项目、gitcode管理项目

Vue学习之使用开发工具创建项目、gitcode管理项目

Vue学习之使用开发工具创建项目、gitcode管理项目 翻阅与学习了vue的开发工具&#xff0c;通过对比最终采用HBuilderX作为开发工具&#xff0c;以下章节对HBuilder安装与基础使用介绍 1. HBuilder 下载 从HbuildX官网&#xff08;http://www.dcloud.io/hbuilderx.html&#…
阅读更多...
HarmonyOS模拟器启动失败，电脑蓝屏解决办法

HarmonyOS模拟器启动失败，电脑蓝屏解决办法

1、在Tool->Device Manager管理界面中&#xff0c;通过Wipe User Data清理模拟器用户数据&#xff0c;然后重启模拟器&#xff1b;如果该方法无效&#xff0c;需要Delete删除已创建的Local Emulater。 2、在Tool->SDK Manager管理界面的PlatForm选项卡中&#xff0c;取消…
阅读更多...
Redis面试（三）

Redis面试（三）

1.Redis报内存不足怎么处理 Redis内存不足的集中处理方式&#xff1a; 修改配置文件redis.cof的maxmemory参数&#xff0c;增加Redis的可用内存通过命令修改set maxmemory动态设置内存上限修改内存淘汰策略&#xff0c;及时释放内存使用Redis集群&#xff0c;及时进行扩容 2…
阅读更多...
【MySQL】双写、重做日志对宕机时脏页数据落盘的作用的疑问及浅析

【MySQL】双写、重做日志对宕机时脏页数据落盘的作用的疑问及浅析

众所周知&#xff0c;双写机制、重做日志文件是mysql的InnoDB引擎的几个重要特性之二。其中两者的作用都是什么&#xff0c;很多文章都有分析&#xff0c;如&#xff0c;双写机制&#xff08;Double Write&#xff09;是mysql在crash后恢复的机制&#xff0c;而重做日志文件&am…
阅读更多...
Java 集合 05 综合练习-返回多个数据

Java 集合 05 综合练习-返回多个数据

代码&#xff1a; import java.util.ArrayList; import java.util.Arrays;public class practice{public static void main(String[] args) {ArrayList<Phone> list new ArrayList<>();Phone p1 new Phone("小米",1000);Phone p2 new Phone("苹…
阅读更多...
51单片机通过级联74HC595实现倒计时秒表Protues仿真设计

51单片机通过级联74HC595实现倒计时秒表Protues仿真设计

一、设计背景 近年来随着科技的飞速发展&#xff0c;单片机的应用正在不断的走向深入。本文阐述了51单片机通过级联74HC595实现倒计时秒表设计&#xff0c;倒计时精度达0.05s&#xff0c;解决了传统的由于倒计时精度不够造成的误差和不公平性&#xff0c;是各种体育竞赛的必备设…
阅读更多...
数据结构.栈

数据结构.栈

一、栈的定义 二、初始化 #include<iostream> using namespace std; const int N 10; typedef struct {int data[N];int top; }SqStack; void InitSqStack(SqStack &S)//初始化 {S.top -1; } 三、进栈 void Push(SqStack& S, int x)//入栈 {S.data[S.top] x; …
阅读更多...
深入了解Matplotlib中的子图创建方法

深入了解Matplotlib中的子图创建方法

深入了解Matplotlib中的子图创建方法 一 add_axes( **kwargs):1.1 函数介绍1.2 示例一 创建第一张子图1.2 示例二 polar参数的运用1.3 示例三 创建多张子图 二 add_subplot(*args, **kwargs):2.1 函数介绍2.2 示例一 三 两种方法的区别3.1 参数形式3.2 布局灵活性3.3 适用场景3…
阅读更多...
机器学习：多项式回归（Python）

机器学习：多项式回归（Python）

多元线性回归闭式解&#xff1a; closed_form_sol.py import numpy as np import matplotlib.pyplot as pltclass LRClosedFormSol:def __init__(self, fit_interceptTrue, normalizeTrue):""":param fit_intercept: 是否训练bias:param normalize: 是否标准化…
阅读更多...
verdaccio搭建npm私服

verdaccio搭建npm私服

一、安装verdaccio 注&#xff1a;加上–unsafe-perm的原因是防止报grywarn权限的错 npm install -g verdaccio --unsafe-perm 二、启动verdaccio verdaccio 三、配置文件 找到config.yml一般情况下都在用户下的这个文件夹下面 注&#xff1a;首次启动后才会生成 C:\Users\h…
阅读更多...
/etc/profile错误，命令失效

/etc/profile错误，命令失效

source /etc/profile后所有命令失效 执行 export PATH/usr/local/sbin:/usr/local/bin:/sbin:/bin:/usr/sbin:/usr/bin:/root/bin 修改后 执行:wq! 执行:w !sudo tee %
阅读更多...
怎么控制Element的数据树形表格展开所有行；递归操作，打造万能数据表格折叠。

怎么控制Element的数据树形表格展开所有行；递归操作，打造万能数据表格折叠。

HTML <el-button type"success" size"small" click"expandStatusFun"> <span v-show"expandStatusfalse"><i class"el-icon-folder-opened"></i>展开全部</span><span v-show"expan…
阅读更多...
鸿蒙原生应用开发已全面启动，你还在等什么？

鸿蒙原生应用开发已全面启动，你还在等什么？

2019年&#xff0c;鸿蒙系统首次公开亮相&#xff0c;你们说&#xff0c;等等看&#xff0c;还不成熟&#xff1b; 2021年&#xff0c;鸿蒙系统首次在手机端升级&#xff0c;你们说&#xff0c;等等看&#xff0c;还不完善&#xff1b; 2024年&#xff0c;鸿飞计划发布&#…
阅读更多...
STM32以太网接口在TCP/IP通信中的应用案例

STM32以太网接口在TCP/IP通信中的应用案例

在STM32的以太网通信中&#xff0c;TCP/IP协议广泛应用于各种领域&#xff0c;如远程监控、物联网、工业控制等。下面以一个STM32基于TCP/IP协议的以太网通信的应用案例为例进行介绍。 ✅作者简介&#xff1a;热爱科研的嵌入式开发者&#xff0c;修心和技术同步精进 ❤欢迎关注…
阅读更多...
C#颜色拾取器

C#颜色拾取器

1&#xff0c;目的&#xff1a; 获取屏幕上任意位置像素的色值。 2&#xff0c;知识点: 热键的注册与注销。 /// <summary>/// 热键注册/// </summary>/// <param name"hWnd">要定义热键的窗口的句柄 </param>/// <param name"id…
阅读更多...
如何使用Python Flask搭建一个web页面并实现远程访问

如何使用Python Flask搭建一个web页面并实现远程访问

文章目录 前言1. 安装部署Flask并制作SayHello问答界面2. 安装Cpolar内网穿透3. 配置Flask的问答界面公网访问地址4. 公网远程访问Flask的问答界面 前言 Flask是一个Python编写的Web微框架&#xff0c;让我们可以使用Python语言快速实现一个网站或Web服务&#xff0c;本期教程…
阅读更多...
最新文章

Copyright @ 2022~2023