hdoj4407题解
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了hdoj4407题解相关的知识,希望对你有一定的参考价值。
题目描述:1,2,3.....n排一行(1<=n<=400000),接下来有m(1<=m<=1000)操作,操作分为两种,一种是查询,一种是改值。
查询:输入x,y,p,求第x个数到第y个数之间与p互质的数的和
改值:输入x,c,将第x个数的值改为c
这题用到容斥原理,先认为第x到第y个数是连续的,对于连续的数求与p互质的数的和可以用第x个数到第y个数的和减去与p有公约数的和。
对于改值,用个数组记录下来要修改的数的下标和修改之后的值
连续的数列求完之后再遍历这个数组,若修改的数的下标落在x到y这个区间之内,则进行替换。
400000以内的数的质约数的最多只有7个,所以求连续的数列与p互质的和的时间复杂度是2*7*2^7,不超过1000次,遍历数组也不会超过1000次,所以对于任意一步操作计算的次数不超过2000次,再乘以1000步操作,所以计算量不会超过2000000次。
#include <iostream> #include <cstdio> #include <cmath> using namespace std; int n; int m; int yueshu[20]; int nyueshu; struct Node { int index; int value; }; Node node[1010]; int nnode; void Initialize() { nnode=0; } int IsCoprime(int a) { int i; for(i=1;i<=nyueshu;i++) { if(a%yueshu[i]==0) { return 0; } } return 1; } void GetYueshu(int p) { nyueshu=0; int a=(int)sqrt((double)p); int i; for(i=2;i<=a;i++) { if(p%i==0) { nyueshu++; yueshu[nyueshu]=i; } while(p%i==0) { p=p/i; } } if(p!=1) { nyueshu++; yueshu[nyueshu]=p; } } __int64 Get(__int64 a) { if(a==0) { return 0; } __int64 sum=(a+1)*a/2; int b[20]; int num=(int)pow((double)2,(double)nyueshu)-1; int i; int j; int temp; __int64 cj; int geshu; for(i=1;i<=num;i++) { temp=i; for(j=1;j<=nyueshu;j++) { b[j]=temp%2; temp=temp/2; } cj=1; geshu=0; for(j=1;j<=nyueshu;j++) { if(b[j]==1) { cj*=yueshu[j]; geshu++; } } if(geshu%2==1) { sum-=((1+a/cj)*(a/cj)/2*cj); } else { sum+=((1+a/cj)*(a/cj)/2*cj); } } return sum; } __int64 GetValue(int l,int r,int p) { GetYueshu(p); __int64 sum=Get(r)-Get(l-1); int i; for(i=1;i<=nnode;i++) { if(l<=node[i].index&&r>=node[i].index) { if(IsCoprime(node[i].index)==1) { sum-=node[i].index; } if(IsCoprime(node[i].value)==1) { sum+=node[i].value; } } } return sum; } void Insert(int index,int value) { int i; for(i=1;i<=nnode;i++) { if(node[i].index==index) { node[i].value=value; return ; } } nnode++; node[nnode].index=index; node[nnode].value=value; } int main() { int t; int T; int i; int operate; int l; int r; int p; int index; int value; scanf("%d",&T); for(t=1;t<=T;t++) { scanf("%d%d",&n,&m); Initialize(); for(i=1;i<=m;i++) { scanf("%d",&operate); if(operate==1) { scanf("%d%d%d",&l,&r,&p); printf("%I64d\n",GetValue(l,r,p)); } else { scanf("%d%d",&index,&value); Insert(index,value); } } } return 0; }
以上是关于hdoj4407题解的主要内容,如果未能解决你的问题,请参考以下文章