线性莫队算法 BZOJ 2038 小Z的袜子

        又是一种能拿到满分的乱搞算法……

        对于一个序列，假设[L,R]转移到[L+1,R],[L,R-1],[L-1,R],[L,R+1]四种状态的时间复杂度不是很大（一般小于log n什么的），便可以使用莫队

        对于本题来说，对于区间[l,r],ans = sigma(num[k]^2/(r-l+1)*(r-l))，很显然转移时间复杂度为O（1）

        于是就可以乱搞了~

        说一下具体实现，将序列分块，对于一个查询[Li,Ri]，优先级：Li所在块的序号，Ri所在位置，也就是

    return l/block_size < rhs.l/block_size || (l/block_size == rhs.l/block_size && r < rhs.r);

        均摊不是很会算，不过YY一下也知道时间复杂度应该是O(n*sqrt(n))的

#include <cstdio>
#include <cmath>
#include <algorithm>

using namespace std;

const int MAXN = 50001;
const int MAXM = 50001;
const int MAXNUM = 50001;

int block_size;

struct point
{
	int l,r,id;
	bool operator < (const point &rhs) const
	{
		return l/block_size < rhs.l/block_size || (l/block_size == rhs.l/block_size && r < rhs.r);
	}
};

int n,m;
long long data[MAXN];
long long ansA[MAXM],ansB[MAXN];
point P[MAXM];
long long num[MAXNUM];

long long gcd(long long a,long long b)
{
	return b ? gcd(b,a%b) : a;
}

int main()
{
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++) scanf("%lld",&data[i]);
	for (int i=1;i<=m;i++) scanf("%d%d",&P[i].l,&P[i].r),P[i].id = i;
	block_size = ceil(sqrt(n));
	sort(P+1,P+m+1);
	int l = 1, r = 0;
	long long temp = 0;
	for (int i=1;i<=m;i++)
	{
		while (r < P[i].r)
		{
			r++;
			temp -= num[data[r]]*num[data[r]];
			num[data[r]]++;
			temp += num[data[r]]*num[data[r]];
		}
		while (r > P[i].r)
		{
			temp -= num[data[r]]*num[data[r]];
			num[data[r]]--;
			temp += num[data[r]]*num[data[r]];
			r--;
		}
		while (l < P[i].l)
		{
       		temp -= num[data[l]]*num[data[l]];
			num[data[l]]--;
			temp += num[data[l]]*num[data[l]];
			l++;
		}
		while (l > P[i].l)
		{
			l--;
			temp -= num[data[l]]*num[data[l]];
			num[data[l]]++;
			temp += num[data[l]]*num[data[l]];
		}
		ansA[P[i].id] = temp-(r-l+1);
		ansB[P[i].id] = (long long) (r-l+1)*(r-l);
		long long d = gcd(ansA[P[i].id],ansB[P[i].id]);
		ansA[P[i].id] /= d; ansB[P[i].id] /= d;
	}
	for (int i=1;i<=m;i++) printf("%lld/%lld\n",ansA[i],ansB[i]);
	return 0;
}