BZOJ 3944 Sum 数论

题目大意:求∑ni=1φ(i)和∑ni=1μ(i)。

n≤231?1

令F(n)为f(n)的前缀和,G(n)为g(n)的前缀和,且满足g(n)=∑i|nf(i),则有:

G(n)=∑ni=1g(i)

=∑ni=1∑j|if(j)

=∑nj=1∑j|if(j)

=∑nj=1?nj?f(j)

=∑nj=1F(?nj?)

∴F(n)=G(n)?∑ni=2F(?ni?)

若G(n)可以在O(1)时间内算出,则F(n)可以在O(n34)的时间内算出(别问我复杂度是咋算的)

如果预处理前O(n23)的部分,则F(n)可以在O(n23)的时间内算出(别问我复杂度是咋算的*2)

我作死多写了个log出来。。。卡死了。。。

#include <map>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 5000000
using namespace std;
int n;
int prime[M/5],tot;
long long phi[M],mu[M];
bool not_prime[M];
map<int,long long> _phi,_mu;
void Linear_Shaker()
{
    long long i,j;
    phi[1]=1;mu[1]=1;
    for(i=2;i<M;i++)
    {
        if(!not_prime[i])
        {
            phi[i]=i-1;
            mu[i]=-1;
            prime[++tot]=i;
        }
        for(j=1;prime[j]*i<M;j++)
        {
            not_prime[prime[j]*i]=true;
            if(i%prime[j]==0)
            {
                phi[prime[j]*i]=phi[i]*prime[j];
                mu[prime[j]*i]=0;
                break;
            }
            phi[prime[j]*i]=phi[i]*(prime[j]-1);
            mu[prime[j]*i]=-mu[i];
        }
    }
    for(i=1;i<M;i++)
    {
        phi[i]+=phi[i-1];
        mu[i]+=mu[i-1];
    }
}
long long Calculate_Phi(long long n)
{
    map<int,long long>::iterator it;
    if(n<M)
        return phi[n];
    if((it=_phi.find(n))!=_phi.end())
        return it->second;
    long long i,last,re=(long long)n*(n+1)>>1;
    for(i=2;i<=n;i=last+1)
    {
        last=n/(n/i);
        re-=(last-i+1)*Calculate_Phi(n/i);
    }
    return _phi[n]=re;
}
long long Calculate_Mu(long long n)
{
    map<int,long long>::iterator it;
    if(n<M)
        return mu[n];
    if((it=_mu.find(n))!=_mu.end())
        return it->second;
    long long i,last,re=1;
    for(i=2;i<=n;i=last+1)
    {
        last=n/(n/i);
        re-=(last-i+1)*Calculate_Mu(n/i);
    }
    return _mu[n]=re;
}
int main()
{
    int T;
    Linear_Shaker();
    for(cin>>T;T;T--)
    {
        scanf("%d",&n);
        printf("%lld %lld\n",Calculate_Phi(n),Calculate_Mu(n));
    }
}
时间: 2024-10-16 23:31:42

BZOJ 3944 Sum 数论的相关文章

●杜教筛入门(BZOJ 3944 Sum)

入门杜教筛啦. http://blog.csdn.net/skywalkert/article/details/50500009(好文!) 可以在$O(N^{\frac{2}{3}})或O(N^{\frac{3}{4}})$的复杂度内解决求某些数论函数f(n)(或f的前缀和S(n)$)的值. 先来看看原理是什么.(接下来推导如何求数论函数f(n)的前缀和S(n)) 现在有两个数论函数$f( )和g( )$ (同时定义f的前缀和函数$S(n)=\sum_{i=1}^{n}f(i)$) 有狄利克雷乘

BZOJ.3944.Sum(杜教筛)

题目链接 又写个模板题用了半晚上.. 卡常写法 还非常短.. //35332 kb 7608 ms //跟Kelin dalao学一波卡常.怎么还是很慢QAQ //phi[]要longlong! #include <cstdio> #include <cstring> #include <algorithm> typedef long long LL; const int N=2e6; int cnt,Max,P[N>>3]; LL phi[N],mu[N]

bzoj 3944 Sum —— 杜教筛

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3944 杜教筛入门题! 看博客:https://www.cnblogs.com/zjp-shadow/p/8491542.html 写法模仿其他博客的,但很慢啊... 代码如下: #include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<ma

3944: Sum[杜教筛]

3944: Sum Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 3471  Solved: 946[Submit][Status][Discuss] Description Input 一共T+1行 第1行为数据组数T(T<=10) 第2~T+1行每行一个非负整数N,代表一组询问 Output 一共T行,每行两个用空格分隔的数ans1,ans2 Sample Input 6 1 2 8 13 30 2333 Sample Output 1 1

3944: Sum(杜教筛)

3944: Sum Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 4930  Solved: 1313[Submit][Status][Discuss] Description Input 一共T+1行 第1行为数据组数T(T<=10) 第2~T+1行每行一个非负整数N,代表一组询问 Output 一共T行,每行两个用空格分隔的数ans1,ans2 Sample Input 6 1 2 8 13 30 2333 Sample Output 1 1

BZOJ 1257: [CQOI2007]余数之和sum( 数论 )

n >= k 部分对答案的贡献为 k * (n - k) n < k 部分贡献为 ∑ (k - ⌊k / i⌋ * i)  = ∑  , ⌊k / i⌋ 相等的数是连续的一段, 此时这段连续的数对答案的贡献成等差数列, 可以O(1)求出..然后就分⌊k / i⌋ 相等的一块一块来就行了. 分出来大概是sqrt(k)块.这个sqrt(k)我并不会证Orz...写了个程序验证了一下, 分出来的块数和2 * sqrt(k)非常接近. 所以时间复杂度为O(sqrt(k)) --------------

Sum BZOJ 3944

Sum [问题描述] 给定一个正整数 N ( N <= 231 - 1 ) 求: [输入格式] 一共T+1行 第1行为数据组数T(T<=10) 第2~T+1行每行一个非负整数N,代表一组询问 [输出格式] 一共T行,每行两个用空格分隔的数ans1,ans2 [样例输入] 6 1 2 8 13 30 2333 [样例输出] 1 1 2 0 22 -2 58 -3 278 -3 1655470 2 题解: 主要算法:杜教筛: 首先推一波式子

BZOJ 2142 礼物 数论

这道题是求组合数终极版. C(n,m) mod P n>=1e9 m>=1e9 P>=1e9且为合数且piqi<=1e5 拓展lucas定理. 实际上就是一点数论小知识的应用. 这篇文章对于CRT和lucas定理的学习非常不错. #include<bits/stdc++.h> using namespace std; #define ll long long #define FILE "dealing" #define up(i,j,n) for(i

bzoj 2818 GCD 数论 欧拉函数

bzoj[2818]Gcd Description 给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对. Input 一个整数N Output 如题 Sample Input 4 Sample Output 4 HINT hint对于样例(2,2),(2,4),(3,3),(4,2) 1<=N<=10^7 题解一(自己yy) phi[i]表示与x互质的数的个数 即gcd(x,y)=1 1<=y<x ∴对于x,y 若a为素数 则gcd(xa,