[BZOJ4542][HNOI2016]大数(莫队)

4542: [Hnoi2016]大数

Time Limit: 20 Sec  Memory Limit: 128 MB
Submit: 1691  Solved: 589
[Submit][Status][Discuss]

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。

HINT

 2016.4.19新加数据一组

代码用时：1h

10W级的数据跑莫队比较正常吧。注意特判P==2 || P==5的情况（只需要判断个位数即可）。

WA了两次，struct D中x变量应该是long long的。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define rep(i,l,r) for (int i=l; i<=r; i++)
typedef long long ll;
using namespace std;

const int N=200100;
int n,m,x,y,cnt,bl[N],num[N],c[N];
ll P,ans,res[N];
char S[N];
struct D{ ll x; int y; }s[N];
struct T{ int l,r,id; }a[N];

bool cmp(D a,D b){ return a.x<b.x; }
bool cmp1(T a,T b){ return (bl[a.l]==bl[b.l]) ? a.r<b.r : bl[a.l]<bl[b.l]; }
void work(int x,int k){
    if (k==1) ans+=num[c[x]],num[c[x]]++;
            else num[c[x]]--,ans-=num[c[x]];
}

void solve(){
    rep(i,1,n){
        s[i]=s[i-1];
        if (!((S[i]-'0')%P)) s[i].x+=i,s[i].y++;
    }
    while (m--) scanf("%d%d",&x,&y),printf("%lld
",s[y].x-s[x-1].x-1ll*(x-1)*(s[y].y-s[x-1].y));
}

int main(){
    freopen("bzoj4542.in","r",stdin);
    freopen("bzoj4542.out","w",stdout);
    scanf("%lld",&P); scanf("%s%d",S+1,&m);
    n=strlen(S+1); ll now=1;
    if (P==2 || P==5) { solve(); return 0; }
    rep(i,1,m) scanf("%d%d",&a[i].l,&a[i].r),a[i].r++,a[i].id=i;
    for (int i=n; i; i--,now=now*10%P) s[i].x=(s[i+1].x+(S[i]-'0')*now)%P,s[i].y=i;
    s[++n].x=0; s[n].y=n; sort(s+1,s+n+1,cmp);
    rep(i,1,n){
        if (i==1 || s[i].x!=s[i-1].x) cnt++;
        c[s[i].y]=cnt;
    }
    int sz=(int)sqrt(n)+1;
    rep(i,1,n) bl[i]=(i-1)/sz;
    sort(a+1,a+m+1,cmp1);
    int L=1,R=1; num[c[1]]++;
    rep(i,1,m){
         while (R<a[i].r) work(++R,1);
         while (L>a[i].l) work(--L,1);
         while (R>a[i].r) work(R--,-1);
         while (L<a[i].l) work(L++,-1);
         res[a[i].id]=ans;
    }
    rep(i,1,m) printf("%lld
",res[i]);
    return 0;
}