[NOI2018]屠龙勇士

描述

小 D 最近在网上发现了一款小游戏。游戏的规则如下:

游戏的目标是按照编号 1∼n 顺序杀掉 n 条巨龙,每条巨龙拥有一个初始的生命值 ai 。同时每条巨龙拥有恢复能力,当其使用恢复能力时,它的生命值就会每次增加 pi,直至生命值非负。只有在攻击结束后且当生命值恰好为 0 时它才会死去。

游戏开始时玩家拥有 m 把攻击力已知的剑,每次面对巨龙时,玩家只能选择一把剑,当杀死巨龙后这把剑就会消失,但作为奖励,玩家会获得全新的一把剑。

小 D 觉得这款游戏十分无聊,但最快通关的玩家可以获得 ION2018 的参赛资格, 于是小 D 决定写一个笨笨的机器人帮她通关这款游戏,她写的机器人遵循以下规则:

每次面对巨龙时,机器人会选择当前拥有的,攻击力不高于巨龙初始生命值中攻击力最大的一把剑作为武器。如果没有这样的剑,则选择攻击力最低的一把剑作为武器。

机器人面对每条巨龙,它都会使用上一步中选择的剑攻击巨龙固定的 x 次,使巨龙的生命值减少 x×ATK。

之后,巨龙会不断使用恢复能力,每次恢复 pi 生命值。若在使用恢复能力前或某一次恢复后其生命值为 0,则巨龙死亡,玩家通过本关。

那么显然机器人的攻击次数是决定能否最快通关这款游戏的关键。小 D 现在得知了每条巨龙的所有属性,她想考考你,你知道应该将机器人的攻击次数 x 设置为多少,才能用最少的攻击次数通关游戏吗?

当然如果无论设置成多少都无法通关游戏,输出 −1 即可。

输入格式

从标准输入读入数据。

第一行一个整数 T,代表数据组数。

接下来 T 组数据,每组数据包含 5 行。

每组数据的第一行包含两个整数,n 和 m,代表巨龙的数量和初始剑的数量;

接下来一行包含 n 个正整数,第 i 个数表示第 i 条巨龙的初始生命值 ai;

接下来一行包含 n 个正整数,第 i 个数表示第 i 条巨龙的恢复能力 pi;

接下来一行包含 n 个正整数,第 i 个数表示杀死第 i 条巨龙后奖励的剑的攻击力;

接下来一行包含 m 个正整数,表示初始拥有的 m 把剑的攻击力。

输出格式

输出到标准输出中。

一共 T 行。

第 i 行一个整数,表示对于第 i 组数据,能够使得机器人通关游戏的最小攻击次数 x,如果答案不存在,输出 −1。

样例一

input

2

3 3

3 5 7

4 6 10

7 3 9

1 9 1000

3 2

3 5 6

4 8 7

1 1 1

1 1

output

59

-1

explanation

第一组数据:

开始时拥有的剑的攻击力为 {1,9,1000},第 1 条龙生命值为 3,故选择攻击力为 1 的剑,攻击 59 次,造成 59 点伤害,此时龙的生命值为 −56,恢复 14 次后生命值恰好为 0,死亡。

攻击力为 1 的剑消失,拾取一把攻击力为 7 的剑,此时拥有的剑的攻击力为 {7,9,1000},第 2 条龙生命值为 5,故选择攻击力为 7 的剑,攻击 59 次,造成 413 点伤害,此时龙的生命值为 −408,恢复 68 次后生命值恰好为 0,死亡。

此时拥有的剑的攻击力为 {3,9,1000},第 3 条龙生命值为 7,故选择攻击力为 3 的剑,攻击 59 次,造成 177 点伤害,此时龙的生命值为 −170,恢复 17 次后生命值恰好为 0,死亡。

没有比 59 次更少的通关方法,故答案为 59。

第二组数据:

不存在既能杀死第一条龙又能杀死第二条龙的方法,故无法通关,输出 −1。





excrt 板子题,就是看你细不细心,注意要乘爆longlong, 所以要做快速乘QAQ。。。


#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 5;
struct lpl{
long long atk, ai, pi, data;
}dra[maxn];
int n, m, tag;
long long mn;
long long sw[maxn];
multiset<long long> s;
multiset<long long>::iterator iter; inline long long fc(long long A, long long B, long long P){
bool tag = false; if(B < 0){B = -B; tag = true;}
long long ret = 0, tmp = A;
while(B){
if(B & 1) ret = (ret + tmp) % P;
tmp = (tmp + tmp) % P; B >>= 1;
}
return (tag) ? -ret : ret;
} long long Gcd(long long a, long long b){return (a % b == 0) ? (b) : (Gcd(b, a % b));} long long exgcd(long long a, long long b, long long &x, long long &y){
if(a == 1 && b == 0){x = 1; y = 0; return a;}
long long g = exgcd(b, a % b, x, y);
long long tmpx = x, tmpy = y;
x = tmpy; y = tmpx - a / b * tmpy; return g;
} inline void putit(){
scanf("%d%d", &n, &m); long long lin; s.clear(); tag = 0; mn = 0; s.insert(0);
for(int i = 1; i <= n; ++i) scanf("%lld", &dra[i].ai);
for(int i = 1; i <= n; ++i) scanf("%lld", &dra[i].pi);
for(int i = 1; i <= n; ++i) scanf("%lld", &sw[i]);
for(int i = 1; i <= m; ++i){scanf("%lld", &lin); s.insert(-lin);}
} inline void prepare(){
long long lin;
for(int i = 1; i <= n; ++i){
iter = s.lower_bound(-dra[i].ai); lin = *iter; if(!lin) iter--; lin = *iter;
dra[i].atk = (-1) * (lin); s.erase(iter); s.insert(-sw[i]);
mn = max(mn, ((dra[i].ai + dra[i].atk - 1) / dra[i].atk) );
}
} inline void workk(){
long long g, inv, lin;
for(int i = 1; i <= n; ++i){
g = Gcd(dra[i].atk, dra[i].pi);
if(dra[i].ai % g != 0){tag = 1; return;}
dra[i].atk /= g; dra[i].pi /= g; dra[i].ai /= g;
exgcd(dra[i].atk, dra[i].pi, inv, lin);
inv = (inv % dra[i].pi + dra[i].pi) % dra[i].pi;
dra[i].data = fc(inv, dra[i].ai, dra[i].pi);
}
} inline void excrt(){
lpl ret, lin = dra[1]; long long a, b, x, y, g;
for(int i = 2; i <= n; ++i){
lpl qwe = dra[i]; if(lin.data > qwe.data) swap(lin, qwe);
g = Gcd(lin.pi, qwe.pi); if((qwe.data - lin.data) % g != 0){tag = 1; return;}
lin.pi /= g; qwe.pi /= g;
exgcd(lin.pi, qwe.pi, x, y); ret.pi = lin.pi * qwe.pi * g;
x = fc(x, ((qwe.data - lin.data) / g), ret.pi); ret.data = fc(lin.pi * g, x, ret.pi); ret.data = (ret.data + lin.data) % ret.pi;
lin = ret;
}
if(lin.data >= mn){printf("%lld\n", lin.data); return;}
long long k = (mn - lin.data) / lin.pi; lin.data += k * lin.pi;
while(lin.data < mn) lin.data += lin.pi;
printf("%lld\n", lin.data);
} int main()
{
int T; scanf("%d", &T);
while(T--){
putit();
prepare();
workk();
if(tag){printf("-1\n"); continue;}
excrt();
if(tag){printf("-1\n"); continue;}
}
return 0;
}

uoj396 [NOI2018]屠龙勇士的更多相关文章

  1. BZOJ5418[Noi2018]屠龙勇士——exgcd+扩展CRT+set

    题目链接: [Noi2018]屠龙勇士 题目大意:有$n$条龙和初始$m$个武器,每个武器有一个攻击力$t_{i}$,每条龙有一个初始血量$a_{i}$和一个回复值$p_{i}$(即只要血量为负数就一 ...

  2. P4774 [NOI2018]屠龙勇士

    P4774 [NOI2018]屠龙勇士 先平衡树跑出打每条龙的atk t[] 然后每条龙有\(xt \equiv a[i](\text{mod }p[i])\) 就是\(xt+kp[i]=a[i]\) ...

  3. [洛谷P4774] [NOI2018]屠龙勇士

    洛谷题目链接:[NOI2018]屠龙勇士 因为markdown复制过来有点炸格式,所以看题目请戳上面. 题解: 因为杀死一条龙的条件是在攻击\(x\)次,龙恢复\(y\)次血量\((y\in N^{* ...

  4. BZOJ_5418_[Noi2018]屠龙勇士_exgcd+excrt

    BZOJ_5418_[Noi2018]屠龙勇士_exgcd+excrt Description www.lydsy.com/JudgeOnline/upload/noi2018day2.pdf 每次用 ...

  5. 洛谷 P4774 [NOI2018] 屠龙勇士

    链接:P4774 前言: 交了18遍最后发现是多组数据没清空/ll 题意: 其实就是个扩中. 分析过程: 首先发现根据题目描述的选择剑的方式,每条龙对应的剑都是固定的,有查询前驱,后继(在该数不存在前 ...

  6. 洛谷P4774 [NOI2018]屠龙勇士 [扩欧,中国剩余定理]

    传送门 思路 首先可以发现打每条龙的攻击值显然是可以提前算出来的,拿multiset模拟一下即可. 一般情况 可以搞出这么一些式子: \[ atk_i\times x=a_i(\text{mod}\ ...

  7. [NOI2018]屠龙勇士

    题目描述 题解 考虑增量法. 假设我们已经做完了前k个条件,前面的模数连乘起来的结果为M,答案为X,当前的攻击力为x,龙的血量为a. 那么我们这一次的答案的表达形式是X+t*M的. 这一次需要满足的是 ...

  8. LOJ2721 [NOI2018] 屠龙勇士 【扩展中国剩余定理】

    好久没写了,写一篇凑个数. 题目分析: 这题不难想,讲一下中国剩余定理怎么扩展. 考虑$$\left\{\begin{matrix}x \equiv a\pmod{b}\\ x \equiv c\pm ...

  9. Luogu P4774 / LOJ2721 【[NOI2018]屠龙勇士】

    真是个简单坑题...++ 前置: exgcd,exCRT,STL-multiset 读完题不难发现,攻击每条龙用的剑都是可以确定的,可以用multiset求.攻击最少显然应该对于每一条龙都操作一次,即 ...

随机推荐

  1. [Luogu1220]关路灯(区间dp)

    [Luogu1220]关路灯 题目描述 某一村庄在一条路线上安装了n盏路灯,每盏灯的功率有大有小(即同一段时间内消耗的电量有多有少).老张就住在这条路中间某一路灯旁,他有一项工作就是每天早上天亮时一盏 ...

  2. JAVA 关于File的使用

    File中常用方法 创建 createNewFile() 在指定位置创建一个空文件,成功就返回true,如果已存在就不创建然后返回false mkdir() 在指定位置创建目录,这只会创建最后一级目录 ...

  3. Java缓冲字符读取

    public class BufferedReaderDemo { public static void main(String[] args) throws IOException { // 创建流 ...

  4. Python 爬虫实战(1):分析豆瓣中最新电影的影评

    目标总览 主要做了三件事: 抓取网页数据 清理数据 用词云进行展示 使用的python版本是3.6 一.抓取网页数据 第一步要对网页进行访问,python中使用的是urllib库.代码如下: from ...

  5. 【HDOJ6628】permutation 1(dfs)

    题意:求1到n的排列中使得其差分序列的字典序为第k大的原排列 n<=20,k<=1e4 思路:爆搜差分序列,dfs时候用上界和下界剪枝 #include<bits/stdc++.h& ...

  6. LNMP环境搭建最好用的两种方法(亲测)

    经历了一个PHP服务器项目,手动编译部署PHP,Swoole环境太让人郁闷了,所以尝试过两种不错的方法,分享出来方便同样经历痛苦的coder. 第一种方式: 安装LNMP按照这里的步骤执行,网址戳我 ...

  7. 123、TensorFlow的Job

    # 如果你在分布式环境中部署TensorFlow # 你或许需要指定job name和task ID # 来将变量放置在参数服务器上 # 将操作放在worker job import tensorfl ...

  8. php中的构造函数与析构函数

    PHP面向对象——构造函数.析构函数 __construct.__destruct__construct 构造方法,当一个对象创建时调用此方法,使用此方法的好处是:可以使构造方法有一个独一无二的名称, ...

  9. Jenkins+GitLab持续集成

    向GitLab提交代码之后自动触发Jenkins构建 https://baijiahao.baidu.com/s?id=1630678692475452408&wfr=spider&f ...

  10. Spring 注解配置 WebApplicationContext

    https://docs.spring.io/spring-framework/docs/4.3.9.RELEASE/spring-framework-reference/htmlsingle/#be ...