bzoj4542: [Hnoi2016]大数（莫队）

　　这题...离散化...$N$和$n$搞错了...查了$2h$...QAQ

　　考虑$s[l...r]$，可以由两个后缀$suf[l]-suf[r+1]$得到$s[l...r]$代表的数乘$10^k$得到的结果，如果$p$不为$2$或$5$，即$gcd(p, 10^k)=1$，那么显然$s[l...r]$乘$10^k$模$p$为$0$的话，$s[l...r]$模p也为$0$，所以我们就可以变成询问$[l,r+1]$里有几个相同的后缀了。

　　如果$p$为$2$或$5$的话，我们还得判断这个数的个位是否是$2$或$5$的倍数。

#include
#include
#include
#include
#include
#include
#define ll long long
using namespace std;
const int maxn=500010, inf=1e9;
struct poi{int l, r, pos;}q[maxn];
int n, N, m, blo;
int bl[maxn], cnt[maxn], cnt2[maxn], b[maxn], sum[maxn];
ll p, ANS;
ll ans[maxn], mi[maxn];
char s[maxn];
inline void read(int &k)
{
    int f=1; k=0; char c=getchar();
    while(c<'0' || c>'9') c=='-'&&(f=-1), c=getchar();
    while(c<='9' && c>='0') k=k*10+c-'0', c=getchar();
    k*=f;    
} 
bool operator < (poi a, poi b)
{return bl[a.l]1)?a.rb.r));} 
inline void update(int x, int delta, int ty)
{
    if(delta==1)
    {
        ANS+=(!ty || p>10 || (s[x-1]-'0')%p==0)*delta*(ty?cnt[sum[x]]:cnt2[sum[x]]);
        cnt[sum[x]]+=delta;
        cnt2[sum[x]]+=(p>10 || (s[x-1]-'0')%p==0)*delta;
    }
    else
    {
        cnt[sum[x]]+=delta;
        cnt2[sum[x]]+=(p>10 || (s[x-1]-'0')%p==0)*delta;
        ANS+=(!ty || p>10 || (s[x-1]-'0')%p==0)*delta*(ty?cnt[sum[x]]:cnt2[sum[x]]);
    }
}
int main()
{
    scanf("%lld", &p); scanf("%s", s+1); n=strlen(s+1);
    blo=sqrt(n); for(int i=1;i<=n;i++) bl[i]=(i-1)/blo+1;
    mi[0]=1; for(int i=1;i<=n;i++) mi[i]=mi[i-1]*10%p;
    for(int i=n;i;i--) sum[i]=(sum[i+1]+mi[n-i]*(s[i]-'0'))%p, b[i]=sum[i]; N=n; b[++N]=0;
    sort(b+1, b+1+N); N=unique(b+1, b+1+N)-b-1;
    for(int i=1;i<=n+1;i++) sum[i]=lower_bound(b+1, b+1+N, sum[i])-b;
    read(m);
    for(int i=1;i<=m;i++) read(q[i].l), read(q[i].r), q[i].r++, q[i].pos=i;
    sort(q+1, q+1+m);
    for(int i=1, l=1, r=0;i<=m;i++)
    {
        while(l1, 0);
        while(l>q[i].l) update(--l, 1, 0);
        while(r1, 1);
        while(r>q[i].r) update(r--, -1, 1);
        ans[q[i].pos]=ANS;
    }
    for(int i=1;i<=m;i++) printf("%lld\n", ans[i]);
}