题目

数据范围

分析

时限5000ms。

我们注意到\(a_{i}初始值以及x小于等于600且非零\)

也就是说,\(a_{i}\)的质因数一定小于600,而600以内的质因数只有109个。

那么考虑常用于区间修改的线段树。

用线段树来维护某个位置的某个质因数的总乘积,以及某个质因数出现的位置的个数。

时间复杂度\(O(QlogN·109)\)

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const long long maxlongint=2147483647;
const long long mo=100000007;
const long long N=50005;
using namespace std;
struct ddx
{
long long a[120],v[120],lazy[120],sum[120];
}tree[50005];
long long ss[80000],n,m,ny[1000005],belong[1000005];
long long ans,re[10005];
bool b[1000005];
long long mi(long long x,long long y)
{
long long sum=1;
while(y)
{
if(y&1) sum=x*sum%mo;
x=x*x%mo;
y/=2;
}
return sum;
}
long long put(long long v,long long l,long long r,long long x,long long y,long long z)
{
if(l==r)
{
tree[v].a[y]=z;
tree[v].v[y]=tree[v].v[y]*mi(ss[y],z)%mo;
tree[v].sum[y]=1;
return 0;
}
long long mid=(l+r)/2;
if(x<=mid)
put(v*2,l,mid,x,y,z);
else put(v*2+1,mid+1,r,x,y,z);
for(long long i=1;i<=109;i++)
{
tree[v].a[i]=(tree[v*2].a[i]+tree[v*2+1].a[i])%mo;
tree[v].v[i]=tree[v*2].v[i]*tree[v*2+1].v[i]%mo;
tree[v].sum[i]=tree[v*2].sum[i]+tree[v*2+1].sum[i];
}
}
long long down(long long v,long long mid,long long l,long long r,long long y)
{
if(!tree[v].lazy[y]) return 0;
long long z=tree[v].lazy[y];
tree[v*2].a[y]+=z*(mid-l+1);
tree[v*2+1].a[y]+=z*(r-mid);
tree[v*2].v[y]=tree[v*2].v[y]*mi(ss[y],z*(mid-l+1))%mo;
tree[v*2+1].v[y]=tree[v*2+1].v[y]*mi(ss[y],z*(r-mid))%mo;
tree[v*2].lazy[y]+=z;
tree[v*2+1].lazy[y]+=z;
tree[v*2].sum[y]=(mid-l+1);
tree[v*2+1].sum[y]=(r-mid);
tree[v].lazy[y]=0;
}
long long change(long long v,long long l,long long r,long long x,long long x1,long long y,long long z)
{
if(l==x && r==x1)
{
tree[v].a[y]=(tree[v].a[y]+z*(r-l+1))%mo;
tree[v].v[y]=tree[v].v[y]*mi(ss[y],z*(r-l+1))%mo;
tree[v].sum[y]=(r-l+1);
tree[v].lazy[y]+=z;
return 0;
}
long long mid=(l+r)/2;
for(long long i=1;i<=109;i++) down(v,mid,l,r,i);
if(x1<=mid)
change(v*2,l,mid,x,x1,y,z);
else
if(x>mid)
change(v*2+1,mid+1,r,x,x1,y,z);
else
change(v*2,l,mid,x,mid,y,z),change(v*2+1,mid+1,r,mid+1,x1,y,z);
for(long long i=1;i<=109;i++)
{
tree[v].a[i]=(tree[v*2].a[i]+tree[v*2+1].a[i])%mo;
tree[v].v[i]=tree[v*2].v[i]*tree[v*2+1].v[i]%mo;
tree[v].sum[i]=tree[v*2].sum[i]+tree[v*2+1].sum[i];
}
}
long long find(long long v,long long l,long long r,long long x,long long x1)
{
if(l==x && r==x1)
{
for(long long i=1;i<=109;i++)
{
if(tree[v].sum[i])
ans=ans*tree[v].v[i]%mo*mi(ny[ss[i]],tree[v].sum[i])%mo*mi(ss[i]-1,tree[v].sum[i])%mo;
}
return 0;
}
long long mid=(l+r)/2;
for(long long i=1;i<=109;i++) down(v,mid,l,r,i);
if(x1<=mid)
find(v*2,l,mid,x,x1);
else
if(x>mid)
find(v*2+1,mid+1,r,x,x1);
else
find(v*2,l,mid,x,mid),find(v*2+1,mid+1,r,mid+1,x1);
for(long long i=1;i<=109;i++)
{
tree[v].a[i]=(tree[v*2].a[i]+tree[v*2+1].a[i])%mo;
tree[v].v[i]=tree[v*2].v[i]*tree[v*2+1].v[i]%mo;
tree[v].sum[i]=tree[v*2].sum[i]+tree[v*2+1].sum[i];
}
}
int main()
{
memset(b,true,sizeof(b));
b[0]=0;
b[1]=0;
for(long long i=2;i<=10000;i++)
{
if(b[i])
{
ss[++ss[0]]=i;
ny[i]=mi(i,mo-2);
belong[i]=ss[0];
}
for(long long j=1;j<=ss[0];j++)
{
if(i*ss[j]<=10000)
b[i*ss[j]]=false;
else break;
if(!(i%ss[j])) break;
}
}
scanf("%lld",&n);
for(long long i=1;i<=50001;i++)
for(long long j=1;j<=110;j++)
tree[i].v[j]=1;
for(long long j=1;j<=n;j++)
{
scanf("%lld",&re[j]);
long long p=re[j];
if(b[p])
{
put(1,1,n,j,belong[p],1);
}
else
{
for(long long k=1;ss[k]<=sqrt(re[j]) && p>1;k++)
{
long long w=0;
while(!(p%ss[k]))
{
p/=ss[k];
w++;
}
if(w)
{
put(1,1,n,j,k,w);
}
}
if(p>1)
{
put(1,1,n,j,belong[p],1);
}
}
}
long long q;
scanf("%lld",&q);
for(long long i=1;i<=q;i++)
{
long long t,x,y,z;
scanf("%lld%lld%lld",&t,&x,&y);
if(t)
{
ans=1;
find(1,1,n,x,y);
printf("%lld\n",ans);
}
else
{
scanf("%lld",&z);
long long p=z;
if(b[p])
{
change(1,1,n,x,y,belong[p],1);
}
else
{
for(long long k=1;ss[k]<=sqrt(z) && p>1;k++)
{
long long w=0;
while(!(p%ss[k]))
{
p/=ss[k];
w++;
}
if(w)
{
change(1,1,n,x,y,k,w);
}
}
if(p>1)
{
change(1,1,n,x,y,belong[p],1);
}
}
}
}
}

【NOIP2016提高组A组7.16】第三条跑道的更多相关文章

  1. JZOJ 【NOIP2016提高A组集训第16场11.15】兔子

    JZOJ [NOIP2016提高A组集训第16场11.15]兔子 题目 Description 在一片草原上有N个兔子窝,每个窝里住着一只兔子,有M条路径连接这些窝.更特殊地是,至多只有一个兔子窝有3 ...

  2. JZOJ 【NOIP2016提高A组集训第16场11.15】SJR的直线

    JZOJ [NOIP2016提高A组集训第16场11.15]SJR的直线 题目 Description Input Output Sample Input 6 0 1 0 -5 3 0 -5 -2 2 ...

  3. NOIP2016提高组解题报告

    NOIP2016提高组解题报告 更正:NOIP day1 T2天天爱跑步 解题思路见代码. NOIP2016代码整合

  4. 5820. 【NOIP提高A组模拟2018.8.16】 非法输入(模拟,字符串)

    5820. [NOIP提高A组模拟2018.8.16] 非法输入 (File IO): input:aplusb.in output:aplusb.out Time Limits: 1000 ms   ...

  5. JZOJ 4732. 【NOIP2016提高A组模拟8.23】函数

    4732. [NOIP2016提高A组模拟8.23]函数 (Standard IO) Time Limits: 1500 ms  Memory Limits: 262144 KB  Detailed ...

  6. 【题解】NOIP2016提高组 复赛

    [题解]NOIP2016提高组 复赛 传送门: 玩具谜题 \(\text{[P1563]}\) 天天爱跑步 \(\text{[P1600]}\) 换教室 \(\text{[P1850]}\) 组合数问 ...

  7. 【题解】NOIP2016 提高组 简要题解

    [题解]NOIP2016 提高组 简要题解 玩具迷题(送分) 用异或实现 //@winlere #include<iostream> #include<cstdio> #inc ...

  8. iOS之UITableView组头组尾视图/标题悬停

    最近笔者在公司的iOS开发中,有一个iOS开发同事跑来问了两个问题:1.给UITableView设置了组头和组尾视图,但是一直显示不出来?2.UITableView的section的header和fo ...

  9. iOS自定义组与组之间的距离以及视图

    iOS自定义组与组之间的距离以及视图 //头视图高度 - (CGFloat)tableView:(UITableView *)tableView heightForHeaderInSection:(N ...

随机推荐

  1. CentOS7_装机软件推荐

    目录 目录 前言 CentOS7中文支持 Install Yum Fastest Mirror Plugin Install Flash Player Install NTFS-3G Install ...

  2. Java中的内存泄露 和 JVM GC(垃圾回收机制)

    一.什么是Java中的内存泄露? 在Java中,内存泄漏就是存在一些被分配的对象,这些对象有下面两个特点, 首先,这些对象是可达的,即在有向图中,存在通路可以与其相连:其次,这些对象是无用的,即程序以 ...

  3. mysql5.6修改密码并授权所有远程用户可登陆

    1.my.ini文件,删除最后一行的"skip-grant-tables 2.执行"use mysql;",使用mysql数据库; 3.执行:update mysql.u ...

  4. 加载动态链接库——dlopen dlsym dlclose

    DLOPEN DLMOPEN DLCLOSE NAME     dlclose, dlopen, dlmopen - 打开/关闭共享对象 SYNOPSIS #include <dlfcn.h&g ...

  5. Go语言入门篇-基本类型排序和 slice 排序

    参见博客:https://blog.csdn.net/u010983881/article/details/52460998 package main import ( "sort" ...

  6. 20191110 Spring Boot官方文档学习(3)

    3.使用Spring Boot 3.1.构建系统 建议选择Maven或Gradle作为构建工具 每个Spring Boot版本都提供了它所支持的依赖关系的精选列表.实际上,您不需要为构建配置中的所有这 ...

  7. dotnet sdk 的镜像tag 相关

    https://hub.docker.com/_/microsoft-dotnet-core-sdk/ 微软的dotnet sdk 的 tag 微软貌似改默认镜像 dockerhub 里面的tag与 ...

  8. C#WebApi自动生成文档

    1.效果图 2.在webApi项目,打开Nuget,搜索WebApiTestClient,安装WebApiTestClient,注意是给HelpPage的 3.打开引入WebApiTestClient ...

  9. typescript是否可以直接编译执行?

    算是个有趣的小问题,由于必须依赖node.js,typescript理论上是不能不转成js直接运行的.

  10. centos7yum安装VirtualBox

    cd 进入目录:/etc/yum.repos.d 新建一个文件virtualbox.repo, 输入如下内容: [virtualbox] name=Oracle Linux / RHEL / Cent ...