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==?a:gcd(b,a%b);
 }

 LL deal(LL x,LL y)
 {
     LL res=x*(x+)/;
     vector<int>p;
     for(int i=; i*i<=y; i++)
     {
         if(y%i==)
         {
             p.push_back(i);
             while(y%i==)
             {
                 y/=i;
             }
         }
     }
     if(y>) p.push_back(y);
     LL ans=;
     for(int i=; i<(<<((int)p.size())); i++)
     {
         LL xx=;
         LL c=;
         for(int j=; j<(int)p.size(); j++)
         {
             if(i&(<<j))
             {
                 xx++;
                 c*=p[j];
             }
         }
         LL t=x/c;
         LL sum=((t+)*t/)*c;
         if(xx%)
         {
             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=; i<=m; i++)
         {
             scanf("%d",&op);
             if(op==)
             {
                 scanf("%d%d%d",&x,&y,&pp);
                 LL sum=deal((LL)y,(LL)pp)-deal((LL)(x-),(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)==)
                         {
                             sum-=it->first;
                         }
                         if(gcd(it->second,pp)==)
                         {
                             sum+=it->second;
                         }
                     }
                     it++;
                 }
                 printf("%I64d\n",sum);
             }
             else
             {
                 scanf("%d%d",&x,&y);
                 q[x]=y;
             }
         }
     }
     return ;
 }