ACM-图论-同余最短路
https://www.cnblogs.com/31415926535x/p/11692422.html
一种没见过的处理模型,,记录一下,,主要是用来处理一个多元一次方程的解的数量的问题,,数据量小时可以用看成背包处理,,数据很大时可以转换成最短路模型+一点数学来处理,,(体积模域下的最短路的问题,,求的一个最简的表示形式有模数来得到所有解
墨墨的等式
因为只是求满足的解的数量,,所以可以将方程转换成一个模方程组,,这样的方程组的解显然也是原来的解的子集,,同时可以利用模数来得到所有解,,
模数的选择是最小的那个系数,,因为如果任意选择,,会出现一些多考虑的情况
弄 mi 个点,表示从0到mi-1的所有数,,建边的方法是 i->(i+a[j])%mi
边权为 a[j] ,,表示从i这个点变成后面一个数的费用,,(因为两边都是取模的,,所以每一个数取几次后的和的余数就是那些经过的点,,也就是说一条路径就是得到一个右边为 i(mod mi) 的一个最小解,,这个最小的解就是费用和,,也就是一条最短路dis[i]
这样我们对于每一个取模的右边的B都计算一下区间里的数量,,,(计算这玩意推错了一次,,emmm
#include <bits/stdc++.h>
#define aaa cout<<233<<endl;
#define endl '\n'
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
// mt19937 rnd(TM(0));
const int inf = 0x3f3f3f3f;//1061109567 > 1e9
const ll linf = 0x3f3f3f3f3f3f3f3f;
const double eps = 1e-6;
const double pi = 3.14159265358979;
const int maxn = 3e6 + 5;
const int maxm = 1e7 + 233;
const int mod = 1e9 + 7;
ll n, l, r, a[maxn];
struct edge
{
int to, nxt; ll w;
}edge[maxn << 1];
int tot, head[maxn << 1];
void init()
{
tot = 0;
memset(head, -1, sizeof head);
}
void addedge(int u, int v, ll w)
{
edge[tot].to = v;
edge[tot].w = w;
edge[tot].nxt = head[u];
head[u] = tot++;
}
struct node
{
int v; ll w;
node(){}
node(int _v, ll _w):v(_v), w(_w){}
const bool operator<(const node &r)const{
return w > r.w;
}
}tmp;
ll dis[maxn];
bool vis[maxn];
void dijkstra(int s)
{
memset(dis, inf, sizeof dis);
memset(vis, false, sizeof vis);
priority_queue<node> q;
while(!q.empty())q.pop();
q.push(node(s, 0));
dis[s] = 0;
while(!q.empty())
{
tmp = q.top(); q.pop();
if(vis[tmp.v])continue;
vis[tmp.v] = true;
for(int i = head[tmp.v]; ~i; i = edge[i].nxt)
{
int v = edge[i].to;
if(dis[v] > dis[tmp.v] + edge[i].w)
{
dis[v] = dis[tmp.v] + edge[i].w;
q.push(node(v, dis[v]));
}
}
}
}
int main()
{
// double pp = clock();
// freopen("233.in", "r", stdin);
// freopen("233.out", "w", stdout);
ios_base::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
cin >> n >> l >> r;
for(int i = 1; i <= n; ++i)cin >> a[i];
sort(a + 1, a + 1 + n);
int mi = a[1];
init();
for(int i = 0; i <= mi - 1; ++i)
for(int j = 1; j <= n; ++j)
addedge(i, (i + a[j]) % mi, a[j]);
dijkstra(0);
ll ans = 0;
for(int i = 0; i <= mi - 1; ++i)
{
if(dis[i] <= r)
{
if(dis[i] == 0)dis[i] = mi;
ans += (r - dis[i]) / mi + 1;
if(l > dis[i])ans -= (l - dis[i] - 1) / mi + 1;
}
}
cout << ans << endl;
return 0;
}
P3403 跳楼机
比上面那个简单些,,就是注意细节,,从1开始,,有一个是1那么值一定是h,,,
#include <bits/stdc++.h>
#define aaa cout<<233<<endl;
#define endl '\n'
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
// mt19937 rnd(TM(0));
const int inf = 0x3f3f3f3f;//1061109567 > 1e9
const ll linf = 0x3f3f3f3f3f3f3f3f;
const double eps = 1e-6;
const double pi = 3.14159265358979;
const int maxn = 3e6 + 5;
const int maxm = 1e7 + 233;
const int mod = 1e9 + 7;
ll n, l, r, a[maxn];
struct edge
{
int to, nxt; ll w;
}edge[maxn << 1];
int tot, head[maxn << 1];
void init()
{
tot = 0;
memset(head, -1, sizeof head);
}
void addedge(int u, int v, ll w)
{
edge[tot].to = v;
edge[tot].w = w;
edge[tot].nxt = head[u];
head[u] = tot++;
}
struct node
{
int v; ll w;
node(){}
node(int _v, ll _w):v(_v), w(_w){}
const bool operator<(const node &r)const{
return w > r.w;
}
}tmp;
ll dis[maxn];
bool vis[maxn];
void dijkstra(int s)
{
memset(dis, inf, sizeof dis);
memset(vis, false, sizeof vis);
priority_queue<node> q;
while(!q.empty())q.pop();
q.push(node(s, 0));
dis[s] = 1;
while(!q.empty())
{
tmp = q.top(); q.pop();
if(vis[tmp.v])continue;
vis[tmp.v] = true;
for(int i = head[tmp.v]; ~i; i = edge[i].nxt)
{
int v = edge[i].to;
if(dis[v] > dis[tmp.v] + edge[i].w)
{
dis[v] = dis[tmp.v] + edge[i].w;
q.push(node(v, dis[v]));
}
}
}
}
int main()
{
// double pp = clock();
// freopen("233.in", "r", stdin);
// freopen("233.out", "w", stdout);
ios_base::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
cin >> n;
for(int i = 1; i <= 3; ++i)cin >> a[i];
sort(a + 1, a + 1 + 3);
if(a[1] == 1){
cout << n << endl;
return 0;
}
ll mi = a[1];
init();
for(int i = 0; i <= mi - 1; ++i)
for(int j = 2; j <= 3; ++j)
addedge(i, (i + a[j]) % mi, a[j]);
dijkstra(1);
ll ans = 0;
for(int i = 0; i <= mi - 1; ++i)
if(dis[i] <= n)
ans += (n - dis[i]) / mi + 1;
cout << ans << endl;
return 0;
}
遥远的旅途
这题的大致思路是将dp问题用最短路来优化,,
设 dp[i][j]
表示从起点走到i时的长度为j的一条路是否存在,,但是空间都会爆掉,,
考虑第二维,假设是通过经过若干个环来达到T,,也就是 len+kw==T
,,这里的w即为环的长度的两倍,,如果取模w就是 len%w==T%w
,,这样子dp方程就变成了到达 i 点时路径长度取模等于j的一条路径的长度,,利用spfa来转移,,只要最后 dp[n][T%w] <= T
就表示存在解,这样子利用模数来压缩了状态,,找等同的就行了,,,参考 参考
#include <bits/stdc++.h>
#define aaa cout<<233<<endl;
#define endl '\n'
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
// mt19937 rnd(TM(0));
const int inf = 0x3f3f3f3f;//1061109567 > 1e9
const ll linf = 0x3f3f3f3f3f3f3f3f;
const double eps = 1e-6;
const double pi = 3.14159265358979;
const int maxn = 1e2 + 5;
const int maxm = 1e7 + 233;
const int mod = 1e9 + 7;
ll n, m, T, a[maxn], ww;
struct edge
{
int to, nxt; ll w;
}edge[maxn << 1];
int tot, head[maxn << 1];
void init()
{
tot = 0;
memset(head, -1, sizeof head);
}
void addedge(int u, int v, ll w)
{
edge[tot].to = v;
edge[tot].w = w;
edge[tot].nxt = head[u];
head[u] = tot++;
}
ll dp[maxn][20005];
bool vis[maxn][20005];
queue<pair<int, ll> > q;
void spfa()
{
memset(dp, inf, sizeof dp);
memset(vis, false, sizeof vis);
while(!q.empty())q.pop();
dp[1][0] = 0; vis[1][0] = true;
q.push(make_pair(1, 0));
while(!q.empty())
{
int u = q.front().first; ll w = q.front().second; q.pop();
vis[u][w] = false;
for(int i = head[u]; ~i; i = edge[i].nxt)
{
int v = edge[i].to; ll vw = edge[i].w;
if(dp[v][(w + vw) % ww] > dp[u][w] + vw)
{
dp[v][(w + vw) % ww] = dp[u][w] + vw;
if(!vis[v][(w + vw) % ww])
{
vis[v][(w + vw) % ww] = true;
q.push(make_pair(v, (w + vw) % ww));
}
}
}
}
}
int main()
{
// double pp = clock();
// freopen("233.in", "r", stdin);
// freopen("233.out", "w", stdout);
ios_base::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
int t; cin >> t;
while(t--)
{
cin >> n >> m >> T;
int u, v, w;
init();
for(int i = 1; i <= m; ++i)
{
cin >> u >> v >> w;
++u, ++v;
addedge(u, v, w);
addedge(v, u, w);
}
bool flag = false;
for(int i = head[n]; ~i; i = edge[i].nxt)
{
ww = edge[i].w << 1;
spfa();
if(dp[n][T % ww] <= T)
{
flag = true;
break;
}
}
if(flag)cout << "Possible" << endl;
else cout << "Impossible" << endl;
}
return 0;
}
(end)
ACM-图论-同余最短路的更多相关文章
- 【66测试20161115】【树】【DP_LIS】【SPFA】【同余最短路】【递推】【矩阵快速幂】
还有3天,今天考试又崩了.状态还没有调整过来... 第一题:小L的二叉树 勤奋又善于思考的小L接触了信息学竞赛,开始的学习十分顺利.但是,小L对数据结构的掌握实在十分渣渣.所以,小L当时卡在了二叉树. ...
- [CF986F]Oppa Funcan Style Remastered[exgcd+同余最短路]
题意 给你 \(n\) 和 \(k\) ,问能否用 \(k\) 的所有 \(>1\) 的因子凑出 \(n\) .多组数据,但保证不同的 \(k\) 不超过 50 个. \(n\leq 10^{1 ...
- HDU 6071 Lazy Running (同余最短路 dij)
Lazy Running Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others)To ...
- BZOJ 2118 墨墨的等式 (同余最短路)
题目大意:已知B的范围,求a1x1+a2x2+...+anxn==B存在非负正整数解的B的数量,N<=12,ai<=1e5,B<=1e12 同余最短路裸题 思想大概是这样的,我们选定 ...
- 洛谷P3403跳楼机(最短路构造/同余最短路)
题目-> 解题思路: 最短路构造很神啊. 先用前两个值跑在第三个值模意义下的同余最短路(这步贪心可以证明,如果第三步长为z,那么如果n+z可以达到,n+2z同样可以达到) 最后计算与楼顶差多少个 ...
- 【CodeChef】LECOINS(同余最短路,背包DP)
题意:给定n个物品,每个物品可以取无限次,每个物品有两种属性:价值v和颜色c 现在有q个询问,每次询问是否能取出价值和为S的方案,如有多解输出不同颜色种数的最大值 题意:看到BZOJ评论区有好心人说C ...
- 【同余最短路】洛谷 P2662 牛场围栏
关于同余最短路的部分 [同余最短路]P3403跳楼机/P2371墨墨的等式 [P2662牛场围栏] 题目背景 小L通过泥萌的帮助,成功解决了二叉树的修改问题,并因此写了一篇论文, 成功报送了叉院(羡慕 ...
- 【同余最短路】【例题集合】洛谷P3403 跳楼机/P2371 墨墨的等式
接触到的新内容,[同余最短路]. 代码很好写,但思路不好理解. 同余最短路,并不是用同余来跑最短路,而是通过同余来构造某些状态,从而达到优化时间空间复杂度的目的.往往这些状态就是最短路中的点,可以类比 ...
- Codeforces 516E - Drazil and His Happy Friends(同余最短路)
Codeforces 题面传送门 & 洛谷题面传送门 首先思考一个非常简单的性质:记 \(d=\gcd(n,m)\),那么每次在一起吃完饭的男女孩编号必定与 \(d\) 同余,而根据斐蜀定理可 ...
- Codeforces 986F - Oppa Funcan Style Remastered(同余最短路)
Codeforces 题面传送门 & 洛谷题面传送门 感谢此题教会我一个东西叫做同余最短路(大雾 首先这个不同 \(k\) 的个数 \(\le 50\) 这个条件显然是让我们对每个 \(k\) ...
随机推荐
- 【LeetCode】DFS 总结
DFS(深度优先搜索) 常用来解决可达性的问题. 两个要点: 栈:用栈来保存当前节点信息,当遍历新节点返回时能够继续遍历当前节点.可以使用递归栈. 标记:和 BFS 一样同样需要对已经遍历过的节点进行 ...
- 关于volatile关键字的最佳解释
直接放原博主链接,真的写得非常好,看懂这个面试官问再多也不怕了: https://www.cnblogs.com/dolphin0520/p/3920373.html Java并发编程:volatil ...
- 剖析nsq消息队列(二) 去中心化代码源码解析
在上一篇帖子剖析nsq消息队列(一) 简介及去中心化实现原理中,我介绍了nsq的两种使用方式,一种是直接连接,还有一种是通过nslookup来实现去中心化的方式使用,并大概说了一下实现原理,没有什么难 ...
- 关于spring boot多张表建立外健的讨论
现在有四张表:student(学生表).blogs(博客表).comment(评论表).reply(回复表) 现在说一下这四张表: student(学生表):学生的信息记录表 blogs(博客表):学 ...
- 设置composer镜像地址为阿里云的方法
所有项目都会使用该镜像地址: composer config -g repo.packagist composer https://mirrors.aliyun.com/composer/ 取消配置: ...
- DBCP
DBCP(DataBase Connection Pool)数据库连接池,由Apache公司开发.连接池的运用避免了反复建立连接造成的资源浪费,预先建立一些连接放在数据库连接池中,需要时取出,不需要时 ...
- C++ std::thread概念介绍
C++ 11新标准中,正式的为该语言引入了多线程概念.新标准提供了一个线程库thread,通过创建一个thread对象来管理C++程序中的多线程. 本文简单聊一下C++多线程相关的一些概念及threa ...
- uni-app开发微信小程序的几天时间
人只有在不断的学习,才能不断的给自己充电,如果我们停止了学习,就像人没有了血脉,就会死亡,近来学习比较忙,压力比较大,整天面对着电脑,敲击代码,从中虽然收获了快乐,但是换来的确实身体的伤痛,最近虽然自 ...
- 【linux】linux固定ip
vi /etc/sysconfig/network-scripts/ifcfg-ens33 ifcfg-ens33为ifconfig显示的网卡名 TYPE="Ethernet"P ...
- 【面试必备】小伙伴栽在了JVM的内存分配策略。。。
周末有小伙伴留言说上周面试时被问到内存分配策略的问题,但回答的不够理想,小伙伴说之前公号里看过这一块的文章的,当时看时很清楚,也知道各个策略是干嘛的,但面试时脑子里清楚,心里很明白,但嘴里就是说不清楚 ...