LOJ.2721.[NOI2018]屠龙勇士(扩展CRT 扩展欧几里得)
题目链接 LOJ
洛谷
rank前3无压力(话说rank1特判打表有意思么)
\(x*atk[i] - y*p[i] = hp[i]\)
对于每条龙可以求一个满足条件的\(x_0\),然后得到其通解\(x[i] = x_0 + p[i]/gcd*k\)。
怎么合并所有龙的通解?可以直接写成 \(Ans\equiv x_0(mod\ p[i]/gcd)\),用扩展中国剩余定理合并即可。
所有\(p[i]=1\)时要特判。(为啥呢...反正我知道它不对...)
所有\(p[i]=hp[i]\)时同余方程结果会是0,这显然不对。这种情况对于每条龙造成的伤害应该是p[i]的倍数(\(p[i]\mid x*atk[i]\))。算出击杀每条龙的最小的\(x\)求个\(lcm\)即可。
查了半上午发现一个问题。。
就算set里全是int,查XX_bound(long long)的话还是要用set<long long>
。。
好歹是把今年NOI两天的T1 A了(也没有多难)。不知明年会怎么样呢...
#include <set>
#include <cstdio>
#include <cctype>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define Failed return (void)puts("-1")
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define set_It std::multiset<LL>::iterator
typedef long long LL;
const int N=1e5+7;
int n,m,rew[N]/*reward*/;
LL hp[N],p[N],a[N],md[N];//Ans = a[i] (mod md[i])
std::multiset<LL> st;//!
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline LL readll()
{
LL now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline LL Mult(LL a,LL b,LL p)
{
LL tmp=a*b-(LL)((long double)a/p*b+1e-8)*p;
return tmp<0?tmp+p:tmp;
}
LL Gcd(LL a,LL b){
return b?Gcd(b,a%b):a;
}
void Exgcd(LL a,LL b,LL &g,LL &x,LL &y)
{
if(!b) g=a,x=1,y=0;
else Exgcd(b,a%b,g,y,x),y-=a/b*x;
}
void Spec1()
{
int ans=0; set_It it;
for(int i=1; i<=n; ++i)
{
it=st.upper_bound(hp[i]);
if(it!=st.begin()) --it;
if(!(*it)) Failed;
ans=std::max(ans,(int)((hp[i]+(*it)-1)/(*it)));//(int)ceil(1.0*hp[i]/(*it)));
st.erase(it), st.insert(rew[i]);
}
printf("%d\n",ans);
}
//void Spec2()
//{
// int ans=1; set_It it;
// for(int i=1,need; i<=n; ++i)
// {
// it=st.upper_bound(hp[i]);
// if(it!=st.begin()) --it;
// if(!(*it)) Failed;
// need=(hp[i]+(*it)-1)/(*it), ans=ans/Gcd(ans,need)*need;
// st.erase(it), st.insert(rew[i]);
// }
// printf("%d\n",ans);
//}
void Solve()
{
bool f1=1;//!
for(int i=1; i<=n; ++i) if(p[i]!=1) {f1=0; break;}
if(f1) {Spec1(); return;}
// f1=1;//然而数据没有这种情况...还是要有的
// for(int i=1; i<=n; ++i) if(p[i]!=hp[i]) {f1=0; break;}
// if(f1) {Spec2(); return;}
set_It it; int cnt=0;
for(int i=1; i<=n; ++i)
{
it=st.upper_bound(hp[i]);
if(it!=st.begin()) --it;
if(hp[i]&&!(*it)) Failed;
int atk=*it;
LL gcd,x0,y0,P;
Exgcd(atk,p[i],gcd,x0,y0);
if(hp[i]%gcd) Failed;
P=p[i]/gcd, x0=(x0%P+P)%P;
a[++cnt]=Mult(x0,hp[i]/gcd,P), md[cnt]=P;
st.erase(it), st.insert(rew[i]);
}
LL A=a[1],M=md[1],g,x,y,t,Mod;
for(int i=2; i<=cnt; ++i)
{
Exgcd(M,md[i],g,x,y);
if((a[i]-A)%g) Failed;
t=md[i]/g, x=Mult(x,(a[i]-A)/g,t), x=(x%t+t)%t;
Mod=M*t, A+=Mult(x,M,Mod)/*直接x*M会炸?...*/, A%=Mod, M=Mod;
}
printf("%lld\n",A);
}
int main()
{
freopen("dragon.in","r",stdin);
freopen("dragon.out","w",stdout);
for(int Case=read(); Case--; )
{
st.clear();
n=read(), m=read();
for(int i=1; i<=n; ++i) hp[i]=readll();
for(int i=1; i<=n; ++i) p[i]=readll();
for(int i=1; i<=n; ++i) rew[i]=read();
for(int i=1; i<=m; ++i) st.insert(read());
Solve();
}
return 0;
}
LOJ.2721.[NOI2018]屠龙勇士(扩展CRT 扩展欧几里得)的更多相关文章
- 洛谷 P4774 / loj 2721 [NOI2018] 屠龙勇士 题解【同余】【exgcd】【CRT】
推导过程存在漏洞+exCRT板子没打熟于是期望得分÷实际得分=∞? 题目描述 小 D 最近在网上发现了一款小游戏.游戏的规则如下: 游戏的目标是按照编号 \(1\sim n\) 顺序杀掉 \(n\ ...
- 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 每次用 ...
- P4774 [NOI2018]屠龙勇士
P4774 [NOI2018]屠龙勇士 先平衡树跑出打每条龙的atk t[] 然后每条龙有\(xt \equiv a[i](\text{mod }p[i])\) 就是\(xt+kp[i]=a[i]\) ...
- [洛谷P4774] [NOI2018]屠龙勇士
洛谷题目链接:[NOI2018]屠龙勇士 因为markdown复制过来有点炸格式,所以看题目请戳上面. 题解: 因为杀死一条龙的条件是在攻击\(x\)次,龙恢复\(y\)次血量\((y\in N^{* ...
- uoj396 [NOI2018]屠龙勇士
[NOI2018]屠龙勇士 描述 小 D 最近在网上发现了一款小游戏.游戏的规则如下: 游戏的目标是按照编号 1∼n 顺序杀掉 n 条巨龙,每条巨龙拥有一个初始的生命值 ai .同时每条巨龙拥有恢复能 ...
- NOI2018屠龙勇士(扩展CRT + splay(multiset))
QWQ 一到假期就颓废 哎 今年新鲜出炉的NOI题,QwQ同步赛的时候写的,后来交了一发洛谷,竟然过了 首先 根据题目,我们很容易得到,假设对应每一条龙的剑的攻击力是\(atk\)的话 \[a_i-x ...
- 洛谷 P4774 [NOI2018] 屠龙勇士
链接:P4774 前言: 交了18遍最后发现是多组数据没清空/ll 题意: 其实就是个扩中. 分析过程: 首先发现根据题目描述的选择剑的方式,每条龙对应的剑都是固定的,有查询前驱,后继(在该数不存在前 ...
- 扩展CRT +扩展LUCAS
再次感谢zyf2000超强的讲解. 扩展CRT其实就是爆推式子,然后一路合并,只是最后一个式子上我有点小疑惑,但整体还算好理解. #include<iostream> #include&l ...
随机推荐
- iperf测试网络带宽
http://blog.chinaaet.com/telantan/p/30901 https://boke.wsfnk.com/archives/288.html https://www.ibm.c ...
- 设计模式之Mixin模式
介绍 mixin模式就是一些提供能够被一个或者一组子类简单继承功能的类,意在重用其功能.在面向对象的语言中,我们会通过接口继承的方式来实现功能的复用.但是在javascript中,我们没办法通过接口继 ...
- Dubbo学习笔记4:服务消费端泛化调用与异步调用
本文借用dubbo.learn的Dubbo API方式来解释原理. 服务消费端泛化调用 前面我们讲解到,基于Spring和基于Dubbo API方式搭建简单的分布式系统时,服务消费端引入了一个SDK二 ...
- spring框架学习(三)spring与junit整合测试
package cn.mf.b_test; import javax.annotation.Resource; import org.junit.Test; import org.junit.runn ...
- 解决Tomcat6解压版在64位windows系统上无法启动服务的问题
解决Tomcat6解压版在64位windows系统上无法启动服务的问题 由于客户环境为64位windows系统,开发环境一直用32位.tomcat使用6.0.20非安装版.部署时发现在 ...
- 关于webpack下热更新?&自动刷新?的小记(非vue-cli)
写本随笔时:webpack4.6.0 为何标题用?号,因为老衲也不知是否用词正确,大概是这样的说法: webpack4.0引入生产模式和开发模式,在开发时使用 webpack 打包后不压缩,所以只需要 ...
- urb传输的代码分析【转】
转自:http://blog.csdn.net/zkami/article/details/2503829 urb传输的代码分析 如需引用,请注明出处blog.csdn.net/zkami 作者Zhe ...
- (mysql)触发器、事件、事务、函数
1.事务操作原理:事务开启之后Start transaction,所有的操作都会临时保存到事务日志.只有在得到commit才会关闭,否则清空:2.设置回滚点: savepoint 回滚点名字: 回到 ...
- Python_oldboy_自动化运维之路_paramiko,mysql(十二)
本节内容: paramiko mysql 1.paramiko http://www.cnblogs.com/wupeiqi/articles/5095821.html paramiko是一个模块,s ...
- Java封装概述
1.封装概述 private public 2.实现封装 例子: package com.java1995; public class Student { private String name; ...