[题解]bzoj3813 奇数国

Description

在一片美丽的大陆上有100000个国家，记为1到100000。这里经济发达，有数不尽的账房，并且每个国家有一个银行。某大公司的领袖在这100000个银行开户时都存了3大洋，他惜财如命，因此会不时地派小弟GFS清点一些银行的存款或者让GFS改变某个银行的存款。该村子在财产上的求和运算等同于我们的乘法运算，也就是说领袖开户时的存款总和为300000。这里发行的软妹面额是最小的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
0 1 3
1 1 5
0 1 3
1 1 7
0 1 3
0 2 3

Sample Output

18
24
36
6

初始化每个国家存款都为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

Solution

用long long压位记录某个数是否含有某个质数（最多有60个质数），用欧拉函数公式+线段树直接做。

代码：

#include<cstdio>
#include<algorithm>
using namespace std;

template<typename T>inline void read(T &x){
    T f=1;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
    for(x=0;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';
    x*=f;
}

typedef long long LL;
const int maxn=100010,mod=19961993;

pair<LL,LL> Merge(pair<LL,LL> a,pair<LL,LL> b){
    return make_pair(a.first*b.first%mod,a.second|b.second);
}
struct Segment_Tree{
    #define lc x<<1
    #define rc x<<1|1
    int L[maxn<<2],R[maxn<<2];
    LL sum[maxn<<2],tag[maxn<<2];
    void Build(int x,int l,int r){
        if((L[x]=l)==(R[x]=r)){
            sum[x]=3;tag[x]=4;
            return;
        }
        int mid=(l+r)>>1;
        Build(lc,l,mid);Build(rc,mid+1,r);
        sum[x]=sum[lc]*sum[rc]%mod;
        tag[x]=tag[lc]|tag[rc];
    }
    void Modify(int x,int pos,LL val,LL t){
        if(L[x]==R[x]){
            sum[x]=val;tag[x]=t;
            return;
        }
        int mid=(L[x]+R[x])>>1;
        if(pos<=mid)Modify(lc,pos,val,t);
        else Modify(rc,pos,val,t);
        sum[x]=sum[lc]*sum[rc]%mod;
        tag[x]=tag[lc]|tag[rc];
    }
    pair<LL,LL> Query(int x,int l,int r){
        if(L[x]>r||R[x]<l)return make_pair(1,0);
        if(L[x]>=l&&R[x]<=r)return make_pair(sum[x],tag[x]);
        int mid=(L[x]+R[x])>>1;
        return Merge(Query(lc,l,r),Query(rc,l,r));
    }
}tree;
int n,prime[maxn];
LL inv[maxn];
bool ok[maxn];

LL pow(LL x,LL y){
    LL res=1;
    for(;y;y>>=1,x=x*x%mod)if(y&1)res=res*x%mod;
    return res;
}
void Make(int mx){
    for(int i=2;i<=mx;i++){
        if(!ok[i])prime[++prime[0]]=i;
        for(int j=1;j<=prime[0]&&i*prime[j]<=mx;j++){
            ok[i*prime[j]]=true;
            if(!(i%prime[j]))break;
        }
    }
    for(int i=1;i<=prime[0];i++)inv[i]=pow(prime[i],mod-2);
}
LL Div(LL x){
    LL res=0;
    for(int i=1;i<=prime[0];i++)
        if(!(x%prime[i]))res|=(1LL<<i);
    return res;
}

int main(){
    read(n);Make(281);
    tree.Build(1,1,100000);
    while(n--){
        int opt,l,r;
        read(opt);read(l);read(r);
        if(!opt){
            pair<LL,LL>temp=tree.Query(1,l,r);
            for(int i=1;i<=prime[0];i++)
                if(temp.second&(1LL<<i))
                    temp.first=temp.first*(prime[i]-1)%mod*inv[i]%mod;
            printf("%lld\n",temp.first);
        }
        else tree.Modify(1,l,r,Div(r));
    }
    return 0;
}