UOJ #32. 【UR #2】跳蚤公路【Floydbellman-ford】
首先看这个范围很夸张但是其实有限制的也就在1e18*n范围里(走完一圈的边权),然后限制一定是有负环
用Floyd传递闭包,然后设f[i][j][k]为从1走了i步到j并且有k个x的最短路,用B-F处理,然后有负环就是kx+f[n][i][k]<jx+f[n-1][i][j]
对每个点求出x的限制
如果1到v的路径上有负环就不合法,所以用传递闭包出来的连通性把对当前v有限制的区间放到一起,求补集即可
判-1就是如果最后剩下的个数比1e18/100大就说明有一边是没限制的
#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=105;
const long long inf=1e18;
int n,m;
long long f[N][N][N<<1];
bool a[N][N];
vector<pair<long long,long long> >b[N],s;
struct qwe
{
int u,v,w,s;
}e[N*N];
int read()
{
int r=0,f=1;
char p=getchar();
while(p>'9'||p<'0')
{
if(p=='-')
f=-1;
p=getchar();
}
while(p>='0'&&p<='9')
{
r=r*10+p-48;
p=getchar();
}
return r*f;
}
void minn(long long &x,long long y)
{
if(x>y)
x=y;
}
void maxx(long long &x,long long y)
{
if(x<y)
x=y;
}
int main()
{
n=read(),m=read();
for(int i=1;i<=n;i++)
a[i][i]=1;
for(int i=1;i<=m;i++)
{
e[i].u=read(),e[i].v=read(),e[i].w=read(),e[i].s=read();
a[e[i].u][e[i].v]=1;
}
for(int k=1;k<=n;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(a[i][k]&&a[k][j])
a[i][j]=1;
memset(f,0x3f,sizeof(f));
// inf=f[0][0][0];
f[0][1][n]=0;
for(int i=1;i<=n;i++)
for(int k=0;k<=2*n;k++)
{
for(int j=1;j<=n;j++)
f[i][j][k]=f[i-1][j][k];
for(int j=1;j<=m;j++)
if(k-e[j].s>=0&&k-e[j].s<=2*n)
minn(f[i][e[j].v][k],f[i-1][e[j].u][k-e[j].s]+e[j].w);
}
for(int i=1;i<=n;i++)
for(int j=0;j<=2*n;j++)
if(f[n][i][j]<inf)
{
long long l=-inf,r=inf;
for(int k=0;k<=2*n;k++)
if(f[n-1][i][k]<inf)
{
if(j==k)
{
if(f[n-1][i][j]==f[n][i][j])
r=-inf,l=inf;
}
else if(j<k)
maxx(l,floor((double)(f[n-1][i][k]-f[n][i][j])/(double)(j-k))+1);
else
minn(r,ceil((double)(f[n-1][i][k]-f[n][i][j])/(double)(j-k))-1);
}
if(l<=r)
b[i].push_back(make_pair(l,r));
}
for(int i=1;i<=n;i++)
{
s.clear();
for(int j=1;j<=n;j++)
if(a[j][i])
for(int k=0,len=b[j].size();k<len;k++)
s.push_back(b[j][k]);
sort(s.begin(),s.end());
long long l=-inf,r=-inf-1,sm=2*inf+1;
for(int j=0,len=s.size();j<len;j++)
{
if(s[j].first<=r)
maxx(r,s[j].second);
else
sm-=r-l+1,l=s[j].first,r=s[j].second;
}
sm-=r-l+1;
if(sm>=inf/100)
puts("-1");
else
printf("%lld\n",sm);
}
return 0;
}
UOJ #32. 【UR #2】跳蚤公路【Floydbellman-ford】的更多相关文章
- 【UOJ#32】【UR #2】跳蚤公路(最短路)
[UOJ#32][UR #2]跳蚤公路(最短路) 题面 UOJ 题解 不难发现要求的就是是否存在负环.也就是我们只需要找到所有的负的简单环,很容易就可以想到维护路径上和\(x\)相关的内容,即维护一下 ...
- 【UR #2】跳蚤公路
[UR #2]跳蚤公路 参照yjc方法.也就是地铁环线那个题. 求每个点不在负环内的x的取值范围.然后所有1到j能到i的j的范围取交.得到答案. 每个边形如kx+b的直线,每个环也是 每个点不在负环内 ...
- Solution -「UR #2」「UOJ #32」跳蚤公路
\(\mathcal{Description}\) Link. 给定一个 \(n\) 个点 \(m\) 条边的带权有向图,每条边还有属性 \(s\in\{-1,0,1\}\).对于每个 \(u ...
- UOJ 【UR #5】怎样跑得更快
[UOJ#62]怎样跑得更快 题面 这个题让人有高斯消元的冲动,但肯定是不行的. 这个题算是莫比乌斯反演的一个非常巧妙的应用(不看题解不会做). 套路1: 因为\(b(i)\)能表达成一系列\(x(i ...
- UOJ #22 UR #1 外星人
LINK:#22. UR #1 外星人 给出n个正整数数 一个初值x x要逐个对这些数字取模 问怎样排列使得最终结果最大 使结果最大的方案数又多少种? n<=1000,x<=5000. 考 ...
- UOJ.52.[UR #4]元旦激光炮(交互 思路)
题目链接 \(Description\) 交互库中有三个排好序的,长度分别为\(n_a,n_b,n_c\)的数组\(a,b,c\).你需要求出所有元素中第\(k\)小的数.你可以调用至多\(100\) ...
- UOJ【UR #12】实验室外的攻防战
题意: 给出一个排列$A$,问是否能够经过以下若干次变换变为排列$B$ 变换:若${A_i> A_i+1}$,可以${swap(A_i,A_i+1)}$ 考虑一个数字从A排列到B排列连出来的路径 ...
- UOJ Round总结
#22. [UR #1]外星人 一开始随便搞出第一问答案,很显然的性质对$x$有变化的$a$一定是递减的,就拿一个桶直接记录可以达到的值 然后我开始想第二问,一开始想直接在这个桶上统计答案,然后发现不 ...
- 【胡策篇】题解 (UOJ 192 + CF938G + SPOJ DIVCNT2)
和泉纱雾与烟花大会 题目来源: UOJ 192 最强跳蚤 (只改了数据范围) 官方题解: 在这里哦~(说的很详细了 我都没啥好说的了) 题目大意: 求树上各边权乘积是完全平方数的路径数量. 这种从\( ...
随机推荐
- DLL中导出ANSI和UNICODE函数
模仿window中的DLL导出ANSI和UNICODE版本的函数,使用UNICODE宏来控制使用哪个版本: 在函数实际的执行代码UNICODE版本中,在ANSI函数的版本中只做参数的转换,及ANSI字 ...
- LeetCode:奇偶链表【328】
LeetCode:奇偶链表[328] 题目描述 给定一个单链表,把所有的奇数节点和偶数节点分别排在一起.请注意,这里的奇数节点和偶数节点指的是节点编号的奇偶性,而不是节点的值的奇偶性. 请尝试使用原地 ...
- MySQL 中事务的实现
在关系型数据库中,事务的重要性不言而喻,只要对数据库稍有了解的人都知道事务具有 ACID 四个基本属性,而我们不知道的可能就是数据库是如何实现这四个属性的: 在这篇文章中,我们将对事务的实现进行分析, ...
- em、pt、px和百分比
浏览器默认的字体大小为100%=16px=12pt=1em px像素(Pixel):是固定大小的单元.相对长度单位.像素px是相对于显示器屏幕分辨率而言的.一个像素等于电脑屏幕上的一个点(是你屏幕分辨 ...
- P4474 王者之剑
P4474 王者之剑 题目大意 n*m的带权网格,任意选择起点开始时刻为0秒.以下操作,每秒按顺序执行 在第i秒开始的时候,在方格(x,y)上,获得(x,y)的值 在偶数秒,周围四格的值清零 每秒可选 ...
- debian7 amd64版本添加对x86包的支持
dpkg --add-architecture i386apt-get updateapt-get install ia32-libs ia32-libs-gtk
- pip 安装模块报错解决
系统版本 ubuntu Kylin 16.04 LTS 报错1:安装pip3 安装 Django 总是提示time out,无法安装. 改用国内源: 豆瓣源: sudo pip3 install - ...
- Spring Boot2.0之 整合JDBC
很入门的知识,大家了解下就OK maven配置文件pom: spring: datasource: url: jdbc:mysql://localhost:3306/test username: ro ...
- Contiki Etimer 模块
一.Etimer概述 Etimer提供产生时间事件(timed event)的机制,当设定好的timer到期时,将会给设定etimer的process发送一个PROCESS_EVENT_TIMER 事 ...
- smokeping 出现的问题
Global symbol "%Config" requires explicit package name at /usr/lib64/perl5/lib.pm line 10. ...