P4774 [NOI2018]屠龙勇士

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

游戏的目标是按照编号 $\rightarrow n1→n 顺序杀掉 nnn 条巨龙，每条巨龙拥有一个初始的生命值 aia_iai 。同时每条巨龙拥有恢复能力，当其使用恢复能力时，它的生命值就会每次增加 pip_ipi ，直至生命值非负。只有在攻击结束后且当生命值恰好为 000 时它才会死去。$
游戏开始时玩家拥有
每次面对巨龙时，机器人会选择当前拥有的，攻击力不高于巨龙初始生命值中攻击力最大的一把剑作为武器。如果没有这样的剑，则选择 攻击力最低 的一把剑作为武器。
机器人面对每条巨龙，它都会使用上一步中选择的剑攻击巨龙固定的 $\times ATKx×ATK 。$
之后，巨龙会不断使用恢复能力，每次恢复 $pip_ipi 生命值。若在使用恢复能力前或某一次恢复后其生命值为 000 ，则巨龙死亡，玩家通过本关。$

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

当然如果无论设置成多少都无法通关游戏，输出

第一行一个整数

接下来

每组数据的第一行包含两个整数，
接下来一行包含 $aia_iai ;$
接下来一行包含 $pip_ipi ;$
接下来一行包含
接下来一行包含

一共

第

第一组数据：

开始时拥有的剑的攻击力为 ${1,9,10}\{1,9,10\}{1,9,10}，第 111 条龙生命值为 333，故选择攻击力为 111 的剑，攻击 595959 次，造成 595959 点伤害，此时龙的生命值为 −56-56−56，恢复 14 次后生命值恰好为 000，死亡。$
攻击力为 ${7,9,10}\{7,9,10\}{7,9,10}，第 2 条龙生命值为 555，故选择攻击力为 777 的剑，攻击 595959 次，造成 413413413 点伤害，此时龙的生命值为 −408-408−408，恢复 686868 次后生命值恰好为 000，死亡。$
此时拥有的剑的攻击力为 ${3,9,10}\{3,9,10\}{3,9,10}，第 333 条龙生命值为 777，故选择攻击力为 333 的剑，攻击 595959 次，造成 177177177 点伤害，此时龙的生命值为 −170-170−170，恢复 171717 次后生命值恰好为 0，死亡。$
没有比

第二组数据：不存在既能杀死第一条龙又能杀死第二条龙的方法，故无法通关，输出

【子任务】

测试点编号			$pip_ipi$	$aia_iai$	攻击力	其他限制
1	$≤105\le 10^5≤105$			$≤105\le 10^5≤105$		无
2	$≤105\le 10^5≤105$			$≤105\le 10^5≤105$		无
3	$≤105\le 10^5≤105$			$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	无
4	$≤105\le 10^5≤105$			$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	无
5	$≤103\le 10^3≤103$	$≤103\le 10^3≤103$	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	特性 1、特性 2
6	$≤103\le 10^3≤103$	$≤103\le 10^3≤103$	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	特性 1、特性 2
7	$≤103\le 10^3≤103$	$≤103\le 10^3≤103$	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	特性 1、特性 2
8			$≤108\le 10^8≤108$	$≤108\le 10^8≤108$	$≤106\le 10^6≤106$	特性 1
9			$≤108\le 10^8≤108$	$≤108\le 10^8≤108$	$≤106\le 10^6≤106$	特性 1
10			$≤108\le 10^8≤108$	$≤108\le 10^8≤108$	$≤106\le 10^6≤106$	特性 1
11			$≤108\le 10^8≤108$	$≤108\le 10^8≤108$	$≤106\le 10^6≤106$	特性 1
12			$≤108\le 10^8≤108$	$≤108\le 10^8≤108$	$≤106\le 10^6≤106$	特性 1
13			$≤108\le 10^8≤108$	$≤108\le 10^8≤108$	$≤106\le 10^6≤106$	特性 1
14	$=105=10^5=105$	$=105=10^5=105$		$≤108\le 10^8≤108$	$≤106\le 10^6≤106$	无特殊限制
15	$=105=10^5=105$	$=105=10^5=105$		$≤108\le 10^8≤108$	$≤106\le 10^6≤106$	无特殊限制
16	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	所有 $pip_ipi 是质数$	$≤1012\le 10^{12}≤1012$	$≤106\le 10^6≤106$	特性 1
17	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	所有 $pip_ipi 是质数$	$≤1012\le 10^{12}≤1012$	$≤106\le 10^6≤106$	特性 1
18	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	无特殊限制	$≤1012\le 10^{12}≤1012$	$≤106\le 10^6≤106$	特性 1
19	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	无特殊限制	$≤1012\le 10^{12}≤1012$	$≤106\le 10^6≤106$	特性 1
20	$≤105\le 10^5≤105$	$≤105\le 10^5≤105$	无特殊限制	$≤1012\le 10^{12}≤1012$	$≤106\le 10^6≤106$	特性 1

特性 1 是指：对于任意的 $iii，ai≤pia_i \le p_iai≤pi。$

特性 2 是指： $lcm⁡(pi)≤106\operatorname{lcm}(p_i) \le 10^6lcm(pi)≤106，即所有 pip_ipi 的最小公倍数不大于 10610^6106。$

对于所有的测试点， $\le 5T≤5，所有武器的攻击力 ≤106\le 10^6≤106，所有 pip_ipi 的最小公倍数 ≤1012\le 10^{12}≤1012。$

保证

【提示】

你所用到的中间结果可能很大，注意保存中间结果的变量类型。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<set>
#pragma GCC optimize(2)
using namespace std;
typedef long long LL;
const int maxn = 1e5 + 10 ;
LL a[maxn],p[maxn],atk[maxn],x,y,d;
multiset<LL> sw;

inline LL read()
{
    LL w=1,s=0; char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')w=-1; ch=getchar();}
    while(ch>='0' && ch<='9'){s=s*10+ch-'0'; ch=getchar();}
    return w*s;
}
void ext_gcd(LL a,LL b,LL &d,LL &x,LL &y)
{
    if (b==0)
    {
        d=a;x=1;y=0;
        return;
     }
     ext_gcd(b,a%b,d,y,x);  
     y-=a/b*x;
}

inline LL ksc(LL x,LL y,LL p){
    LL z=(long double)x/p*y;
    LL res=(unsigned long long)x*y-(unsigned long long)z*p;
    return (res+p)%p;
}

int main()
{
    LL T,n,m,k,mx,c;
    T=read();
    begin:while(T--){
        n=read();m=read();
        sw.clear();
        register multiset<LL>::iterator it;
        for(int i=1;i<=n;i++) a[i]=read();
        for(int i=1;i<=n;i++) p[i]=read();
        for(int i=1;i<=n;i++) atk[i]=read();
        for(int i=1;i<=m;i++)
            sw.insert(read());
        mx=c=0;m=1;
        for(int i=1;i<=n;i++){
            it=sw.upper_bound(a[i]);
            if(sw.begin()!=it) it--;
            k=*it;sw.erase(it);sw.insert(atk[i]);
            mx=max(mx,(a[i]-1)/k+1);//更新限制
            k%=p[i];a[i]%=p[i];
            if(!k && !p[i]) continue;
            ext_gcd(k,p[i],d,x,y);
            p[i]/=d;
            a[i]=ksc(a[i]/d,(x%p[i]+p[i])%p[i],p[i]);
            ext_gcd(m,p[i],d,x,y);
            if((a[i]-c)%d){
                puts("-1");goto begin;
            }
            m=m/d*p[i];
            c=(c+ksc( ksc(m/p[i],((a[i]-c)%m+m)%m,m) , (x%m+m)%m , m))%m;
        }
        if(c>=mx) printf("%lld\n",c);
        else       printf("%lld\n",c+m*((mx-c-1)/m+1));
    }
    return 0;
}

P4774 [NOI2018]屠龙勇士

题目描述

输入格式

输出格式

输入输出样例

说明/提示

猜你喜欢