【题解】[NOI2019Route](70分)
占坑
做法是拆掉所有式子,拆完式子看一下,如果A=0,发现边被分为了终点走向n的边和不走向n的边。所以边就有了新的边权,并且可以相加。然后通过网络流建模的套路建模使得满足时间的限制,然后由于有负边,所以跑spfa,这里60分。然后还有两个点是一条链,送10分。
先放代码,还没测,可能会爆零。交上去的爆零了因为没写freopenQAQ
//@winlere
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define ERP(t,a) for(register int t=head[a];t;t=e[t].nx)
using namespace std; typedef long long ll;
inline int qr(){
register int ret=0,f=0;
register char c=getchar();
while(c<48||c>57)f|=c==45,c=getchar();
while(c>=48&&c<=57) ret=ret*10+c-48,c=getchar();
return f?-ret:ret;
}
namespace sol1{
const int maxn=3e5+5;
struct E{
int fr,to,st,ed;
E(){to=fr=st=ed=0;}
E(const int&a,const int&b,const int&c,const int&d){
fr=a;to=b;st=c;ed=d;
}
}data[maxn];
ll w[maxn];
int arc[maxn];
int cnt;
int nodecnt;
typedef pair < ll , int > P;
//priority_queue< P , vector<P>,greater <P> > q;
queue<P> q;
bool in[maxn<<2];
vector < int > fr[maxn],to[maxn];
vector < P > e[maxn<<2];
inline bool cmp(const int&a,const int&b){
return data[a].ed>data[b].ed;
}
int n,m,A,B,C;
inline bool cmp2(const int&a,const int&b){
return data[a].st>data[b].st;
}
ll d[maxn<<2];
inline ll spfa(){
for(register int t=1;t<maxn<<2;++t) d[t]=2e18;
d[0]=0;in[0]=1;
q.push(make_pair(0,0));
while(!q.empty()){
P cur=q.front();
q.pop();in[cur.second]=0;
for(register int t=0,edd=e[cur.second].size();t<edd;++t){
//cout<<d[e[cur.second][t].second]<<' '<<d[cur.second]+e[cur.second][t].first<<endl;
if(d[e[cur.second][t].second]>d[cur.second]+e[cur.second][t].first){
d[e[cur.second][t].second]=d[cur.second]+e[cur.second][t].first;
if(!in[e[cur.second][t].second])
q.push(make_pair(d[e[cur.second][t].second],e[cur.second][t].second)),in[e[cur.second][t].second]=1;
}
}
}
return d[m+1];
}
inline void add(const int&fr,const int&to,const ll&w){
//cout<<"fr="<<fr<<' '<<"to="<<to<<" w="<<w<<endl;
e[fr].push_back(make_pair(w,to));
}
inline int init(const int&a,const int&b,const int&c,const int&d1,const int&e1){
n=a;m=b;A=c;B=d1;C=e1;
nodecnt=n;
for(register int t=1,x,y,p,q;t<=m;++t){
x=qr();y=qr();p=qr();q=qr();
//cout<<x<<' '<<y<<' '<<p<<' '<<q<<endl;
data[t]=E(x,y,p,q);
w[t]=C+1ll*B*p+1ll*A*p*p;
if(y!=n) w[t]=w[t]+1ll*A*q*q-1ll*B*q;
if(y==n) w[t]=w[t]+0ll+q;
fr[x].push_back(t);to[y].push_back(t);
if(x==1) add(0,t,w[t]);
if(y==n) add(t,m+1,0);
}
for(register int t=1;t<=n;++t){
//if(fr[t].empty()||to[t].empty())continue;
sort(to[t].begin(),to[t].end(),cmp);
sort(fr[t].begin(),fr[t].end(),cmp2);
int k=0;/*
cout<<t<<' '<<fr[t].size()<<' '<<to[t].size()<<endl;
cout<<"fr=";
for(auto f:fr[t])
cout<<f<<" ";
putchar('\n');
cout<<"to=";
for(auto f:to[t])
cout<<f<<" ";
putchar('\n');
*/
for(register int i=0,edd=to[t].size(),ed2=fr[t].size();i<edd;++i){
if(i<edd-1)add(to[t][i+1],to[t][i],0);
while(k<ed2&&data[fr[t][k]].st>=data[to[t][i]].ed) add(to[t][i],fr[t][k],w[fr[t][k]]),++k;
}
}
cout<<spfa()<<endl;
return 0;
}
}
int main(){
freopen("route.in","r",stdin);
freopen("route.out","w",stdout);
int a=qr(),b=qr(),c=qr(),d=qr(),e=qr();
if(a<=100&&b==a-1&&c!=0){
ll A=c,B=d,C=e;
ll ans=1ll*a*C;
ll t1=0,t2=0,t3=0;
for(register int t=1;t<=b;++t){
qr();qr();t1=qr();t2=qr();
ans+=B*(t1-t3)+A*(t1-t3)*(t1-t3);
t3=t2;
}
cout<<ans+t3<<endl;
return 0;
}
return sol1::init(a,b,c,d,e);
return 0;
}
【题解】[NOI2019Route](70分)的更多相关文章
- 洛谷P4559 [JSOI2018]列队 【70分二分 + 主席树】
题目链接 洛谷P4559 题解 只会做\(70\)分的\(O(nlog^2n)\) 如果本来就在区间内的人是不用动的,区间右边的人往区间最右的那些空位跑,区间左边的人往区间最左的那些空位跑 找到这些空 ...
- 「PKUSC2018」星际穿越 (70分做法)
5371: [Pkusc2018]星际穿越 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 27 Solved: 11[Submit][Status] ...
- 洛谷P1979 华容道(70分 暴力)
P1979 华容道 题目描述 [问题描述] 小 B 最近迷上了华容道,可是他总是要花很长的时间才能完成一次.于是,他想到用编程来完成华容道:给定一种局面, 华容道是否根本就无法完成,如果能完成, 最少 ...
- 洛谷P1081 开车旅行70分
https://www.luogu.org/problem/show?pid=1081 太遗憾了明明写出来了,却把最小值初始值弄小了,从第二个点开始就不可能对了.70分! #include<io ...
- 华容道 noip2013 70分搜索
题目描述 [问题描述] 小 B 最近迷上了华容道,可是他总是要花很长的时间才能完成一次.于是,他想到用编程来完成华容道:给定一种局面, 华容道是否根本就无法完成,如果能完成, 最少需要多少时间. 小 ...
- ccf 201712-4 行车路线(70分)
ccf 201712-4 行车路线 解题思路: 首先Dijkstra是基于贪心算法的,即每一次作出的选择都具有贪心选择性.此题由于有“如果连续走小道,小明的疲劳值会快速增加,连续走s公里小明会增加s2 ...
- Cogs 1264. [NOIP2012] 开车旅行(70分 暴力)
1264. [NOIP2012] 开车旅行 ★★☆ 输入文件:drive.in 输出文件:drive.out 简单对比时间限制:2 s 内存限制:128 MB [题目描述] 小A 和小 ...
- 【NOI2016】优秀的拆分 题解(95分)
题目大意: 求一个字符串中形如AABB的子串个数. 思路: 用哈希做到O(1)判断字符串是否相同,O($n^2$)预处理,ans[i]为开头位置为i的形如AA的子串个数.再用O($n^2$)枚举出AA ...
- ACM: 限时训练题解-Runtime Error-二分查找
Runtime Error Bahosain was trying to solve this simple problem, but he got a Runtime Error on one ...
随机推荐
- 创建我的flask第一个应用(二)
继上一篇创建我的flask第一个应用(一),继续学习配置flask 在myproject未提供flask默认运行的主程序文件"wsgi.py"或"app.py" ...
- Twitter 宣布抛弃 Mesos,全面转向Kubernetes
摘要: 从最早Mesos“代言人”到如今的全面转向“Kubernetes Native”,Twitter的举动再一次佐证了‘Kuberentes已经成为容器编排事实标准’这一断言. 本文作者:张磊 阿 ...
- Xib设计UITableViewCell然后动态加载
转自: http://www.2cto.com/kf/201202/120764.html (注:环境Mac OS X Lion 10.7.3 + Xcode 4.2.1 + iOS SDK 5.0. ...
- thinkphp3.2配置redis缓存和文件缓存
如果把一些常用但又不容易变的数据存缓存,而不是每次查数据库,这样能很大减轻数据库压力 最近由于项目需要,就尝试了一把redis,但是后面又用了tp3.2的文件缓存,直接进入主题: 在config.ph ...
- 注意 Laravel 清除缓存 php artisan cache:clear 的一个坑
Laravel 的命令 php artisan cache:clear 用来清除各种缓存,如页面,Redis,配置文件等缓存,它会清空 Redis 数据库的全部数据,比如默认使用的 Redis 的 数 ...
- pytorch更新
Pytorch如何更新版本与卸载,使用pip,conda更新卸载Pytorch 2018年05月22日 07:33:52 醉雨轩Y 阅读数 19047 今天我们主要汇总如何使用使用ubuntu,C ...
- Python--day68--Django ORM的字段参数、元信息
字段参数 null 用于表示某个字段可以为空. unique 如果设置为unique=True 则该字段在此表中必须是唯一的 . db_index 如果db_index=True 则代表着为此字段设置 ...
- 三种查看MySQL数据库版本的方法
https://blog.csdn.net/hj7jay/article/details/51921504 1.使用-V参数 首先我们想到的肯定就是查看版本号的参数命令,参数为-V(大写字母)或者-- ...
- uni-app 快速认识
uni-app 是一个使用 Vue.js 开发所有前端应用的框架,开发者编写一套代码,可发布到iOS.Android.H5.以及各种小程序(微信/阿里/百度/头条/QQ)等多个平台. 即使不跨端,un ...
- Oracle备库宕机启动解决方案
简介 ORA-10458: standby database requires recovery ORA-01196: 文件 1 由于介质恢复会话失败而不一致 ORA-01110: 数据文件 1: ' ...