0707-线段树+打表-好数

今天来讲讲T3吧

好数

 

题目描述

我们定义一个非负整数是“好数”，当且仅当它符合以下条件之一：
    1. 这个数是 0 或 1 。
    2. 所有小于这个数且与它互质的正整数可以排成一个等差数列，例如，8 就是一个好数，因为 1，3，5，7 排成了等差数列。
给出 N 个非负整数，然后进行如下三个操作：
1. 询问区间 [L,R] 有多少个好数。
2. 将区间 [L,R] 内所有数对 S 取余（S≤1000000）。
3. 将第 C 个数更改为 X 。

提示：如果你不知道如何判断一个数是否为好数，你可以打个表找找规律。
 

输入格式

第一行包含两个正整数 N 和 M ，M 表示操作数目。
第二行包含 N 个非负整数。
接下来的 M 行每行表示 1 个操作：“1 L R”表示第 1 个操作，“2 L R S”表示第 2 个操作，“3 C X”表示第 3 个操作。
 

输出格式

对每个操作1，输出一个非负整数，表示区间内好数的个数。
 

样例数据 1

输入　 
3 6 
4 6 9 
1 1 3 
1 3 3 
2 1 1 10 
1 1 3 
3 2 4 
1 1 3
输出
2 
0 
2 
2

样例数据 2

输入　
8 5 
12 24 17 31 16 21 18 30 
1 2 5 
2 4 7 7 
3 2 13 
1 1 8 
1 3 6
输出
3 
6 
4
备注
【数据规模与约定】

在正式讲这道题之前，我们先储备一些知识：

1.每个数取余后要么不变化，要么小于它的1/2（证明什么的，自己想想）

2.根据提示可知我们要寻找好数的特征（找到后超级有用，在我原来的程序上多过了三组，其他还是TLE）

那先打个表吧

（0 1） 2 3 4 5 6 7 8 11 13 16 17 19 23 29 31 32 37 41 43 47 53 59 61 64 67 71 73 79 83 89 97 

然后按照我下面的顺序列一列（竖着看）

2      3 
4      5 
6      7 
8      11 13 
16    17 19 23 29 31 
32    37 41 43 47 53 59 61 

64    67 71 73 79 83 89 97 

我们可以惊奇的发现只要这个数是2的k次幂、质数或者是6、0、1，那么它就是一个好数，但这只是我们看出来的，还是得证明一下咯，根据好数的定义可以很容易的证明我们结论的正确性（反正我当时打了表也没看出来，可能是因为打的太多了，根本看不出来，那么提醒自己下次打表找规律的时候数据不需要太大）

                                                                                                                                                                         好啦     

这些东西知道后，我们就开始讲讲正解是怎么搞的

这道题一看就感觉像是要用线段树，而且对于操作1、3，区间询问和单点修改简直就是线段树的长处（所以我严重怀疑我TLE是因为操作2——区间修改，我每次都是把区间内的每一个点修改了，然后再合并起来，比暴力还不如，再加上判断好数的时候没有找到它的特性，所以妥妥地TLE）所以的所以，现在问题的关键，就是如何使用线段树来解决区间取模，那么由于我们的知识储备1中提到的，每个数x在经历log (x)次的取余后就变成1了（其实区间开方也是同样的道理），那么对于已经变成1的数我们是不是就可以不再对它取余啦，这样就会节省很多时间。so具体我们怎么操作呢，对于每一个区间我们维护一下它的最大值，若此次模数大于最大值，则不用对其进行修改

下面上代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<iostream>
using namespace std;
int ifgood[1050006],vst[1000009],pri[1000009];
int maxn[5000009];
void init(){
	memset(vst,0,sizeof(vst));
	int i,j,k,tot=0;
	for(i=2;i<=1000000;++i){
		if(!vst[i]) {
			pri[++tot]=i;
			ifgood[i]=1;
		}
		for(j=1;j<=tot&&i*pri[j]<=1000000;++j){
			vst[i*pri[j]]=1;
			if(i%pri[j]==0)	break;
		}
	}
	ifgood[6]=1;
	ifgood[0]=1;
	ifgood[1]=1;
}
int n,m,x;
inline int read(){
	char ch;
	bool flag=0;
	while((ch=getchar())<'0'||ch>'9'){
		if(ch=='-') flag=1;
	}
	int res=0;
	while(ch>='0'&&ch<='9'){
		res=res*10+ch-'0';
		ch=getchar();
	}
	return flag?-res:res;
}
int sum[500000],a[100009];
void build(int k,int l,int r){
	if(l==r){
		if(ifgood[a[l]]){
			sum[k]++;
		} 
		maxn[k]=a[l];
		return ;
	}
	int mid=l+r>>1;
	build(2*k,l,mid);
	build(2*k+1,mid+1,r);
	sum[k]=sum[2*k]+sum[2*k+1];
	maxn[k]=max(maxn[k<<1],maxn[k<<1|1]);
}
int query(int k,int l,int r,int x,int y){
	if(x<=l&&y>=r)	return sum[k];
	int mid=l+r>>1,res=0;
	if(x<=mid) res+=query(2*k,l,mid,x,y);
	if(y>mid) res+=query(2*k+1,mid+1,r,x,y);
	return res;
}
void modify1(int k,int l,int r,int x,int y,int s){
	if(maxn[k]<s) return ;
	if(l==r){
		a[l]=a[l]%s;
		maxn[k]=maxn[k]%s;
		sum[k]=ifgood[a[l]];
		return ;
	}
	int mid=l+r>>1;
	if(x<=mid) modify1(2*k,l,mid,x,y,s);
	if(y>mid) modify1(2*k+1,mid+1,r,x,y,s);
	sum[k]=sum[2*k]+sum[2*k+1];
	maxn[k]=max(maxn[k<<1],maxn[k<<1|1]);
}
void modify2(int k,int l,int r,int c,int x){
	if(l==r&&l==c){
		a[l]=x;
		maxn[k]=x;
		sum[k]=ifgood[x];
		return ;
	}
	int mid=l+r>>1;
	if(c<=mid) modify2(2*k,l,mid,c,x);
	if(c>mid) modify2(2*k+1,mid+1,r,c,x);
	sum[k]=sum[2*k]+sum[2*k+1];
	maxn[k]=max(maxn[k<<1],maxn[k<<1|1]);
}
int main(){
//	freopen("good.in","r",stdin);
//	freopen("good.out","w",stdout);
	int i,j;
	n=read();m=read();
	for(i=1;i<=n;++i)
	{
		a[i]=read();
	}
	memset(ifgood,0,sizeof(ifgood));
	int hh=1;
	for(i=1;i<=20;++i){
		hh=hh*2;
		ifgood[hh]=1;
	}
	init();
	build(1,1,n);
	for(i=1;i<=m;++i){
		int id,l,r,s,c,x;
		id=read();
		if(id==1){
			l=read();r=read();
			printf("%d\n",query(1,1,n,l,r));
		}
		else if(id==2){	
				l=read();r=read();s=read();
				modify1(1,1,n,l,r,s);		
			}
			else{
				c=read();x=read();
				modify2(1,1,n,c,x);				
			}
	}
//	fclose(stdin);fclose(stdout);
	return 0;
}

末尾吐槽一下自己：我单点修改后居然只改了sum的值！！！！区间修改后居然也只改了sum！！！！我居然没有找出好数的特性！！！！！只是把这三个更改之后，就瞬间从0蛋蛋变成了70分！！！！！再加上区间修改的优化，就直接AC了