题目描述

输入

输出

有M行,每个询问一行,输出结果mod 1,000,000,007的值。

样例输入

10 3

3 5 1 2 3 1 3 5 2 1

7 9

3 9

2 3

样例输出

10

19

6

数据范围

对于30%的数据,N,M<=1000

对于50%的数据,N,M<=30000

对于100%的数据,N,M<=100000

解法

离线不修改区间询问,考虑莫队算法。

利用线性筛法预处理出所有要用的逆元后。

显然每次容易O(1)处理。

总的时间复杂度为O(n1.5)。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long
#define sqr(x) ((x)*(x))
#define ln(x,y) ll(log(x)/log(y))
using namespace std;
const char* fin="ex3568.in";
const char* fout="ex3568.out";
const ll inf=0x7fffffff;
const ll maxn=100007,mo=1000000007;
ll n,m,i,j,k,ks,l,r,tmp;
ll a[maxn],ans[maxn],tong[maxn],tx[maxn];
ll ny[maxn];
struct qu{
ll l,r,lk,id;
}b[maxn];
bool cmp(qu a,qu b){
return a.lk<b.lk || a.lk==b.lk && a.r<b.r;
}
ll qpower(ll a,ll b){
ll c=1;
while (b){
if (b&1) c=c*a%mo;
a=a*a%mo;
b>>=1;
}
return c;
}
void add(ll v){
tmp=(tmp+mo-tx[a[v]])%mo;
if (tong[a[v]]==0) tx[a[v]]=1;
tong[a[v]]++;
tx[a[v]]=(tx[a[v]]*a[v])%mo;
tmp=(tmp+tx[a[v]])%mo;
}
void del(ll v){
tmp=(tmp+mo-tx[a[v]])%mo;
tong[a[v]]--;
tx[a[v]]=(tx[a[v]]*ny[a[v]])%mo;
if (tong[a[v]]==0) tx[a[v]]=0;
tmp=(tmp+tx[a[v]])%mo;
}
int main(){
scanf("%d%d",&n,&m);
ks=(ll)(sqrt(n));
for (i=1;i<=n;i++) scanf("%d",&a[i]);
for (i=1;i<=100000;i++) ny[i]=qpower(i,mo-2);
for (i=1;i<=m;i++){
scanf("%d%d",&b[i].l,&b[i].r);
b[i].lk=(b[i].l-1)/ks+1;
b[i].id=i;
}
sort(b+1,b+m+1,cmp);
l=1;
r=0;
tmp=0;
for (i=1;i<=m;i++){
while (r<b[i].r) add(++r);
while (l>b[i].l) add(--l);
while (l<b[i].l) del(l++);
while (r>b[i].r) del(r--);
ans[b[i].id]=tmp;
}
for (i=1;i<=m;i++) printf("%lld\n",ans[i]);
return 0;
}

启发

涉及到预处理素数取模时的逆元,可由这道题对线筛的分析得到预处理逆元的可能性。


莫队的单次扩展一定要保证O(1)的时间。

【时光回溯】【JZOJ3568】【GDKOI2014】小纪的作业题的更多相关文章

  1. 【时光回溯】【JZOJ3567】【GDKOI2014】石油储备计划

    题目描述 输入 输出 对于每组数据,输出一个整数,表示达到"平衡"状态所需的最小代价. 样例输入 2 3 6 1 5 1 2 1 2 3 2 5 4 5 4 3 2 1 3 1 1 ...

  2. 【时光回溯】【JZOJ3566】【GDKOI2014】阶乘

    题目描述 输入 第一行有一个正整数T,表示测试数据的组数. 接下来的T行,每行输入两个十进制整数n和base. 输出 对于每组数据,输出一个十进制整数,表示在base进制下,n!结尾的零的个数. 样例 ...

  3. 【时光回溯】【JZOJ3571】【GDKOI2014】内存分配

    题目描述 输入 输出 输出m行,每行一个整数,代表输入中每次程序变化后系统所需要的空闲内存单位数. 样例输入 2 3 1 4 1 4 2 2 1 2 1 1 1 1 1 样例输出 2 3 1 数据范围 ...

  4. zzulioj - 2623: 小H的作业题

    题目链接:http://acm.zzuli.edu.cn/problem.php?id=2623 题目描述        期末考试即将来临,同学们都积极的在图书馆复习.今天小H也来到图书馆复习,小H掏 ...

  5. java lambda小纪

    一个通俗的说法是 :C#委托和Java中实现了 函数式编程的方法,它是函数式编程中的概念,为了更快的处理集合,在Java,c#等静态类型语言中想要引用一个函数的一种方式,(实现了通过封装匿名方法来达到 ...

  6. TCP/IP小纪

    链 路 层 主 要 有 三 个 目 的 :( 1 )为 I P 模 块 发 送 和 接收 I P 数 据 报 ; ( 2 )为 A R P 模块发送 A R P 请 求 和 接 收 A R P 应 答 ...

  7. Python学习小纪

    1.打包发布*.py文件---"文件路径下打开命令行 d:\python\python.exe setup.py sdist" eg:打包发布f:\C\python\print_l ...

  8. 2017年前小纪(有关http的一些缓存理论知识)

    position的top和bottom的区别:前者基准点定在top,后者基准点定在bottom. for-in 遍历属性的顺序不确定 手机端,line-height对光标大小非常有影响 有些css3属 ...

  9. vivado/FPGA 使用小纪

    1.使用FPGA做为外部控制器的总线译码时,将总线时钟接在全局时钟脚上(MRCC),就算接在了局部时钟(SRCC)上,也要通过BUFG转为全局时钟走线,否则会因为local clk到各部分的时延较大引 ...

随机推荐

  1. keras multi-label classification 多标签分类

    问题:一个数据又多个标签,一个样本数据多个类别中的某几类:比如一个病人的数据有多个疾病,一个文本有多种题材,所以标签就是: [1,0,0,0,1,0,1] 这种高维稀疏类型,如何计算分类准确率? 分类 ...

  2. c++使用优先队列时自定义优先出队顺序(和sort)

    优先队列也是一种先进先出的数据结构,元素从队尾入队,从队头出队,但是优先队列相较一般队列多了一个判断优先级的功能,在当前队列中,优先级最高的元素将被第一个删除. 先看一下优先队列的定义 templat ...

  3. c++新特性实验(5)声明与定义:属性列表(C++11 起)

    1.初识属性 1.1 实验A: noreturn 属性 [[ noreturn ]] static void thread1(void *data){ cout << "nore ...

  4. 【python之路39】Python 正则表达式

    Python 正则表达式 正则表达式是一个特殊的字符序列,它能帮助你方便的检查一个字符串是否与某种模式匹配. Python 自1.5版本起增加了re 模块,它提供 Perl 风格的正则表达式模式. r ...

  5. LUOGU P1937 [USACO10MAR]仓配置Barn Allocation

    传送门 解题思路 扫了一眼觉得是贪心+线段树,结果贪心的时候刚开始按区间长度排的序..这还有82分,后来叉了自己,换成按右端点排序过了. 代码 #include<iostream> #in ...

  6. webstorm格式化代码快捷键

    ctrl+alt+L 把网易云音乐的快捷键关了就可以了

  7. Linux安装Desktop 和 vncserver

    sudo su - #使用 root 账户 yum grouplist #查看所有可用的group yum groupinstall GNOME Desktop #安装 GNOME 桌面 yum -y ...

  8. 【maven】maven pom文件详解

    maven pom文件详解 最近配置maven中的pom文件,maven中有些属性不太清楚,在这里记录一下 <project xmlns="http://maven.apache.or ...

  9. Color the ball HDU - 1556 (线段树)

    思路:线段树,区间更新 #include<iostream> #include<vector> #include<string> #include<cmath ...

  10. html5的离线存储问题集合

    HTML5的离线存储使用一个manifest文件来标明哪些文件是需要被存储的,使用如 来引入一个manifest文件,这个文件的路径可以是相对的,也可以是绝对的,如果你的web应用很多,而且希望能集中 ...