HDU5145:5145 ( NPY and girls ) (莫队算法+排列组合+逆元)
传送门
题意
给出n个数,m次访问,每次询问[L,R]的数有多少种排列
分析
\(n,m<=30000\),我们采用莫队算法,关键在于区间如何\(O(1)\)转移,由排列组合知识得到,如果加入一个数,\(区间值*区间长度/该数出现次数\),减去一个数则相反操作讲解
trick
1.我的原先莫队写法不能ac,原因是我传入的L,R是全局变量,在insert和erase前就++,--了,而没有达到预期目的,把R++,L--分离即可
代码
//wa
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define F(i,a,b) for(int i=a;i<=b;++i)
#define R(i,a,b) for(int i=a;i<b;++i)
#define mem(a,b) memset(a,b,sizeof(a))
const ll mod = 1e9+7;
int t,n,m,len,L,R;
int a[30030],cnt[30030];
ll ans[30030],ret;
struct node
{
int l,r,id,block;
node(){}
node(int l,int r,int id):l(l),r(r),id(id){block=l/len;}
bool operator<(const node &p)const
{
return block==p.block?r<p.r:block<p.block;
}
}e[30030];
ll inv[30030];
void get_inverse(int n, ll p) {
inv[1] = 1;
for (int i = 2; i <= n; ++i) {
inv[i] = (p - p / i) * inv[p % i] % p;
}
}
void insert(int loc)
{
cnt[a[loc]]++;
//printf("l=%d r=%d\n",L,R);
ret=ret*(R-L+1)%mod;
ret=ret*inv[cnt[a[loc]]]%mod;
//printf("ret=%lld\n",ret);
}
void erase(int loc)
{
ret=ret*cnt[a[loc]]%mod;
ret=ret*inv[R-L+1]%mod;
cnt[a[loc]]--;
}
int main()
{
get_inverse(30000,mod);
for(scanf("%d",&t);t--;)
{
scanf("%d %d",&n,&m);
F(i,1,n) scanf("%d",a+i);
len=sqrt(n);
F(i,1,m)
{
int left,right;
scanf("%d %d",&left,&right);
e[i]=node(left,right,i);
}
sort(e+1,e+1+m);
L=1,R=0;
ret=1;
mem(ans,0);mem(cnt,0);
//F(i,1,n) printf("%lld%c",inv[i],i==n?'\n':' ');
F(i,1,m)
{
//printf("ret1=%lld\n",ret);
while(R<e[i].r) insert(++R);
while(L>e[i].l) insert(--L);
while(R>e[i].r) erase(R--);
while(L<e[i].l) erase(L++);
ans[e[i].id]=ret;
//printf("ret2=%lld\n",ret);
}
F(i,1,m) printf("%lld\n",ans[i]);
}
return 0;
}
//ac
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define F(i,a,b) for(int i=a;i<=b;++i)
#define R(i,a,b) for(int i=a;i<b;++i)
#define mem(a,b) memset(a,b,sizeof(a))
const ll mod = 1000000007;
int t,n,m,len,L,R;
int a[30030],cnt[30030];
ll ans[30030],inv[30030];
ll ret;
struct node
{
int l,r,id,block;
node(){}
node(int _l,int _r,int _id):l(_l),r(_r),id(_id){block=l/len;}
bool operator<(const node &p)const
{
return block==p.block?r<p.r:block<p.block;
}
}e[30030];
ll pow_mod(ll a,ll p)
{
ll ans=1;
for(ll ret=a;p;p>>=1,ret=ret*ret%mod) if(p&1) ans=ans*ret%mod;
return ans%mod;
}
void insert(int loc)
{
cnt[a[loc]]++;
//printf("l=%d r=%d\n",L,R);
ret=ret*(R-L+1)%mod;
ret=ret*inv[cnt[a[loc]]]%mod;
//printf("ret=%lld\n",ret);
}
void erase(int loc)
{
ret=ret*cnt[a[loc]]%mod;
ret=ret*inv[R-L+1]%mod;
cnt[a[loc]]--;
}
int main()
{
F(i,1,30000) inv[i]=pow_mod(i,mod-2);
for(scanf("%d",&t);t--;)
{
scanf("%d %d",&n,&m);
F(i,1,n) scanf("%d",a+i);
len=sqrt(n);
F(i,1,m)
{
int left,right;
scanf("%d %d",&left,&right);
e[i]=node(left,right,i);
}
sort(e+1,e+1+m);
ret=1;
mem(ans,0);mem(cnt,0);
L=1;R=0;
//F(i,1,n) printf("%lld%c",inv[i],i==n?'\n':' ');
F(i,1,m)
{
//printf("ret1=%lld\n",ret);
while(R<e[i].r) {++R;insert(R);}
while(L>e[i].l) {--L;insert(L);}
while(R>e[i].r) {erase(R);R--;}
while(L<e[i].l) {erase(L);L++;}
ans[e[i].id]=ret;
//printf("ret2=%lld\n",ret);
}
F(i,1,m) printf("%lld\n",ans[i]);
}
return 0;
}
HDU5145:5145 ( NPY and girls ) (莫队算法+排列组合+逆元)的更多相关文章
- HDU 5145 NPY and girls 莫队+逆元
NPY and girls Problem Description NPY's girlfriend blew him out!His honey doesn't love him any more! ...
- NPY and girls-HDU5145莫队算法
Time Limit: 8000/4000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Problem Description ...
- hdu_5145_NPY and girls(莫队算法+组合)
题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=5145 题意:给你n,m,共有n个女孩,标号为1—n,n个数xi表示第ith个女孩在第xi个教室,然后下 ...
- HDU 5145 - NPY and girls
题意: cases T(1≤T≤10) (0<n,m≤30000) (0<ai≤30000) n个数ai 表示n个女孩所在教室 m次询问 [L,R](1 <= L <= ...
- HDU 5145 NPY and girls(莫队算法+乘法逆元)
[题目链接] http://acm.hdu.edu.cn/showproblem.php?pid=5145 [题目大意] 给出一个数列,每次求一个区间数字的非重排列数量.答案对1e9+7取模. [题解 ...
- hdu 5145(莫队算法+逆元)
NPY and girls Time Limit: 8000/4000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Tot ...
- NBUT 1457 莫队算法 离散化
Sona Time Limit:5000MS Memory Limit:65535KB 64bit IO Format: Submit Status Practice NBUT 145 ...
- BZOJ 2038: [2009国家集训队]小Z的袜子(hose) [莫队算法]【学习笔记】
2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 7687 Solved: 3516[Subm ...
- Codeforces617 E . XOR and Favorite Number(莫队算法)
XOR and Favorite Number time limit per test: 4 seconds memory limit per test: 256 megabytes input: s ...
随机推荐
- C# 实体类序列化与反序列化一 (XmlSerializer)
/// <summary> /// 实体类序列化的反序列化的类 /// </summary> /// <typeparam name="T">& ...
- 【自用】OI计划安排表一轮
网络流√ 上下界最大流√ 线性规划转费用流√ RMQ优化建图√ 单纯形√ 字符串相关 hash√ 扩展KMP 回文自己主动机 数据结构 平衡树 启示式合并 替罪羊树 LCT 树套树 KD-Tree 二 ...
- leetcode Valid Palindrome C++&python 题解
题目描写叙述 Given a string, determine if it is a palindrome, considering only alphanumeric characters and ...
- VS2012,VS2010无法生成dll程序集的解决办法
在我们做项目的时候总会遇到dll程序集无法生成导致各种问题. 通常我们的做法就是清理项目,然后重新生成,或者直接到bin目录下删除所有dll,然后重新生成. 有时候某几个dll就是不生成. 这时候就需 ...
- Using Virtual Serial Ports on Linux (Ubuntu)
http://www.xappsoftware.com/wordpress/2013/10/07/using-virtual-serial-ports-on-linux-ubuntu/?goback= ...
- UWP 新手教程2——怎样实现自适应用户界面
系列文章 UWP新手教程1--UWP的前世今生 如上文所说的,布局面板依据可用的屏幕空间.指定界面元素的大小和位置. 比如StackPanel 会水平或垂直排列界面元素.Grid 布局与CSS 中的表 ...
- 嵌入式开发之davinci--- 8168 question about capture PAL on 8168
http://e2e.ti.com/support/dsp/davinci_digital_media_processors/f/717/t/340483.aspx?pi199607=2
- Boost下载安装编译配置使用指南(含Windows和Linux) .
理论上,本文适用于boost的各个版本,尤其是最新版本1.39.0:适用于各种C++编译器,如VC6.0(部分库不支持),VS2003,VS2005,VS2008,gcc,C++ Builder等.先 ...
- base64 hash sha
/*! * Crypto-JS v1.1.0 * http://code.google.com/p/crypto-js/ * Copyright (c) 2009, Jeff Mott. All ri ...
- 20170228 Z_po_send_email
FUNCTION zmm_po_send_email. function zmm_po_send_email. *"------------------------------------- ...