题目链接

大数除法是很麻烦的,考虑能不能将其条件化简

一段区间[l,r]|p,即num[l,r]|p,类似前缀,记后缀suf[i]表示[i,n]的这段区间代表的数字

于是有 suf[l]-suf[r+1]|p -> (suf[l]-suf[r+1])%p = 0 -> suf[l] ≡suf[r+1] (mod p)

即若suf[r+1]%p = suf[l]%p,则num[l,r]|p

于是我们可以把范围控制在p以内,查找是否有%p相等的区间 -> 莫队

即小Z的袜子

这样的实际意义是 \((suf[l]-suf[r+1])*10^{n-r}%p = 0\)

后面有个10的幂,如果n-r>0 对于p(p|10)显然无论[l,r]是什么都会满足条件

于是p=2,5时要换一种判断方式。显然若有一位A[i]满足A[i]|p,则之前的Aj都可以以它为结尾并对答案做出贡献

若数列(都加l-1)a,b,c,d,e中,A[a,d,e]|p,则答案为e+d+a-3*(l-1)

这个在线做就可以了 (另外有人知道这个我莫队为什么不对吗。。)

另外要注意p范围不定,要对出现的余数离散化

然后参数p和模数p会重着 要注意!

要注意会有A[n+1],离散化要注意这个变量默认0!

另外p<P时特判为什么不对。。

洛谷数据真的水。。在BZOJ被longlong什么的各种坑。

#include <cmath>
#include <cctype>
#include <cstdio>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=2e5+5,P=1e6; int n,m,size,tm[P+5],Cnt;
LL p/*LL!*/,A[N],suf[N],ref[N],Now,Ans[N],sum[N],num[N];
struct Ask
{
int l,r,id;
bool operator <(const Ask &a)const {
return l/size==a.l/size?r<a.r:l/size<a.l/size;
}
}q[N]; 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 void Update(int pos,int delta){
if(!(A[pos]%p)) Now+=delta*pos, Cnt+=delta;
}
void Spec()
{
for(int i=1; i<=n; ++i)
if(A[i]%p) sum[i]=sum[i-1],num[i]=num[i-1];
else sum[i]=sum[i-1]+i,num[i]=num[i-1]+1;
for(LL l,r,i=1; i<=m; ++i)
l=read(),r=read(), printf("%lld\n",sum[r]-sum[l-1]-(num[r]-num[l-1])*(l-1));
// for(int l=1,r=0,i=1; i<=m; ++i)
// {
// while(l<q[i].l) Update(l,-1),++l;
// while(l>q[i].l) --l,Update(l,1);
// while(r<q[i].r) ++r,Update(r,1);
// while(r>q[i].r) Update(r,-1),--r;
// Ans[q[i].id]=Now-Cnt*(q[i].l-1);
// }
}
inline void Add(int p){
if(!p) return;
Now+=tm[p], ++tm[p];
// if(++tm[p]>1) Now+=tm[p]-1;//不要判什么tm[]>1之类满足才更新答案,因为tm[]在中间过程是可以存在负数的
}
inline void Subd(int p){
if(!p) return;
--tm[p], Now-=tm[p];
// if(--tm[p]) Now-=tm[p];//这一写法效果一样(因为只判了!=0) 但不能这么写
}
void Normal()
{
for(int l=1,r=0,i=1; i<=m; ++i)
{
while(l<q[i].l) Subd(suf[l++]);
while(l>q[i].l) Add(suf[--l]);
while(r<q[i].r) Add(suf[++r]);
while(r>q[i].r) Subd(suf[r--]);
Ans[q[i].id]=Now;
}
}
int Find(int x,int r)
{
int l=1,mid;
while(l<r)
if(ref[mid=l+r>>1]>=x) r=mid;
else l=mid+1;
return l;
}
void Discrete()
{
for(int i=1; i<=n; ++i) ref[i]=suf[i];
std::sort(ref+1,ref+1+n);
int cnt=1;
for(int i=2; i<=n; ++i)
if(ref[i]!=ref[i-1]) ref[++cnt]=ref[i];
for(int i=1; i<=n; ++i) suf[i]=Find(suf[i],cnt);
suf[n+1]=ref[1]?0:1;
} int main()
{
scanf("%lld",&p);
char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);A[++n]=c-'0',c=gc());
m=read(), size=sqrt(n);
if(p==2||p==5) {Spec(); return 0;}
for(int i=1; i<=m; ++i)
q[i].l=read(),q[i].r=read()+1,q[i].id=i;//q:r+1
LL pw10=1;//longlong!
for(int i=n; i; --i)
suf[i]=(A[i]*pw10%p+suf[i+1])%p, pw10=pw10*10%p;
// if(p>P) Discrete();//这个if?
Discrete();
std::sort(q+1,q+1+m);
// if(p==2||p==5) Spec();
// else Normal();
Normal();
for(int i=1; i<=m; ++i) printf("%lld\n",Ans[i]); return 0;
}

BZOJ.4542.[HNOI2016]大数(莫队)的更多相关文章

  1. bzoj 4542: [Hnoi2016]大数 (莫队)

    Description 小 B 有一个很大的数 S,长度达到了 N 位:这个数可以看成是一个串,它可能有前导 0,例如00009312345.小B还有一个素数P.现在,小 B 提出了 M 个询问,每个 ...

  2. 【BZOJ4542】[Hnoi2016]大数 莫队

    [BZOJ4542][Hnoi2016]大数 Description 小 B 有一个很大的数 S,长度达到了 N 位:这个数可以看成是一个串,它可能有前导 0,例如00009312345.小B还有一个 ...

  3. BZOj 4540: [Hnoi2016]序列 [莫队 st表 预处理]

    4540: [Hnoi2016]序列 题意:询问区间所有子串的最小值的和 不强制在线当然上莫队啦 但是没想出来,因为不知道该维护当前区间的什么信息,维护前后缀最小值的话不好做 想到单调栈求一下,但是对 ...

  4. [BZOJ4542] [Hnoi2016] 大数 (莫队)

    Description 小 B 有一个很大的数 S,长度达到了 N 位:这个数可以看成是一个串,它可能有前导 0,例如00009312345.小B还有一个素数P.现在,小 B 提出了 M 个询问,每个 ...

  5. BZOJ.4540.[HNOI2016]序列(莫队/前缀和/线段树 单调栈 RMQ)

    BZOJ 洛谷 ST表的一二维顺序一定要改过来. 改了就rank1了哈哈哈哈.自带小常数没办法. \(Description\) 给定长为\(n\)的序列\(A_i\).\(q\)次询问,每次给定\( ...

  6. BZOJ 4540 [Hnoi2016]序列 | 莫队 详细题解

    传送门 BZOJ 4540 题解 --怎么说呢--本来想写线段树+矩阵乘法的-- --但是嘛--yali的机房太热了--困--写不出来-- 于是弃疗,写起了莫队.(但是我连莫队都想不出来!) 首先用单 ...

  7. bzoj4542 [Hnoi2016]大数 莫队+同余

    题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4542 题解 我们令 \(f_i\) 表示从 \(i\) 到 \(n\) 位组成的数 \(\bm ...

  8. bzoj 4542: [Hnoi2016]大数

    Description 小 B 有一个很大的数 S,长度达到了 N 位:这个数可以看成是一个串,它可能有前导 0,例如00009312345 小B还有一个素数P.现在,小 B 提出了 M 个询问,每个 ...

  9. 洛谷P3245 [HNOI2016]大数(莫队)

    题意 题目链接 Sol 莫队板子题.. 维护出每个位置开始的字符串\(mod P\)的结果,记为\(S_i\) 两个位置\(l, r\)满足条件当且仅当\(S_l - S_r = 0\),也就是\(S ...

随机推荐

  1. springboot系列七:springboot 集成 MyBatis、事物配置及使用、druid 数据源、druid 监控使用

    一.MyBatis和druid简介 MyBatis 是一款优秀的持久层框架,它支持定制化 SQL.存储过程以及高级映射.MyBatis 避免了几乎所有的 JDBC 代码和手动设置参数以及获取结果集.M ...

  2. tomcat下部署应用helloworld

    部署应用(简单)1.到Tomcat的安装目录的webapps目录,可以看到ROOT,examples, tomcat-docs之类Tomcat自带的的目录.2.在webapps目录下新建一个目录mya ...

  3. [学习笔记]JS 数组Array push相关问题

    前言: 今天用写了一个二维数组,都赋值为零,然后更新其中一个值,结果和预期是不一样,会整列的相同位置都是同一个值. 1.用Chrome的控制台样例如下: arrs[2][2] =1的赋值,竟然是三个数 ...

  4. es6 entries(),keys() 和 values()

    for (let index of ['a', 'b'].keys()) { console.log(index); } // 0 for (let elem of ['a', 'b'].values ...

  5. 自适应电脑、手机和iPad的网页设计方法

    随着3G的普及,越来越多的人使用手机上网. 移动设备正超过桌面设备,成为访问互联网的最常见终端.于是,网页设计师不得不面对一个难题:如何才能在不同大小的设备上呈现同样的网页? 手机的屏幕比较小,宽度通 ...

  6. 并发之atomicInteger与CAS机制

    并发之atomic与CAS自旋锁 通过前几章的讲解我们知道i++这种类似操作是不安全的.针对这种情况,我们可能会想到利用synchronize关键字实现线程同步,保证++操作的原子性,的确这是一种有效 ...

  7. hdu2476

    /* dp[l][r]表示将任意串的[l,r]刷成s2样子的最小代价 ans[i]表示将s1的前i位刷成s2的代价 按照区间dp的常用做法,dp[l][r]的状态由dp[l][k],dp[k+1][r ...

  8. Hibrenate关系映射(一对一外键关联)

    一.一对一(单向):使用外部索引将其中的一个类作为parent,相对应的一个就是子类,并且参照父 类的主键ID来生成数据库表.(比如:可以将husband中设置一个wife_id对应wife中的主键i ...

  9. eclipse的操作

    IDEA至少在4G内存的电脑才能使用 eclipse中:项目名字小写 close project:关掉项目 删除未尽的项目导入eclipse中的步骤: 左边右键>>>import&l ...

  10. MyBatis查询,返回值Map或List<Map>

    一.返回值Map 1.mapper.xml <select id="selectUserMapLimitOne" resultType="java.util.Has ...