传送门

题意

给出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 ) (莫队算法+排列组合+逆元)的更多相关文章

  1. HDU 5145 NPY and girls 莫队+逆元

    NPY and girls Problem Description NPY's girlfriend blew him out!His honey doesn't love him any more! ...

  2. NPY and girls-HDU5145莫队算法

    Time Limit: 8000/4000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Problem Description ...

  3. hdu_5145_NPY and girls(莫队算法+组合)

    题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=5145 题意:给你n,m,共有n个女孩,标号为1—n,n个数xi表示第ith个女孩在第xi个教室,然后下 ...

  4. 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 <= ...

  5. HDU 5145 NPY and girls(莫队算法+乘法逆元)

    [题目链接] http://acm.hdu.edu.cn/showproblem.php?pid=5145 [题目大意] 给出一个数列,每次求一个区间数字的非重排列数量.答案对1e9+7取模. [题解 ...

  6. hdu 5145(莫队算法+逆元)

    NPY and girls Time Limit: 8000/4000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Tot ...

  7. NBUT 1457 莫队算法 离散化

    Sona Time Limit:5000MS     Memory Limit:65535KB     64bit IO Format: Submit Status Practice NBUT 145 ...

  8. BZOJ 2038: [2009国家集训队]小Z的袜子(hose) [莫队算法]【学习笔记】

    2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec  Memory Limit: 259 MBSubmit: 7687  Solved: 3516[Subm ...

  9. Codeforces617 E . XOR and Favorite Number(莫队算法)

    XOR and Favorite Number time limit per test: 4 seconds memory limit per test: 256 megabytes input: s ...

随机推荐

  1. 【整理】nand相关

    记录nand相关知识.主要是mtd和ubi 什么是UBI 它是一种flash管理方式 flash是一系列连续的物理擦除块组成的. UBI卷是一系列连续的逻辑擦除块(eraseblock),每一块都能够 ...

  2. SpringMVC框架下使用jfreechart绘制折线图,柱状图,饼状图

    java代码 @Controller public class CityAction { @Autowired private CityBiz cityBiz; //柱状图 @RequestMappi ...

  3. 云打印-Beta-凡事预则立

    凡事预则立 课程名称:软件工程1916|W(福州大学) 团队名称: 云打印 作业要求: 项目Beta冲刺(团队) 作业目标:Beta冲刺 团队队员 队员学号 队员姓名 个人博客地址 备注 221600 ...

  4. Python开发【2.1 面向对象】

    1.面向对象概述 类(Class): 用来描述具有相同的属性和方法的对象的集合.它定义了该集合中每个对象所共有的属性和方法.对象是类的实例. 类变量:类变量在整个实例化的对象中是公用的.类变量定义在类 ...

  5. Kafka知识点汇总

    整体结构 watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvZXJpY19zdW5haA==/font/5a6L5L2T/fontsize/400/fill/I ...

  6. MUI-折叠面板效果accordion

    在做开发的过程中我们经经常使用到折叠面板. 那我们来看下折叠面板到底是怎么使用. 废话不多说. 代码粘下来: <!DOCTYPE html> <html> <head&g ...

  7. XFire WebService demo

    XFire创建WebService实例应用 XFire使得在JavaEE应用中发布Web服务变得轻而易举.和其他Web服务引擎相比,  XFire的配置非常简单,可以非常容易地和Spring集成.   ...

  8. C# 给窗体添加事件

    1.https://zhidao.baidu.com/question/588485101.html

  9. swing_AbstractTableModel 创建表格

    import javax.swing.table.AbstractTableModel; public class MyTable extends AbstractTableModel { /** * ...

  10. Java 8 新的时间日期 API

    1. 概述 1.1 简介 Java 8 引入了一套全新的时间日期API,操作起来更简便.简单介绍下,LocalDate和LocalTime和LocalDateTime的使用: java.util.Da ...