hdu 4407 Sum

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

题意：给定初始n个数1..n，两个操作，①1 x y p  询问第x个数到第y个数中与p互质的数的和; ②：2 x y  把第x个数变成y；

思路：

把p分解质因子，然后找出(1,pos)内与p不互质的，然后用的减去就是互质的和，第二个操作用到map映射，记录在那个位置改变之后的数。

#include <cstdio>
#include <cstring>
#include <map>
#include <vector>
#include <algorithm>
#define LL __int64
#define maxn 400001
using namespace std;

int t;
int n,m;

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


LL deal(LL x,LL y)
{
    LL res=x*(x+1)/2;
    vector<int>p;
    for(int i=2; i*i<=y; i++)
    {
        if(y%i==0)
        {
            p.push_back(i);
            while(y%i==0)
            {
                y/=i;
            }
        }
    }
    if(y>1) p.push_back(y);
    LL ans=0;
    for(int i=1; i<(1<<((int)p.size())); i++)
    {
        LL xx=0;
        LL c=1;
        for(int j=0; j<(int)p.size(); j++)
        {
            if(i&(1<<j))
            {
                xx++;
                c*=p[j];
            }
        }
        LL t=x/c;
        LL sum=((t+1)*t/2)*c;
        if(xx%2)
        {
            ans+=sum;
        }
        else
            ans-=sum;
    }
    return res-ans;
}

int main()
{
    map<int,int>q;
    scanf("%d",&t);
    while(t--)
    {
        q.clear();
        scanf("%d%d",&n,&m);
        int x,y,op,pp;
        for(int i=1; i<=m; i++)
        {
            scanf("%d",&op);
            if(op==1)
            {
                scanf("%d%d%d",&x,&y,&pp);
                LL sum=deal((LL)y,(LL)pp)-deal((LL)(x-1),(LL)pp);
                map<int,int>::iterator it=q.begin();
                while(it!=q.end())
                {
                    if(it->first>=x&&it->first<=y)
                    {
                        if(gcd(it->first,pp)==1)
                        {
                            sum-=it->first;
                        }
                        if(gcd(it->second,pp)==1)
                        {
                            sum+=it->second;
                        }
                    }
                    it++;
                }
                printf("%I64d
",sum);
            }
            else
            {
                scanf("%d%d",&x,&y);
                q[x]=y;
            }
        }
    }
    return 0;
}