题意:求这个式子 $\sum \limits_{i=1}^{n} \sum \limits_{j=1}^{m} f(i, j) mod (10^9 + 7)$ 的值

就是对每个区间[i, j]枚举区间中的每个数$a_i$到$a_j$, 判断这个$a$是否对[i, j]这个区间内所有数取模都不等于0, 若是,则这个区间满足条件

问有多少个满足条件的区间

比如案例是这样跑的

    int ans=;
for(int i=;i<=;i++)
for(int j=i;j<=;j++)
{
for(int k=i;k<=j;k++) // 注意要枚举[i, j]中的每个数
{
bool flag=;
for(int l=i;l<=j;l++)
if(k!=l && k%l==)
flag=;
if(!flag) // 对区间内所有数取模都不等于0
ans++;
}
}

跟省赛某题很像, 计算每个数a[i]对ans的贡献

比如对于案例 1 2 3 4 5

1这个数字对于答案的贡献是{1}, {1, 2}, {1, 2, 3}, {1, 2, 3, 4}, {1, 2, 3, 4, 5} 这5个区间

2 这个数字对于答案的贡献是{2}, {2, 3}, {2, 3, 4}, {2, 3, 4, 5} 这4个区间   ( {1, 2}区间不满足,因为2%1==0 )

... ...

来看4 这个数

它往左取区间 {4}, {3, 4} 当取到2的时候, 发现4%2==0了,那么就不必再往左了(对于连续的区间, 再往左则必定会经过2,那么该区间就不合法了)

同理,可以想到往右取, 当找到一个数被它取模等于0, 那么就不必再往右了

好了,现在对于一个数a[i], 它左边 到 不合法的数(被它取模等于0)为止 之间有x个数, 右边到不合法的数为止 有y个数

那么a[i]这个数对答案的贡献就是(x+1)*(y+1)

为什么呢?

因为是连续的区间, 所以这个区间的左端点可以取a[i]左边0个数、1个数、2个数... ...x个数;右边0个数、1个数、2个数... ...y个数

左边有(x+1)种取法,右边有(y+1)种, 相乘就是总取法数

那么我们只要找到离a[i]最近的一左一右两个不合法数的位置$l$和$r$, 那么$(i-l)*(r-i)$ 就是a[i]的贡献($i-l-1$就是上面所讲的x)

之后只要遍历每个a, 将每个a的贡献累加起来即是最后答案。

那么现在问题就转化成了如何求一个 离它最近的 能被它整除的数 的位置

我们是这样做的:

开个数组将每个a[i]的倍数都记为i  (比如 a[0]=2, 那么就将2、4、6... ...10000 都记下0号位置;a[x]=y, 就将y、2y、3y... ... 都记下x位置)

就跟筛因子一样   (因为1比较特殊,会退化到$n^2$, 因此特殊处理)   复杂度为O(NlogN)  (N为10000, 因为数最大为10000)

        for(int i=;i<n;i++)
{
if(a[i]==)
one.push_back(i);
else
for(int j=a[i];j<=;j+=a[i])
b[j].push_back(i);
}

因为是按顺序遍历了a数组, 所以记下的位置(比如2 2 3 6 12  对于12记下的是0 1 2 3 4 )一定是递增的

那么就可以二分来寻找离$i$最近的位置

p.s. lower_bound 找的是大于等于x的数位置

upper_bound找的是大于x的数的位置

 const LL mod=1e9+;
int a[];
vector<int> b[];
vector<int> one;
int read()
{
char ch=' ';
int ans=;
while(ch<'' || ch>'')
ch=getchar();
while(ch<='' && ch>='')
{
ans=ans*+ch-'';
ch=getchar();
}
return ans;
} int main()
{
// freopen("1001.in", "r", stdin);
// freopen("out.txt", "w", stdout);
int n;
while(~scanf("%d", &n))
{
for(int i=;i<n;i++)a[i]=read();
// scanf("%d", &a[i]);
for(int i=;i<=;i++)
{
b[i].clear();
b[i].push_back(-);
}
one.clear();
one.push_back(-);
for(int i=;i<n;i++)
{
if(a[i]==)
one.push_back(i);
else
for(int j=a[i];j<=;j+=a[i])
b[j].push_back(i);
}
for(int i=;i<=;i++)
b[i].push_back(n);
one.push_back(n);
LL ans=;
for(int i=;i<n;i++)
{
int p1=lower_bound(b[a[i]].begin(), b[a[i]].end(), i)-b[a[i]].begin()-;
int p2=lower_bound(one.begin(), one.end(), i)-one.begin()-;
int l=max(b[a[i]][p1], one[p2]);
p1=upper_bound(b[a[i]].begin(), b[a[i]].end(), i)-b[a[i]].begin();
p2=upper_bound(one.begin(), one.end(), i)-one.begin();
int r=min(b[a[i]][p1], one[p2]);
// printf("%d %d\n", l, r);
ans=(ans+((i-l)*(r-i))%mod)%mod;
}
printf("%I64d\n", ans%mod);
}
return ;
}

HDOJ 5288

[YY题]HDOJ5288 OO’s Sequence的更多相关文章

  1. HDU 5288 OO’s Sequence 水题

    OO's Sequence 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5288 Description OO has got a array A ...

  2. HDU 5288——OO’s Sequence——————【技巧题】

    OO’s Sequence Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)T ...

  3. Hdu 5288 OO’s Sequence 2015多小联赛A题

    OO's Sequence Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others) ...

  4. 思维题 HDOJ 5288 OO’s Sequence

    题目传送门 /* 定义两个数组,l[i]和r[i]表示第i个数左侧右侧接近它且值是a[i]因子的位置, 第i个数被选择后贡献的值是(r[i]-i)*(i-l[i]),每个数都枚举它的因子,更新l[i] ...

  5. HDU 5288 OO‘s sequence (技巧)

    题目链接:http://acm.hdu.edu.cn/showproblem.php? pid=5288 题面: OO's Sequence Time Limit: 4000/2000 MS (Jav ...

  6. HDOJ 5288 OO’s Sequence 水

    预处理出每一个数字的左右两边能够整除它的近期的数的位置 OO's Sequence Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 13 ...

  7. HDU 5288 OO’s Sequence [数学]

     HDU 5288 OO’s Sequence http://acm.hdu.edu.cn/showproblem.php?pid=5288 OO has got a array A of size ...

  8. 2015 Multi-University Training Contest 1 - 1001 OO’s Sequence

    OO’s Sequence Problem's Link: http://acm.hdu.edu.cn/showproblem.php?pid=5288 Mean: 给定一个数列,让你求所有区间上满足 ...

  9. hdu 5288 OO’s Sequence(2015 Multi-University Training Contest 1)

    OO's Sequence                                                          Time Limit: 4000/2000 MS (Jav ...

随机推荐

  1. 什么是WEB服务器?

    简单的理解:  一个公用的服务器,你可以放置网站文件,让全世界浏览,可以放文件,全世界下载.web服务器一般指网站服务器.  web服务器是基于网站架设的服务器,我们平时可以浏览的网页都是在别人的服务 ...

  2. CentOS学习笔记--程序管理

    程序管理 一个程序被加载到内存当中运行,那么在内存内的那个数据就被称为程序(process).程序是操作系统上非常重要的概念, 所有系统上面跑的数据都会以程序的型态存在.那么系统的程序有哪些状态?不同 ...

  3. linux 关机方式

    linux 关机命令: 1-  init 0 关机. 具体详情接受可以 init --help 查询 如下: init [OPTIONS...] {COMMAND} Send control comm ...

  4. js 如何获取文本框中光标索引位置

    function getTxt1CursorPosition(){ var oTxt1 = document.getElementById("txt1"); var cursurP ...

  5. 使用 rem 实现 适配各种屏幕布局

    年数已久!!!技术更新太快,谨慎观看,不过作为了解一下思路还是可以的 在此呢,我只大略的谈一下在研究了px,em,rem,和手淘的做法之后,我所采用的做法及硬性规则: 我就不再过多的将上面三种单位的区 ...

  6. linux服务器报No space left on device错误的解决过程记录

    起因 今天在本地提交了点代码,但到服务器上git pull的时候提示No space left on device,第一反应是猜想可能硬盘满了(很有可能是log导致的),不过想想又觉得不太可能,这台服 ...

  7. liunx下vi命令详解

    vi编辑器是所有Unix及Linux系统下标准的编辑器,它的强大不逊色于任何最新的文本编辑器,这里只是简单地介绍一下它的用法和一小部分指令.由于对Unix及Linux系统的任何版本,vi编辑器是完全相 ...

  8. 修改UI中的值,无反应解决办法

    var targetObj = $("<input name='mydate' class='easyui-datebox'>").appendTo("#id ...

  9. Boa练习程序2

    做一个地址簿的gui. #Boa:Frame:AddressEntry import wx def create(parent): return AddressEntry(parent) [wxID_ ...

  10. WPF学习笔记 控件篇 属性整理【1】FrameworkElement

    最近在做WPF方面的内容,由于好多属性不太了解,经常想当然的设置,经常出现自己未意料的问题,所以感觉得梳理下. ps:先补下常用控件的类结构,免得乱了 .NET Framework 4.5 Using ...