[BZOJ3529] [SDOI2014] 数表 - 莫比乌斯反演 - 树状数组

原创 2016年05月30日 20:05:17

具体怎么反演的请看POPOQQQ大爷的莫比乌斯反演论文TATATAT

然后注意,当用到函数F(i)的时候再将函数筛入树状数组里。。

时间复杂度O(n log^2 n+q√n log n)

如果不自然溢出会慢很多TAT

<span style="font-size:18px;">#include "stdio.h"
#include "iostream"
#include "algorithm"
#include "stdlib.h"
#include "string.h"
#include "memory.h"
#include "math.h"

#define rep(f,a,b) for(f=a;f<=b;f++)
#define run(f,to) for(f=1;f<=to;f++)

typedef long long ll;
using namespace std;

const int N=100005,Q=20005;
inline int read(){
	int v=0; char ch=getchar();
	while(ch>'9'||ch<'0') ch=getchar();
	while(ch<='9'&&ch>='0') 
	v=v*10+ch-'0',ch=getchar();
	return v;
}

struct Query{
	int n,m,a,f,id;
	bool operator <
	(const Query&q)const{
		return a<q.a;
	}
} ask[Q];
inline bool cmp(const Query&p,const Query&q){
	return p.id<q.id;
}
struct Ans{
	int f,id;
	bool operator <
	(const Ans&q)const{
		return f<q.f;
	}
} save[N];
    
int mobius[N],pri[N],prl;
bool flag[N]; int F[N];
int val[N],sum[N],q;
const int n=100000;

void build(){
	mobius[1]=1; int i,j;
	rep(i,2,n){
		if(flag[i]^1){
			pri[++prl]=i;
			mobius[i]=-1;
		}
		for(j=1;i*pri[j]<=n;j++){
			flag[i*pri[j]]=true;
			if(i%pri[j]==0){
				mobius[i*pri[j]]=0;
				break;
			}
			mobius[i*pri[j]]=-mobius[i];
		}
	}
	rep(i,1,n){
		for (j=i;j<=n;j+=i){
			F[j]+=i;
		}
	}
	rep(i,1,n){
		save[i].f=F[i]&0x7fffffff;
		save[i].id=i;
	}
	sort(save+1,save+n+1);
}

void init(){
	int i; q=read(); rep(i,1,q){
		ask[i].n=read();
		ask[i].m=read();
		ask[i].a=read();
		ask[i].id=i;
	} sort(ask+1,ask+q+1);
}

inline int lowbit(int x){
	return x&-x;
}

inline int asksum(int l,int r){
	int ret=0;
	while(l<=r){
		ret+=val[r];
		for(r--;r-lowbit(r)>=l;r-=lowbit(r)){
			ret+=sum[r];
		}
	}
	return ret;
}

inline void add(int p,int value){
	val[p]+=value;
	for(int k=p;k<=n;k+=lowbit(k)){
		sum[k]+=value;
	}
}

void work(){
	int i,j=1,k;
	rep(i,1,q){
		while(ask[i].a>=save[j].f){
			for (k=1;k*save[j].id<=n;k++){
				add(k*save[j].id,save[j].f*mobius[k]);
			}
			j++;
		}
		int x=ask[i].n,y=ask[i].m,a=0;
		if(x>y) swap(x,y); int last;
		for (k=1;k<=x;k=last+1){
			last=min(x/(x/k),y/(y/k));
			ask[i].f+=(x/k)*(y/k)*asksum(k,last);
		}
	}
	sort(ask+1,ask+q+1,cmp);
	rep(i,1,q){
		printf("%d\n",ask[i].f&0x7fffffff);
	}
}

int main(){
	build(); init(); work();
	return 0;
}</span><span style="font-size:24px;">
</span>


版权声明:本文为博主原创文章,未经博主允许不得转载。

BZOJ 3529 SDOI2014 数表 莫比乌斯反演+树状数组

题目大意:令F(i)为i的约数和,多次询问对于1 n,m 首先如果不考虑a的限制 令g(i)为1 那么显然有 利用线性筛处理出F(i) 那么答案显然是 治好了我多年的公式恐惧...
  • PoPoQQQ
  • PoPoQQQ
  • 2014年12月22日 09:28
  • 2261

BZOJ 3930 CQOI2015 选数 莫比乌斯反演

题目见 http://pan.baidu.com/s/1o6zajc2 此外不知道H-L #include #include #include #include #includ...
  • PoPoQQQ
  • PoPoQQQ
  • 2015年04月07日 11:16
  • 18566

【莫比乌斯反演】SDOI2014 数表

【SDOI2014】数表
  • hedongnike
  • hedongnike
  • 2014年04月17日 17:38
  • 1071

bzoj 3529: [Sdoi2014]数表 莫比乌斯反演&树状数组

然后按照f[i]和a分别排序之后离线插入用树状数组维护前缀和即可。 AC代码如下: #include #include #include #define N 100005 using namespac...
  • lych_cys
  • lych_cys
  • 2016年04月04日 20:36
  • 289

【BZOJ3529】【SDOI2014】数表(莫比乌斯反演+树状数组)(好题)

Description有一张N×M的数表,其第i行第j列(1
  • mousemjq
  • mousemjq
  • 2017年03月27日 16:51
  • 322

[BZOJ3529][Sdoi2014]数表 && 莫比乌斯反演+树状数组

因为这个好好地复习了一下莫比乌斯反演 结果才发现这是自己的第一篇莫比乌斯反演 具体内容太长了 写了两页草稿纸 想看就去orz一下PoPoQQQ大爷吧 传送门http://blog.csdn.net/p...
  • shiyukun1998
  • shiyukun1998
  • 2015年04月26日 18:36
  • 346

bzoj 3529: [Sdoi2014]数表 莫比乌斯反演+树状数组

题意有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和。给定a,计算数表中不大于a的数之和。 1 < =N.m < =10^...
  • qq_33229466
  • qq_33229466
  • 2017年01月16日 11:10
  • 138

3529: [Sdoi2014]数表 莫比乌斯反演+树状数组

Dragonite神犇的题。。不知道省选现场怎么样。。我们定义F(i)F(i)为ii的约束和,题目要求的是 ∑i=1n∑j=1mF(gcd(i,j)) mod 231,F(gcd(i,j))...
  • Phenix_2015
  • Phenix_2015
  • 2016年03月03日 10:05
  • 470

[莫比乌斯反演+数状数组] BZOJ3529: [Sdoi2014]数表

题意有一张N×m的数表,其第i行第j列(1
  • CHHNZ
  • CHHNZ
  • 2017年07月03日 11:49
  • 148

【BZOJ3529】数表(莫比乌斯反演)(树状数组)(数学)

题目大意:当1≤i≤n,1≤j≤m1≤i≤n,1≤j≤m时gcd(i,j)gcd(i,j)的因数之和中不大于a的值之和...
  • can919
  • can919
  • 2017年07月11日 20:24
  • 75
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:[BZOJ3529] [SDOI2014] 数表 - 莫比乌斯反演 - 树状数组
举报原因:
原因补充:

(最多只允许输入30个字)