2038: [2009国家集训队]小Z的袜子(hose)

2038: [2009国家集训队]小Z的袜子(hose)

Time Limit: 20 Sec  Memory Limit: 259 MB
Submit: 15879  Solved: 7205
[Submit][Status][Discuss]

Description

作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……
具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双,甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬。
你的任务便是告诉小Z,他有多大的概率抽到两只颜色相同的袜子。当然,小Z希望这个概率尽量高,所以他可能会询问多个(L,R)以方便自己选择。

Input

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

Output

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

Sample Input

6 4
1 2 3 3 3 2
2 6
1 3
3 5
1 6

Sample Output

2/5
0/1
1/1
4/15
【样例解释】
询问1:共C(5,2)=10种可能,其中抽出两个2有1种可能,抽出两个3有3种可能,概率为(1+3)/10=4/10=2/5。
询问2:共C(3,2)=3种可能,无法抽到颜色相同的袜子,概率为0/3=0/1。
询问3:共C(3,2)=3种可能,均为抽出两个3,概率为3/3=1/1。
注:上述C(a, b)表示组合数,组合数C(a, b)等价于在a个不同的物品中选取b个的选取方案数。
【数据规模和约定】
30%的数据中 N,M ≤ 5000;
60%的数据中 N,M ≤ 25000;
100%的数据中 N,M ≤ 50000,1 ≤ L < R ≤ N,Ci ≤ N。
代码如下:
 #include<bits/stdc++.h>
using namespace std;
#define IO ios::sync_with_stdio(false),cin.tie(0);
const int maxn = 1<<20;
int a[maxn];
int sz;
long long gcd(long long x,long long y){while(y>0){long long tmp=y;y=x%y;x=tmp;}return x;
}
long long comb(long long x){return x*(x-1)/2;
}
struct node{int l,r,id;bool operator < (const node& x) const {return (l-1)/sz==(x.l-1)/sz?r<x.r:(l-1)/sz<(x.l-1)/sz;}
}Q[maxn];
long long ans[maxn][2];
long long flag[maxn];int n,m,k;int L=1,R=0;
long long Ans=0;
void add(int x){Ans-=comb(flag[a[x]]);flag[a[x]]++;Ans+=comb(flag[a[x]]);
}
void del(int x){Ans-=comb(flag[a[x]]);flag[a[x]]--;Ans+=comb(flag[a[x]]);
}int main(){scanf("%d%d",&n,&m);sz=sqrt(n);for(int i=1;i<=n;i++){scanf("%d",&a[i]);}for(int i=1;i<=m;i++){scanf("%d%d",&Q[i].l,&Q[i].r);Q[i].id=i;}sort(Q+1,Q+1+m);flag[0]=1;for(int i=1;i<=m;i++){while(L<Q[i].l){del(L++);}while(L>Q[i].l){add(--L);}while(R<Q[i].r){add(++R);}while(R>Q[i].r){del(R--);}if(Q[i].l==Q[i].r){ans[Q[i].id][0]=0;ans[Q[i].id][1]=1;continue;}long long t=comb((long long)(Q[i].r-Q[i].l+1));long long gg=gcd(t,Ans);ans[Q[i].id][0]=Ans/gg;ans[Q[i].id][1]=t/gg; }for(int i=1;i<=m;i++){if(ans[i][0]==0){printf("0/1\n");}else printf("%lld/%lld\n",ans[i][0],ans[i][1]);}
}
View Code

我跑了3000ms

然后 某屁屁黄就嘲讽我 说他只用了800ms

此处贴上他的代码,也不是很懂他为什么就比我快?(男人那么快真的好吗,,,滴滴滴下车

代码如下:

#include <set>
#include <map>
#include <queue>
#include <stack>
#include <cmath>
#include <bitset>
#include <cstdio>
#include <string>
#include <vector>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;typedef long long ll;
typedef unsigned long long ull;#define bug printf("*********\n");
#define FIN freopen("in.txt", "r", stdin);
#define debug(x) cout<<"["<<x<<"]" <<endl;
#define IO ios::sync_with_stdio(false),cin.tie(0);const double eps = 1e-8;
const int mod = 1e9 + 7;
const int maxn = 5e4 + 7;
const double pi = acos(-1);
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;int n, m, block;
ll sum;
int a[maxn], cnt[maxn];struct node {int l, r, id;ll ans1, ans2;bool operator < (const node& x) const {return (l - 1) / block == (x.l - 1) / block ? r < x.r : (l - 1) / block < (x.l - 1) / block;}
}ask[maxn];void update (int p, int val) {sum -= cnt[p] * cnt[p];cnt[p] += val;sum += cnt[p] * cnt[p];
}ll gcd(ll a, ll b) {return b == 0 ? a : gcd(b, a % b);
}int main() {//FIN;while(~scanf("%d%d", &n, &m)) {block = sqrt(n);for(int i = 1; i <= n; i++) {scanf("%d", &a[i]);}for(int i = 1; i <= m; i++) {scanf("%d%d", &ask[i].l, &ask[i].r);ask[i].id = i;}sort(ask + 1, ask + m + 1);for(int i = 1, l = 1, r = 0; i <= m; i++) {for(; r < ask[i].r; r++) update(a[r + 1], 1);for(; r > ask[i].r; r--) update(a[r], -1);for(; l < ask[i].l; l++) update(a[l], -1);for(; l > ask[i].l; l--) update(a[l - 1], 1);if(ask[i].l == ask[i].r) {ask[ask[i].id].ans1 = 0, ask[ask[i].id].ans2 = 1;continue;}ll k1 = sum - (ask[i].r - ask[i].l + 1);ll k2 = (long long) (ask[i].r - ask[i].l + 1) * (ask[i].r - ask[i].l);ll gg = gcd(k1, k2);ask[ask[i].id].ans1 = k1 / gg;ask[ask[i].id].ans2 = k2 / gg;}for(int i = 1; i <= m; i++) {if(ask[i].ans1 == 0) {printf("0/1\n");} else {printf("%lld/%lld\n", ask[i].ans1, ask[i].ans2);}}}return 0;
}
View Code

 

转载于:https://www.cnblogs.com/buerdepepeqi/p/9452278.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/485421.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Jeff Dean三执笔:一文看尽2020年谷歌AI重大突破

来源&#xff1a;AI前线作者&#xff1a;Jeff Dean译者&#xff1a;核子可乐策划&#xff1a;钰莹 在二十年前刚刚加入谷歌时&#xff0c;我们关注的问题只有一个——如何面向这么多不同种类的联网计算机提供一整套质量出色且涵盖范围全面的网络信息搜索服务。到如今&#xff0…

输出参数

输出参数为resultType 1.简单类型&#xff08;8个基本类型String&#xff09; 2.输出参数为实体对象类型 3. 输出参数为实体对象类型的集合&#xff1a;虽然输出类型是集合&#xff0c;但是resultType依然写集合的元素类型 4. 输出参数类型为HashMap HashMap本身是一个集合…

模板 普通平衡树

Treap 版&#xff1a; 查询操作是最简单的&#xff0c;就是二叉搜索树。左旋右旋画一画就好了。 最近代码能力下降&#xff0c;码一码模板QAQ 1 #include<ctime>2 #include<cstdio>3 #include<cstring>4 #include<cstdlib>5 #include<algorithm>…

信号时频图的横纵坐标对调

博客1&#xff1a;如何翻转信号的时频图论文有时候需要绘制信号的时频图。Matlab有专用的函数spectrogram可以用来帮助我们绘制信号的时频图。默认的图例是横轴为频率&#xff0c;纵轴为时间。如果希望将横轴和纵轴对调&#xff0c;只需要在函数的末尾指定xaxis或yaxis即可。sp…

技术的顶点,是普惠人类还是超越人类

来源&#xff1a;大数据文摘我们很容易忽视我们在宏大的生态系统中的卓越地位&#xff0c;尤其是在一个容易依赖技术引领我们走向我们渴望的未来的时代。我们真的相信技术比人类更复杂和必要吗&#xff1f;被创造的东西真的能取代它的创造者吗&#xff1f;我们必须共同决定&…

HDU 2089 不要62

杭州人称那些傻乎乎粘嗒嗒的人为62&#xff08;音&#xff1a;laoer&#xff09;。 杭州交通管理局经常会扩充一些的士车牌照&#xff0c;新近出来一个好消息&#xff0c;以后上牌照&#xff0c;不再含有不吉利的数字了&#xff0c;这样一来&#xff0c;就可以消除个别的士司机…

Leetcode--149. 直线上最多的点数

给定一个二维平面&#xff0c;平面上有 n 个点&#xff0c;求最多有多少个点在同一条直线上。 示例 1: 输入: [[1,1],[2,2],[3,3]] 输出: 3 解释: ^ | | o | o | o -------------> 0 1 2 3 4 示例 2: 输入: [[1,1],[3,2],[5,3],[4,1],[2,3],[1,4]] 输出…

生态位构建-一种普遍的主动进化

来源&#xff1a;混沌巡洋舰人类正在塑造地球上生命进化的未来。我们不仅造成了大规模的物种灭绝&#xff0c;我们还迫使动物、植物和真菌适应我们的人造世界: 例如&#xff0c;城市鸟类现在唱得更高了&#xff0c;因为这种音调似乎有助于它们的歌声传递交通的声音。但是&#…

[引] 如何在Matlab中得到透明背景的图片?

% you created a figure and it is “current”. % the following, you could have guessed set(gcf,‘color’,‘none’); set(gca,‘color’,‘none’); % but this following little detail took me ages to figure out set(gcf,‘InvertHardCopy’,‘off’); 作者&#xf…

Openstack入门篇(十一)之neutron服务(控制节点)的部署与测试

1.Neutron的介绍 Neutron 为整个 OpenStack 环境提供网络支持&#xff0c;包括二层交换&#xff0c;三层路由&#xff0c;负载均衡&#xff0c;防火墙和 *** 等。Neutron 提供了一个灵活的框架&#xff0c;通过配置&#xff0c;无论是开源还是商业软件都可以被用来实现这些功能…

Leetcode--139. 单词拆分

给定一个非空字符串 s 和一个包含非空单词列表的字典 wordDict&#xff0c;判定 s 是否可以被空格拆分为一个或多个在字典中出现的单词。 说明&#xff1a; 拆分时可以重复使用字典中的单词。 你可以假设字典中没有重复的单词。 示例 1&#xff1a; 输入: s "leetcode…

【物理方程】物理学中最难的方程之一,你知道多少?

来源&#xff1a;数学中国物理学是一门包含许多方程式的学科&#xff0c;这些方程描述了从微观世界粒子的行为到宏观宇宙的演化。在所有的物理方程中&#xff0c;有一组在数学上被认为极具挑战性&#xff0c;还被克莱数学研究所列为七个“千禧年大奖问题”之一&#xff0c;它们…

SSL/TLS 配置

Quick Start 下列说明将使用变量名 $CATALINA_BASE 来表示多数相对路径所基于的基本目录。如果没有为 Tomcat 多个实例设置 CATALINA_BASE 目录&#xff0c;则 $CATALINA_BASE 就会设定为 $CATALINA_HOME 的值&#xff0c;也就是你安装 Tomcat 的目录。 在 Tomcat 中安装并配置…

在LaTex中引用图片失败,出现(??)怎么办

查了一些博客&#xff0c;发现是加图片或是表格标签时出的问题&#xff1a; 一般来说图片会命名&#xff1a;\caption{图片名}. 要想引用成功, 必须要把标签加在图片命名之后&#xff1a;即 \caption{Name} \label{Label_Name}

Leetcode--27. 移除元素

给定一个数组 nums 和一个值 val&#xff0c;你需要原地移除所有数值等于 val 的元素&#xff0c;返回移除后数组的新长度。 不要使用额外的数组空间&#xff0c;你必须在原地修改输入数组并在使用 O(1) 额外空间的条件下完成。 元素的顺序可以改变。你不需要考虑数组中超出新…

2021年10个激动人心的工程里程碑

来源&#xff1a;IEEE电气电子工程师未来智能实验室的主要工作包括&#xff1a;建立AI智能系统智商评测体系&#xff0c;开展世界人工智能智商评测&#xff1b;开展互联网&#xff08;城市&#xff09;云脑研究计划&#xff0c;构建互联网&#xff08;城市&#xff09;云脑技术…

设置 sideload Outlook Add-ins

上期&#xff0c;我们讲到了用前端技术去建立一个outlook add-ins 我们今天来讲解一下怎样测试一个sideload outlook add-ins。 1. 我们需要登录Outlook在Office 365中。 2. 选择右上角设置按钮&#xff0c;然后选择Manage add-ins选项 3. 打开manage add-in的页面之后&#xf…

用户体验已成过去时 AI要从公民乃至社会的角度去思考问题

来源&#xff1a;Venture Beat作者&#xff1a;Jarno M. Koponen编译&#xff1a;科技行者时至今日&#xff0c;科技巨头及其基于AI的数字平台与解决方案&#xff0c;完全有能力影响全球领导人、民族国家、跨国企业、全球股市乃至每位个人的命运。目前&#xff0c;各大主要数字…

Leetcode--523. 连续的子数组和

给定一个包含非负数的数组和一个目标整数 k&#xff0c;编写一个函数来判断该数组是否含有连续的子数组&#xff0c;其大小至少为 2&#xff0c;总和为 k 的倍数&#xff0c;即总和为 n*k&#xff0c;其中 n 也是一个整数。 示例 1: 输入: [23,2,4,6,7], k 6 输出: True 解释…

css比较特殊选择器汇总(持续更新)

1、css选择器中加号&#xff08;&#xff09;是啥意思&#xff1f; 加号&#xff08;&#xff09;为&#xff1a;相邻同胞选择器&#xff0c;如果需要选择紧接在另一个元素后的元素&#xff0c;而且二者有相同的父元素&#xff0c;可以使用相邻兄弟选择器(Adjacent sibling sel…