[2017纪中10-27]查询 离散化+记忆询问+暴力剪枝

题面
其实这题是O(n^2)的。。。
首先离散化，那么询问最多也就n^2种，用一个数组记录每种询问的答案，如果重复询问就O(1)回答。
设f[i][c]表示1~i中出现过多少次c，那么区间[l,r]满足x和y出现次数相同当且仅当f[r][x]-f[l-1][x]=f[r][y]-f[l-1][y]，移项得f[r][x]-f[r][y]=f[l-1][x]-f[l-1][y]。即维护f[i][x]-f[i][y]有多少种取值。暴力扫过去显然是O(nq)的，但我们考虑假如一段区间x,y都没有出现过，那么这一段的f[i][x]-f[i][y]取值都相同，我们直接加一个区间长度，再跳过去就是了。
对于每一种c，设Tc表示c总共出现的次数，就算c和其他所有数都询问一次，复杂度最多Tc *(n-Tc)=Tc *n-Tc^2。总复杂度=∑(Tc *n-Tc^2)<∑ Tc *n=n^2。
代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int z=8010;
int n,q,a[z],c[z],cnt[z<<1],ans[z][z],lst[z],nxt[z],fir[z],tim[z]; 
int read()
{
    int x=0;
    char ch=getchar();
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9') {x=x*10+ch-'0';ch=getchar();}
    return x;
}
int main()
{
    n=read();
    q=read();
    for(int i=1;i<=n;i++) a[i]=read(),c[i]=a[i];
    sort(c+1,c+n+1);
    for(int i=1;i<=n;i++) a[i]=lower_bound(c+1,c+n+1,a[i])-c;
    memset(ans,-1,sizeof(ans)); 
    for(int i=0;i<=n;i++)
        ans[i][i]=n*(n+1)/2;
    for(int i=1;i<=n;i++)
    {

        tim[a[i]]++;
        if(lst[a[i]]==0) fir[a[i]]=i;
        else nxt[lst[a[i]]]=i;
        lst[a[i]]=i;
    }
    for(int i=1;i<=n;i++)
        nxt[lst[i]]=n+1;

    for(int i=1;i<=q;i++)   
    {
        int x=read(),y=read(),px,py,sz;
        px=lower_bound(c+1,c+n+1,x)-c;if(c[px]!=x) x=0;else x=px;
        py=lower_bound(c+1,c+n+1,y)-c;if(c[py]!=y) y=0;else y=py;
        if(ans[x][y]!=-1) {printf("%d\n",ans[x][y]);continue;}
        ans[x][y]=0;
        sz=max(tim[x],tim[y]);
        for(int i=z-sz;i<=z+sz;i++) cnt[i]=0;
        int p=fir[x],q=fir[y],r=0,nx=0,ny=0;
        while(p<=n||q<=n)
        {
            if(p<q)
            {
                cnt[z+nx-ny]+=p-r;
                nx++;
                r=p;
                p=nxt[p];
            }
            else
            {
                cnt[z+nx-ny]+=q-r;
                ny++;
                r=q;
                q=nxt[q];
            }
        }
        cnt[z+nx-ny]+=n-r+1;
        for(int i=z-sz;i<=z+sz;i++) ans[x][y]+=cnt[i]*(cnt[i]-1)/2;
        ans[y][x]=ans[x][y];
        printf("%d\n",ans[x][y]);
    }
    return 0;
}