HNOI2016(BZOJ4542) 大数

HNOI2016 Day2 T3 大数

 

Description

  小 B 有一个很大的数 S,长度达到了 N 位;这个数可以看成是一个串,它可能有前导 0,例如00009312345
。小B还有一个素数P。现在,小 B 提出了 M 个询问,每个询问求 S 的一个子串中有多少子串是 P 的倍数(0 也
是P 的倍数)。例如 S为0077时,其子串 007有6个子串:0,0,7,00,07,007;显然0077的子串007有6个子串都是素
数7的倍数。

Input

  第一行一个整数:P。第二行一个串:S。第三行一个整数:M。接下来M行,每行两个整数 fr,to,表示对S 的
子串S[fr…to]的一次询问。注意:S的最左端的数字的位置序号为 1;例如S为213567,则S[1]为 2,S[1…3]为 2
13。N,M<=100000,P为素数

Output

  输出M行,每行一个整数,第 i行是第 i个询问的答案。

Sample Input

11
121121
3
1 6
1 5
1 4

Sample Output

5
3
2
//第一个询问问的是整个串,满足条件的子串分别有:121121,2112,11,121,121。

 

 

 

正解:莫队算法

解题报告:

   大概题意:给一串数和一个质数,然后每次询问一段区间,问对于这一个区间,有多少个子串组成的大数能够整除这个质数  

  今天为了做这道题先去学了莫队算法,然后A掉了莫队算法的入门题目——小Z的袜子。

  考场上面我傻逼的打了一个高精度,华丽丢掉暴力分。然而我发现只需要取个模就可以了,考场上傻了。

  学完莫队算法之后,发现这道题其实就是一个裸题。

 

  一开始依然是莫队算法的方式,按左端点所在块的编号为第一关键字,右端点编号为第二关键字排序,

  考虑先暴力处理每个块的第一组询问,之后可以发现我们只需要微调一下区间的左右端点就可以了。比如说:我上次处理了1 5,那么1 6就可以只把6加进去就可以了;如果这      一次是2 6,那么我还需要把1给去掉。

  考虑怎么转移:先离散化一下余数,记录[l, r]中这个数出现了几次,区间长度±1时显然答案改变值为这个点的余数的原出现次数(假设第一个数余数为2,用cnt[2]表示余数为2的数在当前区间的出现次数,那么去掉1之后答案改变值为cnt[2],自己yy一下应该想得通的)。

     然后好像很有道理了,AC了?!好吧,其实是善良的出题人没有卡我们,其实p = 2 或 p = 5的情况是可以被卡掉的

  特判一下就好了,两个数组分别表示[1, i]中2或5的倍数时有多少种情况及有多少个数末尾是2或5的倍数,用前缀和维护。

  这样就可以完美AC了,可以过掉BZOJ上的新增数据了

  下面这份是没改过的,不要怪我懒QAQ 

 

//It is made by ljh2000
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <string>
using namespace std;
typedef long long LL;
const int MAXN = 200011;
int MOD,m,len,l,r,block,L;
LL ans,A[MAXN],cnt[MAXN],a[MAXN],c[MAXN];
char ch[MAXN];
struct ask{int l,r,belong,id;}Q[MAXN];
inline bool cmp(ask q,ask qq){ if(q.belong==qq.belong) return q.r<qq.r; return q.belong<qq.belong; }
inline int getint(){
    int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
    if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
}

inline void change(int val,int type){
	ans-=cnt[val]*(cnt[val]-1)/2;
	cnt[val]+=type;
	ans+=cnt[val]*(cnt[val]-1)/2;
}

inline void work(){
	MOD=getint(); scanf("%s",ch+1);	len=strlen(ch+1); LL/*!!!*/ now=0,ss=1; now=ch[len]-'0'; a[len]=now%MOD; c[len]=a[len];
	block=sqrt(len); for(int i=len-1;i>=1;i--) { ss*=10; ss%=MOD; now+=ss*(ch[i]-'0'); now%=MOD; a[i]=now; c[i]=a[i]; }
	sort(c+1,c+len+1); L=unique(c+1,c+len+1)-c-1; for(int i=1;i<=len;i++) a[i]=lower_bound(c+1,c+len+1,a[i])-c; 
	if(c[1]==0) a[len+1]=1; else a[len+1]=0;/*!!!*/
	m=getint(); for(int i=1;i<=m;i++) Q[i].l=getint(),Q[i].r=getint(),Q[i].belong=(Q[i].l-1)/block+1,Q[i].id=i;
	sort(Q+1,Q+m+1,cmp); l=1; r=0; 
	for(int i=1;i<=m;i++) {
		while(r<=Q[i].r) r++,change(a[r],1);
		while(r>Q[i].r+1) change(a[r],-1),r--;
		while(l<Q[i].l) change(a[l],-1),l++;
		while(l>Q[i].l) l--,change(a[l],1);
		A[Q[i].id]=ans;
	}
	for(int i=1;i<=m;i++) printf("%lld\n",A[i]);
}

int main()
{
    work();
    return 0;
}

 

  

 

本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ 转载请注明出处,侵权必究,保留最终解释权!
posted @ 2016-05-01 23:05  ljh_2000  阅读(409)  评论(0编辑  收藏  举报