bzoj1968(数论/欧拉筛/筛法)

这题貌似有许多解法。。一个个尝试。。

解法一：欧拉筛。。做这题最主要还是要这个模板。。不然我是不会做数论题的。。看起来貌似是o(n)？

只记板子不作深究。。

解法二：普通筛，欧拉筛或许筛素数可以到o(n)，然而个人筛因子的复杂度感觉也不特别好的样子。。用普通筛试了一下竟然还过了。。虽然差点超时。。复杂度o(nlogn)

解法三：对这题来说如果转念一想参考每个因子的贡献的话，每个因子的贡献就是n/i了，所以答案为sigma(n/i)，不过同样是o(n)和欧拉筛差得挺大的。。

解法四：在解法三的基础上把n/i列出来，发现有许多段是重复的，所以将重复的压缩之后再求，具体就是将值n/i右端点给求出来，为n/(n/i)，emmmm这个为什么。。当t=(double)n/(int)(n/i)时，n/t恰好为n/i，所以如果大于t，用实数除出来是比n/i小的，又取了整，所以比n/i小。。这样求复杂度为 o(sqrt(n))

数学是个有趣的东西。。可惜自己太弱了qaq

解法一：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<list>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define eps 1e-8
#define inf 1e9
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define ls T[i<<1]
#define rs T[i<<1|1]
#define op T[i]
#define mid (x+y>>1)
#define NM 1000005
#define nm 400005
#define pi 3.141592653
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}

int n,d[NM],prime[NM],t[NM],tot,ans;
bool v[NM];
void Euler(){
	d[1]=1;
	inc(i,2,n){
		if(!v[i]){
			prime[++tot]=i;
			d[i]=2;t[i]=1;
		}
		inc(j,1,tot){
			if(i*prime[j]>n)break;
			v[i*prime[j]]=true;
			if(i%prime[j]==0){
				t[i*prime[j]]=t[i]+1;
				d[i*prime[j]]=d[i]/(t[i]+1)*(t[i*prime[j]]+1);
				break;
			}else{
				d[i*prime[j]]=d[i]*2;
				t[i*prime[j]]=1;
			}
		}
	}
}


int main(){
	n=read();
	Euler();
	inc(i,1,n)ans+=d[i];
	return 0*printf("%d\n",ans);
}

解法二：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<list>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define eps 1e-8
#define inf 1e9
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define ls T[i<<1]
#define rs T[i<<1|1]
#define op T[i]
#define mid (x+y>>1)
#define NM 1000005
#define nm 400005
#define pi 3.141592653
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}






int n,a[NM];

int main(){
	n=read();
	inc(i,1,n)for(int j=1;i*j<=n;j++)a[i*j]++;
	inc(i,1,n)a[i]+=a[i-1];
	printf("%d\n",a[n]);
	return 0;
}

解法三：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<list>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define eps 1e-8
#define inf 1e9
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define ls T[i<<1]
#define rs T[i<<1|1]
#define op T[i]
#define mid (x+y>>1)
#define NM 1000005
#define nm 400005
#define pi 3.141592653
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}






int n,ans;
int main(){
	n=read();
	inc(i,1,n)ans+=n/i;
	return 0*printf("%d\n",ans);
}

解法四：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<list>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define eps 1e-8
#define inf 1e9
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define ls T[i<<1]
#define rs T[i<<1|1]
#define op T[i]
#define mid (x+y>>1)
#define NM 1000005
#define nm 400005
#define pi 3.141592653
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}






int n,ans;
int main(){
	n=read();
	for(int i=1,j;i<=n;i=j+1){
		j=n/(n/i);
		ans+=n/i*(j-i+1);
	}
	return 0*printf("%d\n",ans);
}

1968: [Ahoi2005]COMMON 约数研究

Time Limit: 1 Sec   Memory Limit: 64 MB
Submit: 3350   Solved: 2285
[ Submit][ Status][ Discuss]

Description

Input

只有一行一个整数 N（0 < N < 1000000）。

Output

只有一行输出，为整数M，即f(1)到f(N)的累加和。

Sample Input

3

Sample Output

5

HINT

Source

Day2

[ Submit][ Status][ Discuss]