HDU 4407

http://acm.hdu.edu.cn/showproblem.php?pid=4407

把修改和询问分成两部分解决

询问求区间内与p不互素的和，和求个数一样，用容斥原理解决，只不过做容斥的时候把每一段的个数改成每一段的和，这个求和的方式我一下写搓了，导致这道题坑了很久

修改用map记录，每次扫一遍即可，原数和c互素就减掉，修改完的数和c互素就加上去，logn*m^2的复杂度

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <map>
#include <cmath>

using namespace std;

typedef __int64 ll;

ll gcd(ll a,ll b){
    return b==?a:gcd(b,a%b);
}

ll sum(ll l,ll r,ll k){
    ll num=r/k-(l-)/k;
    ll front;
    if(l%k==)front=l;
    else front=(k-l%k)+l;
    ll rear=r-r%k;
    return (front+rear)*num/;
}

ll gao(ll l,ll r,ll n){
    vector <ll> v;
    for(ll i=;i*i<=n;i++){
        if(n%i==){
            v.push_back(i);
            while(n%i==)n/=i;
        }
    }
    if(n>)v.push_back(n);
    int m=v.size();
    ll res=;
    for(int i=;i<(<<m);i++){
        int cnt=;
        ll val=;
        for(int j=;j<m;j++){
            if(i&(<<j)){
                cnt++;
                val*=v[j];
            }
        }
        if(cnt&)res+=sum(l,r,val);
        else res-=sum(l,r,val);
    }
    return (l+r)*(r-l+)/-res;
}

map <ll,ll> mp;
map <ll,ll>::iterator it;

ll gan(ll x,ll y,ll p){
    ll res=gao(x,y,p);
    for(it=mp.begin();it!=mp.end();it++){
        if(it->first>y || it->first<x)continue;
        if(gcd(it->first,p)==)res-=it->first;
        if(gcd(it->second,p)==)res+=it->second;
    }
    return res;
}

int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        int n,m;
        scanf("%d%d",&n,&m);
        mp.clear();
        while(m--){
            int op;
            scanf("%d",&op);
            ll x,y,p;
            if(op==){
                scanf("%I64d%I64d%I64d",&x,&y,&p);
                if(x>y)swap(x,y);
                printf("%I64d\n",gan(x,y,p));
            }
            else{
                scanf("%I64d%I64d",&x,&p);
                mp[x]=p;
            }
        }
    }
    return ;
}