A1958
Magic Girl Haze
T组
n个点,m条有向含权边,可以选择不超过k条边,将其权值变为0.
问点1到点n的最短距离是多少?
1≤T≤5n≤105m≤2×105k≤10wi≤109
1\leq T \leq 5 \\
n \leq 10^5 \\
m \leq 2\times10^5 \\
k \leq 10 \\
w_i \leq 10^9
1≤T≤5n≤105m≤2×105k≤10wi≤109
二维状态Dijkstra
原本dijkstra算法是dis[v],只有点的编号一个维度,dis[v]状态表示的是1到v最短的距离。
现在令dis[k][v]表示恰好k条边变0,点1到v最短的距离,状态变为二维
原本是uuu到vvv有边就可以尝试状态转移
现在是
- (k,u)(k,u)(k,u)到(k,v)(k,v)(k,v)可以转移,代价w;
- (k,u)(k,u)(k,u)到(k+1,u)(k+1,u)(k+1,u)可以转移,代价0;
正确性是显然的:
对于状态k=k0,∀u(k,u)k = k_0,\forall u(k,u)k=k0,∀u(k,u),假设它们的初始值已经设置好了,则固定k=k0k=k_0k=k0,它们之间的转移更新就是传统的一维dijkstra,因此是正确的
k=k0,∀u(k,u)k = k_0,\forall u (k,u)k=k0,∀u(k,u)初始值本是无穷大, 但是k=K0−1k=K_0-1k=K0−1的状态可以优化它们的初始值。
code
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int t;
ll n, m;
int k;
const int maxn = 1e+5 + 5;
const int maxm = 2e+5 + 5;
const int maxk = 10;
vector<int> ev[maxn];
vector<ll> ew[maxn];
// bool solved[maxk + 1][maxn];
bool vis[maxk + 1][maxn];
ll dis[maxk + 1][maxn];
struct vnd
{
int v;
ll dis;
bool operator<(const struct vnd &obj) const
{
return dis > obj.dis;
}
void print() const {
cout << "(v,dis) = "<<v<<","<<dis<<endl;
}
};
priority_queue<vnd> q[maxk + 1];
void init()
{
cin >> n >> m >> k;
for (int u = 1; u <= n; ++u)
{
ev[u].clear();
ew[u].clear();
}
int u, v, w;
for (int i = 1; i <= m; ++i)
{
cin >> u >> v >> w;
ev[u].push_back(v);
ew[u].push_back(w);
}
// memset(solved, false, sizeof(solved));
memset(vis, false, sizeof(vis));
for (int i = 0; i <= k; ++i)
q[i] = priority_queue<vnd>();
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
cin >> t;
while (t--)
{
init();
bool first = true;
ll ans;
vis[0][1] = true;
dis[0][1] = 0;
int u, v;
ll w;
for (int i = 0; i <= k; ++i)
{
// push visited vertices into queue
for (int u = 1; u <= n; ++u)
if (vis[i][u]) {
q[i].emplace(vnd{u,dis[i][u]});
}
while (!q[i].empty())
{
vnd t = q[i].top();
q[i].pop();
u = t.v;
for (size_t j = 0; j < ev[u].size(); ++j)
{
v = ev[u][j];
w = ew[u][j];
// do not change w
if (!vis[i][v])
{
vis[i][v] = true;
dis[i][v] = t.dis + w;
q[i].emplace(vnd{v,dis[i][v]});
}
else if (dis[i][v] > t.dis + w)
{
dis[i][v] = t.dis + w;
q[i].emplace(vnd{v,dis[i][v]});
}
// change w to 0
if (i + 1 > k)
continue;
if (!vis[i + 1][v])
{
vis[i + 1][v] = true;
dis[i + 1][v] = t.dis;
}
else if (dis[i + 1][v] > t.dis)
{
dis[i + 1][v] = t.dis;
}
}
}
if (vis[i][n])
{
if (first)
{
first = false;
ans = dis[i][n];
}
else
{
ans = min(ans, dis[i][n]);
}
}
}
// the input ensure: first == false
cout << ans << endl;
}
return 0;
}
A1958的更多相关文章
- ICPC 2018 南京网络赛 J Magical Girl Haze(多层图最短路)
传送门:https://nanti.jisuanke.com/t/A1958 题意:n个点m条边的路,你有k次机会将某条路上的边权变为0,问你最短路径长度 题解:最短路变形,我们需要在常规的最短路上多 ...
随机推荐
- 函数式编程/lambda表达式入门
函数式编程/lambda表达式入门 本篇主要讲解 lambda表达式的入门,涉及为什么使用函数式编程,以及jdk8提供的函数式接口 和 接口的默认方法 等等 1.什么是命令式编程 命令式编程就是我们去 ...
- java架构之路-(dubbo专题)dubbo的基本使用
今天我们来说一下dubbo的基本使用,没有什么高阶的知识,真的超级简单易懂. Dubbo核心功能解释 dubbo 阿里开源的一个SOA服务治理框架,从目前来看把它称作是一个RPC远程调用框架更为贴切. ...
- 黑龙江大学图书馆抢座python
2019.11.01 增加成功发送QQ邮件. 2019.10.31更新 晚上测试成功 ========================== 昨晚肝了一晚上,今天已经写出成品并且部署到服务器上了,代码很 ...
- window10 cmd 常见命令
AT 计划在计算机上运行的命令和程序. ATTRIB 显示或更改文件属性. BREAK 设置或清除扩展式 CTRL+C 检查. CACLS 显示或修改文件的访问控制列表(ACLs). CALL 从另一 ...
- modbus 协议说明及常用格式
--- 说明: modbus协议一般适用于一个主设备访问多个从设备的硬件开发环境,类似于zigbee网络中的一个路由器多个协调器的一对多模型. modbus常用的寄存器类型为 3X 保持寄存器和 4X ...
- postman设置变量
参数化 变量引用格式:{{username}} , 区别jmeter的 {username} 一.设置与引用环境变量 背景:在不同的环境下跑相同的测试,生产环境或测试环境 二.设置与引用全局变量 ...
- 除了降低成本和加速数字化转型,低代码还能给企业带来什么价值 ZT
翻译自:https://dzone.com/articles/measuring-the-roi-of-low-code-1,有删改 低代码 是一种近些年兴起的企业软件快速开发技术和工具.借助低代码使 ...
- JAVA JDK keytool 生成数字证书
简介: 数字证书作为网络安全数据传输的凭证,web在传输时客户端(浏览器)和 服务端(服务器)先进行会话握手,在握手过程中服务端会验证客户端的是否已经在服务端做了认证,这是单向认证.如果是双向认证的话 ...
- CSS3中新增的对文本和字体的设置
文字阴影 text-shadow: 水平偏移 垂直偏移 模糊 颜色 兼容性:IE10+ <!DOCTYPE html> <html lang="en" mani ...
- Linux配置安装
1. 安装jdk 1.1 卸载:使用java version查看虚拟机是否有jdk环境,存在先卸载: 1. 首先我的系统是CenOS7,安装完成后,先打开终端,切换到管理员账号,命令如 ...