洛谷P4774 BZOJ5418 LOJ2721 [NOI2018]屠龙勇士(扩展中国剩余定理)
题目链接:
题目大意:这么长的题面,就饶了我吧emmm
这题第一眼看上去没法列出同余方程组。为什么?好像不知道用哪把剑杀哪条龙……
仔细一看,要按顺序杀龙,所以获得的剑出现的顺序也是固定的。
那么如果能把所有龙杀死,就能模拟出哪把剑杀那条龙了。
(以下设所有除 $n,m$ 外的数的最大值为 $v$)
$O(nm)$?
不,发现这里用剑的限制实际上是给出一个上界,来用lower_bound的。
插入也不要太暴力。我们想到什么?手写平衡树multiset!
这一部分复杂度是 $O(n\log m)$。
接下来就成了一个同余方程组:$atk_ix\equiv hp_i(\operatorname{mod}\ rec_i)$。
($atk_i$ 表示攻击第 $i$ 条龙的攻击力,$hp_i$ 表示第 $i$ 条龙的血量,$rec_i$ 表示第 $i$ 条龙的恢复能力)
但是好像没那么简单!有一些神奇的情况……要特判!
在考场上的话我肯定想不到特判任何东西
首先大家应该都知道,$x\equiv a_i(\operatorname{mod}\ b_i)$ 等价于 $x\equiv a_i\ \operatorname{mod}\ b_i(\operatorname{mod}\ b_i)$。
但是在这题中不行!你总不能把一头龙的血量凭空减少吧!
所以 $hp_i>rec_i$ 时(题目描述中的特性1不满足时)怎么办?
我们发现不满足特性1的点都满足 $rec_i=1$……也就是只要把他的血量打到小于0就赢了!
此时答案是 $\max\limits_{1\le i\le n}(\lceil\frac{hp_i}{atk_i}\rceil)$。
另外要注意 $x\ge\max(\lceil\dfrac{hp_i}{atk_i}\rceil)$。
至此所有特判结束。
我们一般的EXCRT同余方程组的形式是 $x\equiv a_i(\operatorname{mod}\ b_i)$。
但是这里有讨厌的 $atk_i$!更讨厌的是可能没有逆元,不能直接除过去!
我们来看看如何变形:
$atk_ix\equiv hp_i(\operatorname{mod}\ rec_i)$
$atk_ix+rec_iy=hp_i$
此时可以EXGCD解出 $x$ 的最小正整数解 $x_0$。
根据某定理(什么定理来着?)可得通解为 $x=x_0+\frac{rec_i}{\gcd(atk_i,rec_i)}k(k\in N^+)$。
于是就变成了 $x\equiv x_0(\operatorname{mod}\ \frac{rec_i}{\gcd{atk_i},rec_i})$。
这一部分复杂度为 $O(n\log v)$。
接下来就是裸的EXCRT了!
这一部分复杂度为 $O(n\log v)$。
等一下,还有一句提示没看到:
你所用到的中间结果可能很大,注意保存中间结果的变量类型。
突然慌张
冷静分析一下,我们在各种算法中模数都有可能超过int范围,要是一乘……
那……只能上龟速乘了。时间复杂度是没变,但是常数的话……
幸好不是wys的题,不然时限就是1s了
好吧,都说完了,上代码吧。(之前因为数据水没有注意 $x\ge\max(\lceil\dfrac{hp_i}{atk_i}\rceil)$)
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=;
multiset<ll> mts; //模拟拿剑
int t,n,m;
bool spec,exist,same;
ll hp[maxn],rec[maxn],a[maxn],b[maxn]; //a,b表示最后方程组是x=a_i(mod\ b_i)
int award[maxn],atk[maxn]; //award是奖励的剑
ll qmul(ll a,ll b,ll mod){ //*速乘
ll ans=;
for(;b;b>>=,a=(a+a)%mod) if(b&) ans=(ans+a)%mod;
return ans;
}
ll exgcd(ll a,ll b,ll &x,ll &y){ //模板
if(!b){x=;y=;return a;}
ll d=exgcd(b,a%b,y,x);y-=a/b*x;return d;
}
ll gcd(ll a,ll b){ //(不要问我为什么有这个)
if(!b) return a;
return gcd(b,a%b);
}
int main(){
scanf("%d",&t);
while(t--){
spec=false;same=exist=true; //spec是rec=1,same是hp=rec,exist是有没有解
mts.clear();
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++) scanf("%lld",hp+i);
for(int i=;i<=n;i++){scanf("%lld",rec+i);if(rec[i]<hp[i]) spec=true;} //rec<hp?
for(int i=;i<=n;i++) scanf("%d",award+i);
for(int i=;i<=m;i++){
int x;
scanf("%d",&x);
mts.insert(x); //全部进去
}
for(int i=;i<=n;i++){
multiset<ll>::iterator it;
if(hp[i]<*mts.begin()) it=mts.begin(); //如果没有比血量小的剑,则用第一把(最小的)
else it=mts.upper_bound(hp[i]),it--; //否则用第一个大于它的-1(即最后一个小于等于它的)
atk[i]=*it;
mts.erase(it);
mts.insert(award[i]); //拿走再选
}
if(spec){ //特判hp>rec
ll ans=;
for(int i=;i<=n;i++) ans=max(ans,ll(ceil(1.0*hp[i]/atk[i])));
printf("%lld\n",ans);
continue;
}
for(int i=;i<=n;i++){
if(!atk[i]){puts("-1");exist=false;break;} //无法攻击?无解
if(hp[i]!=rec[i]) same=false; //hp=rec?
ll x,y,d=exgcd(atk[i],rec[i],x,y); //解x0
if(hp[i]%d){puts("-1");exist=false;break;} //无解
x=(x+rec[i])%rec[i];b[i]=rec[i]/d;a[i]=qmul(x,hp[i]/d,b[i]); //真的x0推出b=rec/gcd(atk,rec),a=x0
}
if(same){ //特判hp=rec
ll ans=ll(ceil(1.0*hp[]/atk[]));
for(int i=;i<=n;i++){
ll x=ceil(1.0*hp[i]/atk[i]),d=gcd(ans,x);
ans=ans/d*x;
}
printf("%lld\n",ans);continue;
}
if(!exist) continue;
ll ans=a[],mod=b[]; //开始EXCRT(以下模板,不解释,可以去我的另一个博客看)
for(int i=;i<=n;i++){
ll x,y,d=exgcd(mod,b[i],x,y),r=((a[i]-ans)%b[i]+b[i])%b[i],tmp=mod/d*b[i];
if(r%d){puts("-1");exist=false;break;}
x=(x+b[i])%b[i];x=qmul(x,r/d,b[i]);
ans=(qmul(mod,x,tmp)+ans)%tmp;
mod=tmp;
}
if(!exist) continue;
printf("%lld\n",ans); //终于没了
}
}
NOI2018D2T1_dragon(EXCRT)
终于啊……
洛谷P4774 BZOJ5418 LOJ2721 [NOI2018]屠龙勇士(扩展中国剩余定理)的更多相关文章
- 洛谷 P4774 / loj 2721 [NOI2018] 屠龙勇士 题解【同余】【exgcd】【CRT】
推导过程存在漏洞+exCRT板子没打熟于是期望得分÷实际得分=∞? 题目描述 小 D 最近在网上发现了一款小游戏.游戏的规则如下: 游戏的目标是按照编号 \(1\sim n\) 顺序杀掉 \(n\ ...
- LOJ.2721.[NOI2018]屠龙勇士(扩展CRT 扩展欧几里得)
题目链接 LOJ 洛谷 rank前3无压力(话说rank1特判打表有意思么) \(x*atk[i] - y*p[i] = hp[i]\) 对于每条龙可以求一个满足条件的\(x_0\),然后得到其通解\ ...
- LOJ2721 [NOI2018] 屠龙勇士 【扩展中国剩余定理】
好久没写了,写一篇凑个数. 题目分析: 这题不难想,讲一下中国剩余定理怎么扩展. 考虑$$\left\{\begin{matrix}x \equiv a\pmod{b}\\ x \equiv c\pm ...
- 扩展中国剩余定理学习笔记+模板(洛谷P4777)
题目链接: 洛谷 题目大意:求同余方程组 $x\equiv b_i(mod\ a_i)$ 的最小正整数解. $1\leq n\leq 10^5,1\leq a_i\leq 10^{12},0\leq ...
- [洛谷P4774] [NOI2018]屠龙勇士
洛谷题目链接:[NOI2018]屠龙勇士 因为markdown复制过来有点炸格式,所以看题目请戳上面. 题解: 因为杀死一条龙的条件是在攻击\(x\)次,龙恢复\(y\)次血量\((y\in N^{* ...
- (伪)再扩展中国剩余定理(洛谷P4774 [NOI2018]屠龙勇士)(中国剩余定理,扩展欧几里德,multiset)
前言 我们熟知的中国剩余定理,在使用条件上其实是很苛刻的,要求模线性方程组\(x\equiv c(\mod m)\)的模数两两互质. 于是就有了扩展中国剩余定理,其实现方法大概是通过扩展欧几里德把两个 ...
- P4774 [NOI2018]屠龙勇士
P4774 [NOI2018]屠龙勇士 先平衡树跑出打每条龙的atk t[] 然后每条龙有\(xt \equiv a[i](\text{mod }p[i])\) 就是\(xt+kp[i]=a[i]\) ...
- BZOJ5418[Noi2018]屠龙勇士——exgcd+扩展CRT+set
题目链接: [Noi2018]屠龙勇士 题目大意:有$n$条龙和初始$m$个武器,每个武器有一个攻击力$t_{i}$,每条龙有一个初始血量$a_{i}$和一个回复值$p_{i}$(即只要血量为负数就一 ...
- BZOJ_5418_[Noi2018]屠龙勇士_exgcd+excrt
BZOJ_5418_[Noi2018]屠龙勇士_exgcd+excrt Description www.lydsy.com/JudgeOnline/upload/noi2018day2.pdf 每次用 ...
随机推荐
- python图像处理模块Pillow--Image模块
一.简介 PIL:Python Imaging Library,已经是Python平台事实上的图像处理标准库了.PIL功能非常强大,但API却非常简单易用 由于PIL仅支持到Python 2.7,加上 ...
- VB6 选择文件夹路径
'--------------------------------------------------------------------------------------- ' Module : ...
- angularJs 技巧总结及最佳实践
强烈建议通读官方wiki文档,里面包含了FAQ,最佳实践,深入理解最核心的Directive及Scope等文章, 基础 1. 使用ng-repeat指令,为防止重复值发生的错误.加上track by ...
- libgdx学习记录8——对话框Dialog
Dialog在游戏中也很常用,尤其在设置.退出.商店.暂停等画面.Dialog的使用也可以通过skin实现,也可以自定义. 下面是一个简单的实例: package com.fxb.newtest; i ...
- ajax传参data里面的键是一个变量的解决方法
直接用这种方式来传参,比如bean中有字段 username password,则是 data[username] = "用户名"; data[password] = " ...
- MQ见解
MQ 消息队列是系统级.模块级的通信.RPC是对象级.函数级通信. 1) 什么是推模式,什么是拉模式 2) 有没有消息丢失情况,如何防止 3) MQ用来解决什么问题 4) 你们用的什么MQ,为什么要用 ...
- HTTP协议图--概述
1.计算机网络体系结构分层 2.TCP/IP 通信传输流 利用 TCP/IP 协议族进行网络通信时,会通过分层顺序与对方进行通信.发送端从应用层往下走,接收端则从链路层往上走.如下: 3.TCP/ ...
- From today 2019.02.27
HIT开设软件构造课程,需要在博客上分享记录学习体验,感觉还是挺好的. 以后会不定期更新一些关于学下java的笔记和实验相关的内容.
- 小学生都能写智能语音助手了,我这颗转战AI的心要何去何从?
前言——我是不是老了 前天看了一个关于AI类的综艺节目我感觉整个人都不好了.这个综艺的名字叫<智造将来>上面那个小屁孩自己写了一个智能语音助手,这个小屁孩叫袁翊闳是2018年百度AI开发者 ...
- PAT甲题题解-1041. Be Unique (20)-水题
博主欢迎转载,但请给出本文链接,我尊重你,你尊重我,谢谢~http://www.cnblogs.com/chenxiwenruo/p/6789189.html特别不喜欢那些随便转载别人的原创文章又不给 ...