P1494 [国家集训队]小Z的袜子

题目描述

作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……

具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双,甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬。

你的任务便是告诉小Z,他有多大的概率抽到两只颜色相同的袜子。当然,小Z希望这个概率尽量高,所以他可能会询问多个(L,R)以方便自己选择。

然而数据中有L=R的情况,请特判这种情况,输出0/1。

输入输出格式

输入格式:

输入文件第一行包含两个正整数N和M。N为袜子的数量,M为小Z所提的询问的数量。接下来一行包含N个正整数Ci,其中Ci表示第i只袜子的颜色,相同的颜色用相同的数字表示。再接下来M行,每行两个正整数L,R表示一个询问。

输出格式:

包含M行,对于每个询问在一行中输出分数A/B表示从该询问的区间[L,R]中随机抽出两只袜子颜色相同的概率。若该概率为0则输出0/1,否则输出的A/B必须为最简分数。(详见样例)

输入输出样例

输入样例#1: 复制

6 4
1 2 3 3 3 2
2 6
1 3
3 5
1 6
输出样例#1: 复制

2/5
0/1
1/1
4/15

说明

30%的数据中 N,M ≤ 5000;

60%的数据中 N,M ≤ 25000;

100%的数据中 N,M ≤ 50000,1 ≤ L < R ≤ N,Ci ≤ N。


Solution

今天认认真真理解了一下莫队(不带修!!),发现还是挺简单的??(也许是做的题太简单了.....)

核心要点是离线处理。通过分块排序调整时间复杂度降到$O(n\sqrt{n})$级别,在排序后的左右端点来回跳就好了。

在某些情况下像分块一样要预处理出一些东西(这道题还不需要)

最后再排回来序就好了。

然后这道题主要还是观察怎么更新答案比较重要......

对于$L,R$的询问。

设其中颜色为$x,y,z$的袜子的个数为$a,b,c$...

那么答案即为$(a*(a-1)/2+b*(b-1)/2+c*(c-1)/2....)/((R-L+1)*(R-L)/2)$

      $(C_a^2+C_b^2+C_c^2+....)/(C_{a+b+c+...}^2)$(选袜子的组合)

化简得:$(a^2+b^2+c^2+...x^2-(a+b+c+d+.....))/((R-L+1)*(R-L))$

即:$(a^2+b^2+c^2+...x^2-(R-L+1))/((R-L+1)*(R-L))$

所以只用每次更新时更新每种颜色数量的平方和就行了。

Code

#include<bits/stdc++.h>
#define LL long long
using namespace std; int n, m; int blo[], tot;
struct Node {
int l, r, id;
LL ans1, ans2;
} qus[];
bool cmp(Node a, Node b) { if(blo[a.l] == blo[b.l]) return a.r < b.r; return a.l < b.l; }
bool cmp2(Node a, Node b) { return a.id < b.id; } LL ans;
int cnt[];
void update(int x, int d) {
ans -= cnt[x] * cnt[x];
cnt[x] += d;
ans += cnt[x] * cnt[x];
} LL gcd(LL a, LL b) {
return b == ? a : gcd(b, a % b);
} int a[], ll[], rr[];
LL Ans[];
int main() {
scanf("%d%d", &n, &m);
for(int i = ; i <= n; i ++) scanf("%d", &a[i]);
for(int i = ; i <= m; i ++) {
scanf("%d%d", &ll[i], &rr[i]);
qus[i].l = ll[i], qus[i].r = rr[i];
qus[i].id = i;
}
tot = sqrt(n);
for(int i = ; i <= n; i ++) blo[i] = i / tot + ;
sort(qus + , qus + + m, cmp);
int L = qus[].l, R = qus[].r;
for(int i = L; i <= R; i ++)
update(a[i], );
qus[].ans1 = ans - (qus[].r - qus[].l + );
qus[].ans2 = 1ll * (qus[].r - qus[].l) * (qus[].r - qus[].l + );
if(qus[].l == qus[].r) qus[].ans1 = , qus[].ans2 = ;
for(int i = ; i <= m; i ++) {
if(qus[i].l == qus[i].r) {
qus[i].ans1 = , qus[i].ans2 = ; continue;
}
while(qus[i].l < L) {
L --;
update(a[L], );
}
while(qus[i].l > L) {
update(a[L], -);
L ++;
}
while(qus[i].r < R) {
update(a[R], -);
R --;
}
while(qus[i].r > R) {
R ++;
update(a[R], );
}
qus[i].ans1 = ans - (qus[i].r - qus[i].l + );
qus[i].ans2 = 1ll * (qus[i].r - qus[i].l) * (qus[i].r - qus[i].l + );
}
sort(qus + , qus + + m, cmp2);
for(int i = ; i <= m; i ++) {
LL ans1 = qus[i].ans1, ans2 = qus[i].ans2;
if(ans1 == ) ans2 = ;
else {
LL d = gcd(ans1, ans2);
ans1 /= d;
ans2 /= d;
}
printf("%lld/%lld\n", ans1, ans2);
}
return ;
}

【洛谷】1494:[国家集训队]小Z的袜子【莫队】的更多相关文章

  1. BZOJ 2038 [2009国家集训队]小Z的袜子 莫队

    2038: [2009国家集训队]小Z的袜子(hose) 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=2038 Descriptionw ...

  2. P1494 [国家集训队]小Z的袜子/莫队学习笔记(误

    P1494 [国家集训队]小Z的袜子 题目描述 作为一个生活散漫的人,小\(Z\)每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小\(Z\)再也无法忍受这恼人的找袜子过程,于是他 ...

  3. 洛谷 P1494 [国家集训队] 小Z的袜子

    题目概述: 小Z把N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双,甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬. 你的任务 ...

  4. BZOJ2038 [2009国家集训队]小Z的袜子 莫队+分块

    作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命…… 具体来说,小Z把这N只袜子从1到N编号,然后从 ...

  5. 洛谷P1494 [国家集训队]小Z的袜子

    Code: #include<cstdio> #include<iostream> #include<algorithm> #include<cstring& ...

  6. 洛谷 P1494 [国家集训队]小Z的袜子(莫队)

    题目链接:https://www.luogu.com.cn/problem/P1494 一道很经典的莫队模板题,然而每道莫队题的大体轮廓都差不多. 首先莫队是一种基于分块的算法,它的显著特点就是: 能 ...

  7. [日常摸鱼]bzoj2038[2009国家集训队]小Z的袜子-莫队算法

    今天来学了下莫队-这题应该就是这个算法的出处了 一篇别人的blog:https://www.cnblogs.com/Paul-Guderian/p/6933799.html 题意:一个序列,$m$次询 ...

  8. BZOJ2120/洛谷P1903 [国家集训队] 数颜色 [带修改莫队]

    BZOJ传送门:洛谷传送门 数颜色 题目描述 墨墨购买了一套N支彩色画笔(其中有些颜色可能相同),摆成一排,你需要回答墨墨的提问.墨墨会向你发布如下指令: 1. Q L R代表询问你从第L支画笔到第R ...

  9. BZOJ 2038: [2009国家集训队]小Z的袜子

    二次联通门 : BZOJ 2038: [2009国家集训队]小Z的袜子 /* BZOJ 2038: [2009国家集训队]小Z的袜子 莫队经典题 但是我并不认为此题适合入门.. Answer = ∑ ...

随机推荐

  1. JSON数据生成树——(四)

    1.页面中准备树的div <div class="user_left_tree_info"> <div class="user_left_tree_in ...

  2. MySQL修改datadir目录

    更改MySQL的数据目录由/usr/local/mysql/data/更换至/data/mysql/ 第一种方法修改my.cnf 1)关闭MySQL [root@test2 ~]# service m ...

  3. HDU 6061 RXD and functions

    题目链接:HDU-6061 题意:给定f(x),求f(x-A)各项系数. 思路:推导公式有如下结论: 然后用NTT解决即可. 代码: #include <set> #include < ...

  4. poj1102

    模拟 #include <iostream> #include <string> using namespace std; ][][] = { { ' ', '-', ' ', ...

  5. 移动端默认meta标签

    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><met ...

  6. win7下weblogic安装与部署项目调试记录

    下载 weblogic12c, 官网下载通用版本 安装 略 下载jdk   jdk所有版本连接 http://www.oracle.com/technetwork/java/archive-13921 ...

  7. MySQL学习笔记:三种组内排序方法

    由于MySQ没有提供像Oracle的dense_rank()或者row_number() over(partition by)等函数,来实现组内排序,想实现这个功能,还是得自己想想办法,最终通过创建行 ...

  8. Ubuntu Touch On Nexus4 Manual Install (手动安装) under Gentoo

    Table of Contents 1. 准备工作: 2. Saucy Salamander 3. 刷入 最新 版Touch 最近手里的 Nexus 4 手机一直闲置,它的配置要比我六年前买的笔记本还 ...

  9. CF 576A 猜数

    A给出一个数x,B每次猜一个y,A回答B,x是否可以被y整除,求出要猜的最小次数和需要猜的数. 枚举每个素数p,可以知道如果p^k<=n,则p^k一定需要选 Sample test(s)inpu ...

  10. windows svn post-commit 报错解决 error resolving case

    在svn仓库目录下有个hooks目录,下面建一个 post-commit.cmd 文件,有代码提交到仓库,自动checkout到指定目录.   @echo onSET REPOS=%1SET USER ...