【CodeForces】576 D. Flights for Regular Customers
【题目】D. Flights for Regular Customers
【题意】给定n个点m条边的有向图,每条边有di表示在经过该边前必须先经过di条边,边可重复经过,求1到n的最小经过边数。n,m<=150,di<=10^9,time=4s。
【算法】floyd+矩阵快速幂
【题解】需要计算步数,很容易联想到将floyd中每一步拆成矩阵乘法的经典做法。
令a[d][i][j]表示恰好d步能否从 i 走到 j(邻接矩阵) ,令b[d][i][j]表示当前已走d步时允许通过的边(连边矩阵)。
当前d步时,走一步相当于a[d+1][i][j]=max(a[d][i][k]&&b[d][k][j]),k=1~n(取max 是 或运算),将其定义为一次矩阵乘法就可以用矩阵快速幂优化了。
根据d的不同,连边矩阵b会发生至多m次变化。所以将所有边按d排序后顺序处理,每次快速幂的幂为d[i]-d[i-1],就能处理出走恰好d步能到达的点(作为起点集)。
最后,每次用1~i的所有边跑一次朴素的最短路(可以用floyd),从起点集到n的最短路+d就是每次的答案,取最小答案。
因为矩阵快速幂中使用的均为01矩阵,所以可以用bitset优化,从而满足时间限制。
复杂度O(m*n^3*log(di)/32+m*n^3),其中第二部分的最短路还可以用dijkstra算法优化,但已经不必要了。
#include<cstdio>
#include<cstring>
#include<bitset>
#include<algorithm>
using namespace std;
const int N=,inf=0x3f3f3f3f;
struct cyc{int x,y,d;}e[N];
int n,m,mp[N][N];
bitset<N>a[N],ans[N],c[N],z[N];
bool cmp(cyc a,cyc b){return a.d<b.d;}
void mul(bitset<N>a[N],bitset<N>b[N]){
for(int i=;i<=n;i++)c[i].reset();
for(int k=;k<=n;k++)
for(int i=;i<=n;i++)
if(a[i][k])c[i]|=b[k];
for(int i=;i<=n;i++)a[i]=c[i];
}
int main(){
scanf("%d %d",&n,&m);
for(int i=;i<=m;i++)scanf("%d %d %d",&e[i].x,&e[i].y,&e[i].d);
sort(e+,e+m+,cmp);
for(int i=;i<=n;i++)ans[i].reset();
for(int i=;i<=n;i++)ans[i][i]=;
int ANS=inf;
for(int b=;b<=m;b++){
for(int i=;i<=n;i++)a[i].reset();
for(int i=;i<b;i++)a[e[i].x][e[i].y]=;
int k=e[b].d-e[b-].d;
while(k){
if(k&)mul(ans,a);
mul(a,a);
k>>=;
}
memset(mp,0x3f,sizeof(mp));
for(int i=;i<=n;i++)mp[i][i]=;
for(int j=;j<=b;j++)mp[e[j].x][e[j].y]=;
for(int k=;k<=n;k++)
for(int i=;i<=n;i++)
for(int j=;j<=n;j++)
mp[i][j]=min(mp[i][j],mp[i][k]+mp[k][j]);
for(int i=;i<=n;i++)if(ans[][i])ANS=min(ANS,e[b].d+mp[i][n]);
}
if(ANS==inf)printf("Impossible");else printf("%d",ANS);
return ;
}
代码中要注意的是,因为这里矩阵乘法并不是通常意义上定义的,所以初始值的设置不能沿袭“单位矩阵”等概念,而必须考虑其实际意义。
初始邻接矩阵:全为0,对角线为1,表示初始只能到达自己。(第二轮后可能就不能到达自己了)
连边矩阵:对角线为0,表示不能通过连边到达自己,这样才能使邻接矩阵表示恰好d步的可达信息。
如果连边矩阵的对角线1,即能到达自己,那么表示的就是<=d步的可达信息。
【CodeForces】576 D. Flights for Regular Customers的更多相关文章
- 【CodeForces】576 B. Invariance of Tree
[题目]B. Invariance of Tree [题意]给定n个数的置换,要求使n个点连成1棵树,满足u,v有边当且仅当a[u],a[v]有边,求一种方案或无解.n<=10^5. [算法]数 ...
- 【CodeForces】576 C. Points on Plane
[题目]C. Points on Plane [题意]给定坐标系中n个点的坐标(范围[0,10^6]),求一种 [ 连边形成链后总长度<=2.5*10^9 ] 的方案.n<=10^6. [ ...
- 【Codeforces】CF 5 C Longest Regular Bracket Sequence(dp)
题目 传送门:QWQ 分析 洛谷题解里有一位大佬讲的很好. 就是先用栈预处理出可以匹配的左右括号在数组中设为1 其他为0 最后求一下最长连续1的数量. 代码 #include <bits/std ...
- 【Codeforces】Round #491 (Div. 2) 总结
[Codeforces]Round #491 (Div. 2) 总结 这次尴尬了,D题fst,E没有做出来.... 不过还好,rating只掉了30,总体来说比较不稳,下次加油 A:If at fir ...
- 【Codeforces】Round #488 (Div. 2) 总结
[Codeforces]Round #488 (Div. 2) 总结 比较僵硬的一场,还是手速不够,但是作为正式成为竞赛生的第一场比赛还是比较圆满的,起码没有FST,A掉ABCD,总排82,怒涨rat ...
- Codeforces 576D Flights for Regular Customers(矩阵加速DP)
题目链接 Flights for Regular Customers 首先按照$d$的大小升序排序 然后分成$m$个时刻,每条路径一次处理过来. $can[i][j]$表示当前时刻$i$能否走到$j ...
- CF576D Flights for Regular Customers 矩阵乘法 + Bitset优化
%%%cxhscst2's blog Codeforces 576D Flights for Regular Customers(矩阵加速DP) 代码非常优美 + 简洁,学习到了 Code: #inc ...
- 「CF576D」 Flights for Regular Customers
「CF576D」 Flights for Regular Customers 对不起我又想网络流去了 你看这长得多像啊,走过至少多少条边就是流量下界,然后没上界 但是这个题求的最少走多少条边啊...完 ...
- 【CodeForces】601 D. Acyclic Organic Compounds
[题目]D. Acyclic Organic Compounds [题意]给定一棵带点权树,每个点有一个字符,定义一个结点的字符串数为往下延伸能得到的不重复字符串数,求min(点权+字符串数),n&l ...
随机推荐
- Learn Docker(一)—软件安装与常规操作
一.安装Docker Windows平台 在Windows10 X64专业版上可以直接下载Docker原生应用进行安装,在控制面板的程序与功能里启用Hyper-v,之后就可以运行docker程序啦. ...
- OSG学习:使用OSG中预定义的几何体
常用的内嵌几何体包括: osg::Box //正方体 osg::Capsule //太空舱 osg::Cone //椎体 osg::Cylinder //柱体 osg::HeightField //高 ...
- docker 开启远程
# vi /etc/init.d/docker 在start()中加入:$exec -H tcp://0.0.0.0:2375 -H unix://var/run/docker.sock -d & ...
- 在linux下编译线程程序undefined reference to `pthread_create'
由于是Linux新手,所以现在才开始接触线程编程,照着GUN/Linux编程指南中的一个例子输入编译,结果出现如下错误:undefined reference to 'pthread_create'u ...
- poj3107 Godfather 求树的重心
Description Last years Chicago was full of gangster fights and strange murders. The chief of the pol ...
- python 内置函数02
1. lambda 匿名函数 lambda 参数: 返回值 #常规计算两个数相加的函数 def func(a,b): return a+b print(func(1,9)) #lambda函数 my_ ...
- Splitter Control for Dialog
原文链接地址:https://www.codeproject.com/Articles/595602/Splitter-Control-for-Dialog Introduction Yes, tha ...
- WlanGetAvailableNetworkList
原文msdn链接地址:https://docs.microsoft.com/zh-cn/windows/desktop/api/wlanapi/nf-wlanapi-wlangetavailablen ...
- bzoj5055: 膜法师(BIT)
大水题WA了两发T T 记录一下a[i]的前缀和,a[i]*a[j]就是sigma(a[j]*sumi[j-1]) 记录一下a[i]*a[j]的前缀和,a[i]*a[j]*a[k]就是sigma(a[ ...
- 公告:开通csdn博客,敬请关注!
公告:开通csdn博客,敬请关注!地址:https://blog.csdn.net/cyjch