bzoj 3813: 奇数国 （线段树+质因数分解+线性筛+线性逆元）

3813: 奇数国

Time Limit: 10 Sec   Memory Limit: 256 MB
Submit: 502   Solved: 295
[ Submit][ Status][ Discuss]

Description

在一片美丽的大陆上有100000个国家，记为1到100000。这里经济发达，有数不尽的账房，并且每个国家有一个银行。某大公司的领袖在这100000个银行开户时都存了3大洋，他惜财如命，因此会不时地派小弟GFS清点一些银行的存款或者让GFS改变某个银行的存款。该村子在财产上的求和运算等同于我们的乘法运算，也就是说领袖开户时的存款总和为3100000。这里发行的软妹面额是最小的60个素数（p1=2,p2=3,…,p60=281），任何人的财产都只能由这60个基本面额表示，即设某个人的财产为fortune（正整数），则fortune=p1^k1*p2^k2*......p60^K60。

领袖习惯将一段编号连续的银行里的存款拿到一个账房去清点，为了避免GFS串通账房叛变，所以他不会每次都选择同一个账房。GFS跟随领袖多年已经摸清了门路,知道领袖选择账房的方式。如果领袖选择清点编号在[a,b]内的银行财产，他会先对[a,b]的财产求和（计为product），然后在编号属于[1,product]的账房中选择一个去清点存款，检验自己计算是否正确同时也检验账房与GFS是否有勾结。GFS发现如果某个账房的编号number与product相冲，领袖绝对不会选择这个账房。怎样才算与product不相冲呢？若存在整数x,y使得number*x+product*y=1，那么我们称number与product不相冲，即该账房有可能被领袖相中。当领袖又赚大钱了的时候，他会在某个银行改变存款，这样一来相同区间的银行在不同的时候算出来的product可能是不一样的，而且领袖不会在某个银行的存款总数超过1000000。

现在GFS预先知道了领袖的清点存款与变动存款的计划，想请你告诉他，每次清点存款时领袖有多少个账房可以供他选择，当然这个值可能非常大，GFS只想知道对19961993取模后的答案。

Input

第一行一个整数x表示领袖清点和变动存款的总次数。

接下来x行，每行3个整数ai,bi,ci。ai为0时表示该条记录是清点计划，领袖会清点bi到ci的银行存款，你需要对该条记录计算出GFS想要的答案。ai为1时表示该条记录是存款变动，你要把银行bi的存款改为ci，不需要对该记录进行计算。

Output

输出若干行，每行一个数，表示那些年的答案。

Sample Input

6
013
115
013
117
013
023

Sample Output

18
24
36
6

explanation

初始化每个国家存款都为3;

1到3的product为27，[1,27]与27不相冲的有18个数;
1的存款变为5;
1到3的product为45，[1,45]与45不相冲的有24个数;
1的存款变为7;
1到3的product为63，[1,63]与63不相冲的有36个数;
2到3的product为9，[1,9]与9不相冲的有6个数。

HINT

x≤100000，当ai=0时0≤ci−bi≤100000

Source

2015年国家集训队测试

[ Submit][ Status][ Discuss] 

题解：线段树+质因数分解+线性筛+线性逆元

对于所有位置的数质因数分解，然后维护区间中质因数的出现次数。

number*x+product*y=1 => number*x=1 (mod product) 

同余方程有解，当且仅当number,product互质。

所有我们要求的就是[1,product]中与product互质的数的个数，即phi(product)

然后利用公式求解即可。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 100003
#define p 19961993
#define LL long long
using namespace std;
LL mp[303],pd[303],prime[303],inv[303];
struct data {
	int cnt[63];
	data() {
		memset(cnt,0,sizeof(cnt));
	}
}tr[N*4];
int n;
void init()
{
	for (int i=2;i<=300;i++) {
		if (!pd[i]) prime[++prime[0]]=i,mp[i]=prime[0];
		for (int j=1;j<=prime[0];j++) {
			if (prime[j]*i>300) break;
			pd[prime[j]*i]=1;
			if (i%prime[j]==0) break;
		}
	}
	prime[0]=60;
	inv[0]=inv[1]=1;
	for (int i=2;i<=281;i++) 
	 inv[i]=(p-p/i)*inv[p%i]%p;
	//for (int i=2;i<=261;i++) cout<<inv[i]<<" ";
	//cout<<endl;
}
void update(int now)
{
	for (int i=1;i<=60;i++) tr[now].cnt[i]=tr[now<<1].cnt[i]+tr[now<<1|1].cnt[i];
}
void build(int now,int l,int r)
{
	if (l==r) {
		memset(tr[now].cnt,0,sizeof(tr[now].cnt));
		tr[now].cnt[2]=1;
		return;
	}
	int mid=(l+r)/2;
	build(now<<1,l,mid);
	build(now<<1|1,mid+1,r);
	update(now);
}
void add(data &a,data b)
{
	for (int i=1;i<=60;i++)
	 a.cnt[i]+=b.cnt[i];
}
data query(int now,int l,int r,int ll,int rr)
{
	if (ll<=l&&r<=rr) return tr[now];
	int mid=(l+r)/2; data a; 
	if (ll<=mid) add(a,query(now<<1,l,mid,ll,rr));
	if (rr>mid) add(a,query(now<<1|1,mid+1,r,ll,rr));
	return a;
}
void pointchange(int now,int l,int r,int x,int val)
{
	if (l==r) {
		memset(tr[now].cnt,0,sizeof(tr[now].cnt));
		int k=val;
		for (int i=1;prime[i]*prime[i]<=k;i++) 
		 if (k%prime[i]==0) {
		 	while (k%prime[i]==0) k/=prime[i],tr[now].cnt[i]++;
		 }
		if (k>1) tr[now].cnt[mp[k]]++;
		return;
	}
	int mid=(l+r)/2;
	if (x<=mid) pointchange(now<<1,l,mid,x,val);
	else pointchange(now<<1|1,mid+1,r,x,val);
	update(now);
}
LL quickpow(LL num,int x)
{
	LL base=num%p; LL ans=1;
	while (x) {
		if (x&1) ans=ans*base%p;
		x>>=1;
		base=base*base%p;
	}
	return ans;
}
int main()
{
	freopen("a.in","r",stdin);
	freopen("my.out","w",stdout);
	init(); int T;
	scanf("%d",&T); n=100000;
	build(1,1,n);
	for (int i=1;i<=T;i++) {
		int opt,x,y;
		scanf("%d%d%d",&opt,&x,&y);
		if (opt==0) {
		  data a=query(1,1,n,x,y);
		  LL ans=1;
		  for (int j=1;j<=60;j++)
		   if (a.cnt[j]) 
		    ans=ans*quickpow(prime[j],a.cnt[j])%p*(LL)(prime[j]-1)%p*inv[prime[j]]%p;
		    //cout<<quickpow(prime[j],a.cnt[j])%p<<" "<<inv[prime[j]]<<endl;
		  printf("%I64d\n",ans%p);
	    }
	    else pointchange(1,1,n,x,y);
	}
}