HDU 4407 Sum(12年金华 容斥原理)

转载请注明出处,谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526       by---cxlove 

题目:查询区间内与互质的数的和

比赛的时候SB了没看清题目,忽视了一个重要条件,初始的时候a[i]=i,而且操作不多

其实看到了也不一定会。

初始做法是:维护一段区间的LCM,然后判断LCM与P的GCD是否为1,如果是则是区间和,否则 是到子节点。

但是LCM必要太大,就只能用SET维护区间的素因子,但是素因子又多,最后又无奈的MLE。

然后又用JAVA维护前者的LCM,然后又是各种错误。

然后就没有然后了,苦逼啊

其实那个a[i]=i的条件很有用,而且操作不多,修改的位置不多

先求出l-r这个区间内与P互质的数的和,然后针对修改进行处理一下就行了。

MAP保存一下操作2就OK了。

区间和用容斥

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define inf 1<<28
#define M 100005
#define N 400005
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define pb(a) push_back(a)
#define mem(a,b) memset(a,b,sizeof(a))
#define LL long long
#define MOD 1000000007
using namespace std;
mapmp;
map::iterator it;
int n,q;
int prime[N][15]={0},flag[N]={0};
int gcd(int a,int b){
	return b==0?a:gcd(b,a%b);
}
void Prime(){
    for(int i=2;iprime[p][0]) return;
	if(num>n) return ;
	dfs(idx+1,num,cnt,m,n,p);
	dfs(idx+1,num*prime[p][idx],cnt+1,m,n,p);
}
LL slove(int n,int p){
	if(n<=0) return 0;
	ans=(LL)n*(n+1)/2;
	for(int i=1;i<=prime[p][0];i++){
	//	cout<=x&&y>=(*it).first){				
						if(gcd(c,(*it).first)==1) ret-=(*it).first;	
						if(gcd(c,(*it).second)==1) ret+=(*it).second;
					}
				}
				printf("%I64d\n",ret);;
			}
			else{
				scanf("%d%d",&x,&y);
				mp[x]=y;
			}
		}
	}
	return 0;
}



你可能感兴趣的:(ACM_数学类)