传送门

题意

给出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. mui + vue 模板

    示例代码: (function(mui, doc) { // 定义全局变量(计时器) var timer; // mui初始化 mui.init(); // 创建vue的实例 var app = ne ...

  2. HDU 1003 Max Sum (动规)

    Max Sum Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Su ...

  3. C#.net word excel powerpoint (ppt) 转换成 pdf 文件

    using System;using System.Collections.Generic;using System.ComponentModel;using System.Data;using Sy ...

  4. 6.非关系型数据库(Nosql)之mongodb:集群(主从复制)

     1.主从复制是MongoDB最经常使用的复制方式.这样的方式很灵活,可用于备份.故障恢复.读扩展等 2最主要的设置方式就是建立一个主节点和一个或多个从节点,每一个从节点要知道主节点的地址. 执行 ...

  5. INAPP登陆调用的FB接口

    public function login_get (){ $this->load->helper ( 'auth' ); $redirectUrl = $this->input-& ...

  6. long long , __int64 范围

    VC的64位整数 分别叫做__int64与unsigned __int64,其范 围分别是[-2^63, 2^63)与[0,2^64),即-9223372036854775808~9223372036 ...

  7. 转载:用python爬虫抓站的一些技巧总结

    原文链接:http://www.pythonclub.org/python-network-application/observer-spider 原文的名称虽然用了<用python爬虫抓站的一 ...

  8. Jmeter使用Http代理服务器报DNSName components must begin with a letter的错

    最近了解到JMeter可以实现app的性能测试,需要借助JMeter的Http代理服务器来录制脚本. 于是,就按着网上的教程来进行操作,然而出师不利啊,刚启动就报错

  9. STM32 DMA中断只进入一次的解决办法

    问题解决参见:http://bbs.ednchina.com/BLOG_ARTICLE_3014819.HTM 经过我验证,这个说的是对的.

  10. rails find find_by where

    find根据id进行查询,像Product.find(3),查询语句是Product Load (0.1ms) SELECT "products".* FROM "pro ...