莫比乌斯&&莫比乌斯反演

最新推荐文章于 2023-11-06 15:45:12 发布
最新推荐文章于 2023-11-06 15:45:12 发布
阅读量251 收藏 5
点赞数
分类专栏： 训练日常 文章标签： 莫比乌斯
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hzaukotete/article/details/95938856
版权

莫比乌斯 懵比乌斯 懵比尤斯 梦比优斯。。。

首先入门博客：https://blog.csdn.net/codeswarrior/article/details/81541972
https://www.cnblogs.com/peng-ym/p/8647856.html

hdu 1695

欧拉＋容斥

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+10;
typedef long long ll;
ll phi[maxn];

void init(){
    phi[1] = 1;
    for(ll i = 2; i < maxn; i++){
        phi[i] = i;
    }
    for(ll i = 2; i < maxn; i++){
        if(phi[i] == i){
            for(ll j = i; j < maxn; j += i){
                phi[j] = phi[j] / i * (i - 1);
            }
        }
    }
       for(ll i = 2; i < maxn; i++){
        phi[i] += phi[i-1];
    }
}

ll owprime[maxn];
ll ownum=0;
ll solve(ll d,ll n) //求1-d中有多少个与n互质的数
{
    ownum=0;
    for(ll i = 2; i * i <= n; i++){
        if(n % i == 0){
            owprime[ownum++] = i;
            while(n % i == 0) n /= i;
        }
    }
    if(n > 1) owprime[ownum++] = n;
    ll sum=0;
    for(ll i=1; i<(1<<ownum) ;i++)
    {
        ll tmp=1,cnt=0;
        for(ll j=0;j<ownum;j++)
        {
            if(i & (1<<j))
            {
                tmp*=owprime[j];
                cnt++;
            }
        }
        if(cnt&1)
           sum+=d/tmp;
        else
           sum-=d/tmp;
    }
    return d-sum;
}

int main()
{
    ll T;
    scanf("%lld",&T);
    ll t=0;
    ll a,b,c,d,k;
    init();
    while(T--)
    {
      ll ans=0;
       scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);
      //b大
      if(k == 0 || k > b || k > d)
      {
          printf("Case %lld: 0\n",++t);
          continue;
      }
      if(b<d) swap(b,d);
      b=b/k;
      d=d/k;
      for(ll i=d+1;i<=b;i++)
      {
          ans+=solve(d,i);
      }
      printf("Case %lld: %lld\n",++t,ans+phi[d]);
    }
    return 0;
}

懵比乌斯

#include<bits/stdc++.h>
using namespace std;
#define maxn 100005
typedef long long ll;
ll isprime[maxn];
ll prime[maxn];
ll num;
ll mu[maxn];

void init()
{
    memset(isprime,false,sizeof isprime);
    mu[1]=1;
    for(ll i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            mu[i]=-1;
        }
        for(ll j=0;j<num && i*prime[j]<maxn;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                 mu[i*prime[j]]=0;
                 break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
}

int main()
{
    init();
    int T;
    int cas=0;
    scanf("%d",&T);
    while(T--)
    {
        ll a,b,c,d,k;
        scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);
        if(k==0 || k>b || k>d)
        {
           printf("Case %d: 0\n",++cas);
           continue;
        }
        ll minn=min(b,d);
        ll ans=0;
        for(ll i=1;i*k<=minn;i++)
        {
            ans+=mu[i]*(b/(i*k))*(d/(i*k));
        }
        ll tmp=0;
        for(ll i=1;i*k<=minn;i++)
        {
            tmp+=mu[i]*(minn/(i*k))*(minn/(i*k));
        }
        ans-=tmp/2;
        printf("Case %d: %lld\n",++cas,ans);
    }
    return 0;
}






#include<bits/stdc++.h>
using namespace std;
#define maxn 100005
int a,b,c,d,k;
bool isprime[maxn];
int prime[maxn];
int num;
int mu[maxn];

void init()
{
    memset(isprime,false,sizeof isprime);
    mu[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            mu[i]=-1;
        }
        for(int j=0;j<num&&i*prime[j]<maxn;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                mu[i*prime[j]]=0;
                break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
}

int main()
{
   int T;
   int cas=0;
   init();
   scanf("%d",&T);
   while(T--)
   {
       scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
       if(k==0||k>b||k>d)
       {
           printf("Case %d: 0\n",++cas);
           continue;
       }
       if(b>d) swap(b,d);
       int m=b/k;
       int n=d/k;
       long long ans=0;
       for(int i=1;i<=m;i++)
       {
           ans+=(long long)mu[i]*(m/i)*(n/i);
       }
       long long tmp=0;
       for(int i=1;i<=m;i++)
       {
           tmp+=(long long)mu[i]*(m/i)*(m/i);
       }
       ans-=tmp/2;
       printf("Case %d: %lld\n",++cas,ans);
   }
   return 0;
}

luogu 2568 https://www.luogu.org/problemnew/lists?name=2568&orderitem=pid&tag=&content=0&type=

#include<bits/stdc++.h>
using namespace std;
#define maxn 10000007
typedef long long ll;
bool isprime[maxn];
ll prime[maxn];
ll mu[maxn];
ll num;
ll n;

void init()
{
    memset(isprime,false,sizeof isprime);
    num=0;
    mu[1]=1;
    for(ll i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            mu[i]=-1;
        }
        for(ll j=0;j<num&&i*prime[j]<maxn;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                mu[i*prime[j]]=0;
                break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
}

int main()
{
    init();
    scanf("%lld",&n);
    ll ans=0;
    for(ll i=0;prime[i]<=n;i++)
    {
        ll minn=n/prime[i];
        for(ll j=1;j<=minn;j++)
        {
            ans+=mu[j]*(minn/j)*(minn/j);
        }
    }
    printf("%lld\n",ans);
    return 0;
}

hdu2841

#include<bits/stdc++.h>
using namespace std;
#define maxn 100007
typedef long long ll;
ll isprime[maxn];
ll prime[maxn];
ll num;
ll mu[maxn];

void init()
{
    memset(isprime,false,sizeof isprime);
    mu[1]=1;
    for(ll i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            mu[i]=-1;
        }
        for(ll j=0;j<num && i*prime[j]<maxn;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                mu[i*prime[j]]=0;
                break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
}

int main()
{
    init();
    int T;
    scanf("%d",&T);
    while(T--)
    {
       ll m,n;
       scanf("%lld%lld",&m,&n);
       ll ans=0;
       ll minn=min(m,n);
       for(ll i=1;i<=minn;i++)
       {
           ans+=mu[i]*(m/i)*(n/i);
       }
       printf("%lld\n",ans);
    }
    return 0;
}

spoj7001https://www.spoj.com/problems/VLATTICE/

#include<bits/stdc++.h>
using namespace std;
#define maxn 1000007
typedef long long ll;

bool isprime[maxn];
ll prime[maxn];
ll num;
ll mu[maxn];

 void init()
 {
     memset(isprime,false,sizeof isprime);
     mu[1]=1;
     for(ll i=2;i<maxn;i++)
     {
         if(!isprime[i])
         {
             prime[num++]=i;
             mu[i]=-1;
         }
         for(ll j=0;j<num && i*prime[j]<maxn;j++)
         {
             isprime[i*prime[j]]=true;
             if(i%prime[j]==0)
             {
                 mu[i*prime[j]]=0;
                 break;
             }
             else
                mu[i*prime[j]]=-mu[i];
         }
     }
 }

 int main()
 {
     int T;
     scanf("%d",&T);
     init();
     while(T--)
     {
         ll n;
         scanf("%lld",&n);
         ll ans=3;
         for(ll i=1;i<=n;i++)
         {
             ans+=mu[i]*(n/i)*(n/i)*(n/i+3);
         }
         printf("%lld\n",ans);
     }
     return 0;
 }

luogu3455

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define maxn 50005
ll isprime[maxn];
ll prime[maxn];
ll num;
ll mu[maxn];
ll sum[maxn];

void init()
{
    memset(isprime,false,sizeof isprime);
    mu[1]=1;
    for(ll i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            mu[i]=-1;
        }
        for(ll j=0;j<num && i*prime[j]<maxn;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                mu[i*prime[j]]=0;
                break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
    for(ll i=1;i<maxn;i++)
      sum[i]=sum[i-1]+mu[i];
}

int main()
{
    init();
    int T;
    scanf("%d",&T);
    while(T--)
    {
        ll a,b,d;
        scanf("%lld%lld%lld",&a,&b,&d);
        ll m=a/d;
        ll n=b/d;
        ll minn=min(m,n);
        ll r;
        ll ans=0;
        for(ll l=1;l<=minn;l=r+1)
        {
            r=min(m/(m/l),n/(n/l));
            ans+=(sum[r]-sum[l-1])*(m/l)*(n/l);
        }
        printf("%lld\n",ans);
    }
    return 0;
}

luogu2522

#include<bits/stdc++.h>
using namespace std;
#define maxn 50005
typedef long long ll;
ll isprime[maxn];
ll prime[maxn];
ll num;
ll mu[maxn];
ll sum[maxn];

void init()
{
    memset(isprime,false,sizeof isprime);
    mu[1]=1;
    for(ll i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            mu[i]=-1;
        }
        for(ll j=0;j<num && i*prime[j]<maxn;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                mu[i*prime[j]]=0;
                break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
    for(ll i=1;i<maxn;i++)
        sum[i]=sum[i-1]+mu[i];
}

ll solve(ll m,ll n,ll k)
{
   m=m/k;
   n=n/k;
   ll minn=min(m,n);
   ll ans=0;
   ll r=0;
   for(ll l=1;l<=minn;l=r+1)
   {
       r=min(n/(n/l),m/(m/l));
       ans+=(sum[r]-sum[l-1])*(m/l)*(n/l);
   }
   return ans;
}

int main()
{
    init();
    int t;
    scanf("%d",&t);
    while(t--)
    {
      ll a,b,c,d,k;
      scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);
      ll ans=0;
      ans=solve(b,d,k)-solve(b,c-1,k)-solve(a-1,d,k)+solve(a-1,c-1,k);
      printf("%lld\n",ans);
    }
    return 0;
}

https://www.luogu.org/recordnew/show/20913910

#include<bits/stdc++.h>
using namespace std;
#define maxn 1000005

int main()
{
    int n;
    scanf("%d",&n);
    int r=0;
    long long ans=0;
    for(int l=1;l<=n;l=r+1)
    {
        r=n/(n/l);
        ans+=(r-l+1)*(n/l);
    }
    printf("%lld\n",ans);
    return 0;
}

https://www.luogu.org/problemnew/show/P3935

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define mod 998244353

ll solve(ll n)
{
    ll r=0;
    ll ans=0;
    for(ll l=1;l<=n;l=r+1)
    {
        r=n/(n/l);
        ans+=(r-l+1)*(n/l);
    }
    return ans;
}

int main()
{
    ll l,r;
    scanf("%lld%lld",&l,&r);
    ll ans=0;
    ans=(solve(r)%mod-solve(l-1)%mod+mod)%mod;
    printf("%lld\n",ans);
    return 0;
}

https://www.luogu.org/problemnew/solution/P2257

#include<bits/stdc++.h>
using namespace std;
#define maxn 10000000
typedef long long ll;
bool isprime[maxn];
int prime[maxn];
int num;
int mu[maxn];
int sum[maxn];
int lca[maxn];
inline void read(int &x)
{
    x=0;
    static int p;p=1;
    static char c;c=getchar();
    while(!isdigit(c)){if(c=='-')p=-1;c=getchar();}
    while(isdigit(c)) {x=(x<<1)+(x<<3)+(c-48);c=getchar();}
    x*=p;
}
inline void print(long long x)
{
    static int cnt;
    static int a[15];
    cnt=0;
    do
    {
        a[++cnt]=x%10;
        x/=10;
    }while(x);
    for(int i=cnt;i>=1;i--)putchar(a[i]+'0');
    puts("");
}
void init()
{
    memset(isprime,false,sizeof false);
    mu[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            mu[i]=-1;
        }
        for(int j=0;j<num && i*prime[j]<maxn;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                mu[i*prime[j]]=0;
                break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
    for(int i=0;i<num;i++)
        for(ll j=1;j*prime[i]<maxn;j++)
          lca[j*prime[i]]+=mu[j];
    for(int i=1;i<maxn;i++)
       sum[i]=sum[i-1]+lca[i];
}

int n,m;
int main()
{
    init();
    int T;
    scanf("%d",&T);
    while(T--)
    {
        read(n);
        read(m);
        if(n>m) n^=m^=n^=m;
        int r=0;
        ll ans=0;
        for(int l=1;l<=n;l=r+1)
        {
            r=min(n/(n/l),m/(m/l));
            ans+=1LL*(n/l)*(m/l)*(sum[r]-sum[l-1]);
        }
        print(ans);
    }
    return 0;
}

https://www.luogu.org/problemnew/show/P3327

题解：https://www.cnblogs.com/peng-ym/p/8667321.html

#include<bits/stdc++.h>
using namespace std;
#define maxn 50005
typedef long long ll;
bool isprime[maxn];
int prime[maxn];
int num;
int mu[maxn];
int sum[maxn];
ll g[maxn];

void init()
{
    memset(isprime,false,sizeof isprime);
    mu[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            mu[i]=-1;
        }
        for(int j=0;j<num && i*prime[j]<maxn;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                mu[i*prime[j]]=0;
                break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
    for(int i=1;i<maxn;i++)
        sum[i]=sum[i-1]+mu[i];
    for(int i=1;i<maxn;i++)
    {
        ll ans=0;
        for(int l=1,r;l<=i;l=r+1)
        {
            r=i/(i/l);
            ans+=1ll*(r-l+1)*1ll*(i/l);
        }
        g[i]=ans;
    }
}

int main()
{
    init();
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int n,m;
        scanf("%d%d",&n,&m);
        int minn=min(n,m);
        ll ans=0;
        for(int l=1,r;l<=minn;l=r+1)
        {
            r=min(n/(n/l),m/(m/l));
            ans+=(sum[r]-sum[l-1])*1ll*g[n/l]*1ll*g[m/l];
        }
        printf("%lld\n",ans);
    }
    return 0;
}

P1447 [NOI2010]能量采集 题解

#include<bits/stdc++.h>
using namespace std;
#define maxn 100005
bool isprime[maxn];
int prime[maxn];
int phi[maxn];
int n,m;
int num;
long long sum[maxn];


void init()
{
    phi[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!isprime[i])
        {
            prime[num++]=i;
            phi[i]=i-1;
        }
        for(int j=0;j<num && i*prime[j]<maxn ;j++)
        {
            isprime[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=prime[j]*phi[i];
                break;
            }
            else
                phi[i*prime[j]]=(prime[j]-1)*phi[i];
        }
    }
    for(int i=1;i<maxn;i++)
      sum[i]=sum[i-1]+phi[i];
}

int main()
{
    init();
    scanf("%d%d",&n,&m);
    if(n>m) swap(n,m);
    long long ans=0;
    for(int l=1,r;l<=n;l=r+1)
    {
        r=min(n/(n/l),m/(m/l));
        ans+=1ll*(n/l)*(m/l)*(sum[r]-sum[l-1]);
    }
    ans=2*ans-1ll*n*m;
    printf("%lld\n",ans);
    return 0;
}

 

hzaukotete
关注
专栏目录
莫比乌斯反演入门
mumei的博客
11-06 929
这个算是竞赛中用到的组合数学里最多的知识点，也是比较难理解的，在被将近一个星期的折磨后总算是理解了一些。 其实莫比乌斯反演就只有两个重要的公式，而且用到的最多的也只有一个（其实都差不多），但是需要一些预备知识。 莫比乌斯函数，质数线性筛，整除分块，积性函数以及一些其他的数论知识，下面会将其中一部分重要的算法， 目录 1.莫比乌斯函数 2.莫比乌斯函数的线性筛 3.整除分块 4....
欧拉函数&莫比乌斯反演
Zeyu_King的专栏
02-01 1990
近几天做了几道有关反演的问题，再次集合一下吧。 1、[BZOJ 2301]HAOI2011 Problem b 2、[BZOJ 2440]中山市选2011 完全平方数 3、gcd 4、[BZOJ 2186]SDOI2008 莎拉公主的困惑 5、[BZOJ 3529]SDOI2014 数表 （蒟蒻自认为反演一类的的题目重要的就是记住两个重要的公式：1、sigma(mu[i] ,
HDU 5608 莫比乌斯反演 + 莫比乌斯函数前缀和
退役ACMer
08-06 1249
传送门: HDU 5608 题解: 令 G（n） = n ^ 2 - 3 * n + 2 先反演得:f(n)=∑d|nG(d)∗u(nd)f(n) = \sum_{d|n}G(d) * u(\frac{n}{d})令: A(n)=∑i=1nf(i)A(n) = \sum_{i = 1}^nf(i) 则:A(n)=∑i=1nf(i)=∑i=1n∑d|iG(d)∗u(id)A(n) = \sum_
莫比乌斯函数&莫比乌斯反演
qq_15647239的博客
05-09 284
莫比乌斯反演：莫比乌斯反演：莫比乌斯反演： f(n)=∑d∣ng(d)⇒g(n)=∑d∣nf(nd)∗μ(d)f(n)=\sum^{}_{d|n}g(d)\Rightarrow g(n)=\sum^{}_{d|n}f(\frac{n}{d})*\mu(d)f(n)=d∣n∑​g(d)⇒g(n)=d∣n∑​f(dn​)∗μ(d) 莫比乌斯函数：莫比乌斯函数：莫比乌斯函数： (1)若d=1,则μ(d)...
莫比乌斯反演 & 二项式反演
cqbzlydd的博客
01-12 289
感觉莫比乌斯反演的可操作性不大。 一般地： ∑d∣nμ(d)=[n=1]\sum_{d|n}\mu(d)=[n=1]∑d∣n​μ(d)=[n=1] - 公式 1 欧拉反演： ∑d∣nϕ(d)=n\sum_{d|n}\phi(d)=n∑d∣n​ϕ(d)=n - 公式 2 例一：求 ∑i=1ngcd⁡(i,n)\sum_{i=1}^n\gcd(i,n)∑i=1n​gcd(i,n) sol：原式 = ∑j∣nnjϕ(j)\sum_{j|n}\frac{n}{j}\phi(j)∑j∣n​jn​ϕ(j) （代入公式
莫比乌斯反演基本类型
欠我半块小饼干的博客
11-01 287
背景 ​ 遇到过的莫比乌斯反演，除了狄雷克利卷积没看懂没做过外，我将遇到过的莫比乌斯反演大致分成以下三类。 三大莫比乌斯套路类型 1、i∈[1,n],j∈[1,m],求gcd(i,j)==1的数对个数。(n<=m) ​ ∑i=1n∑j=1m[gcd(i,j)==1]=∑i=1nμ(i)⌊ni⌋⌊mi⌋\sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)==1]=\sum_{i=1}^nμ(i)\lfloor \frac{n}{i} \rfloor \lfloor \frac{m}{i}
莫比乌斯反演&整数分块
I_have_a_world的博客
09-17 275
目录前言参考资料莫比乌斯反演预备知识：整数分块（感觉就是一个思想）总结莫比乌斯函数μ\muμ（数论函数也可以视作一个数列）总结莫比乌斯反演题目训练 前言 参考资料 peng-ym-莫比乌斯反演 莫比乌斯反演 预备知识：整数分块（感觉就是一个思想） 参考博客：peng-ym-整除分块 总结 一般形式：∑i=1n[n/i]\sum_{i=1}^{n}[n/i]∑i=1n​[n/i] 过程：r=n/(n/l)，然后不断枚举l，求得r，[l,r]部分的值是一样的，不必遍历，直接得到长度*数。 复杂度：由
莫比乌斯函数&莫比乌斯反演证明
weixin_30872499的博客
09-13 461
前言？ 终于放假了~~感觉再不趁机颓会儿我博客就废了…… 赶紧写点东西刷刷存在感（骗点积分） 莫比乌斯函数 定义一种函数$\mu(d)$，满足： 1.若$d=1$，则$\mu(d)=1$。 2.若$d=p_1p_2p_3\cdots p_k$且$p_i$为互异素数时，$\mu(d)=(-1)^k$。 3.其他情况$\mu(d)=0$。 那么我们将函数$\mu(d)$称作...
容斥原理，莫比乌斯反演
最新发布
wanglianda2007的博客
11-06 163
浅谈莫比乌斯反演在计数问题中的一些应用。
莫比乌斯反演—详解
XianHaoMing的博客
05-27 8356
介绍 1、莫比乌斯反演是组合数学中很重要的内容，可以用于解决很多组合数学的问题。 2、莫比乌斯反演是数论中的重要内容，在许多情况下能够简化运算。 3、是个个很神奇的东西。 引入 考虑以下求和函数 fn=∑d|ngdf_n=\sum_{d|n}g_d 那么根据定义我们可以知道 f1f_1=g1g_1 f2f_2=g1g_1+g2g_2 f3f_3=g1g_1+g3g_3
莫比乌斯反演
05-18
### 莫比乌斯反演详解 #### 一、莫比乌斯反演概念与基本公式 **莫比乌斯反演**是一种在数论领域内常用的技巧，主要用于简化某些涉及数论函数的计算问题。它通过一种特定的方式转换问题的形式，从而达到简化计算的...
莫比乌斯反演公式 Möbius inversion formula
10-16
### 莫比乌斯反演公式的详细解析 #### 引言 在数学领域，尤其是在数论中，莫比乌斯反演公式扮演着极其重要的角色。它由德国数学家奥古斯特·费迪南德·莫比乌斯于1832年引入，并在后续的发展中得到了广泛的推广与...
spoj4491 莫比乌斯反演
08-22
### spoj4491 莫比乌斯反演 #### 题目解析 题目要求计算在给定的两个正整数\( n \)和\( m \)范围内，有多少对\((a, b)\)使得它们的最大公约数\(\text{gcd}(a, b) = d\)，其中\(d\)为素数，并且\(1 \leq a \leq n\)...
搜索进阶-迭代加深搜索
热门推荐
hzaukotete的博客
07-26 1万+
我个人对这个搜索的理解就是以BFS的思想写DFS。 具体来说就是，首先深度优先搜索k层，若没有找到可行解，再深度优先搜索k+1层，直到找到可行解为止。由于深度是从小到大逐渐增大的，所以当搜索到结果时可以保证搜索深度是最小的。这也是迭代加深搜索在一部分情况下可以代替广度优先搜索的原（还比广搜省空间）。   前提： 题目一定要有解，否则会无限循环下去。   好处： 1.时间复杂度只比BF...
DFS（深度优先搜索）+剪枝
hzaukotete的博客
07-23 7657
                                       DFS(深度优先搜索) 引入：为什么要搜索解决问题？因为很多问题并不可以直接计算出结果，需要遍历所有情况才能发找到最优解或者可行解。 深度优先搜索： 采用递归的方法，先沿一条路搜到底，在递归回上一个节点，沿另一个方向搜索，以此类推。 类似于树的先根遍历; 思想：一直往深处走，直到找到解或者走不下去为止 模板：...
二分法/二分思想（二分答案法）
hzaukotete的博客
07-24 2462
解题方法： 二分答案法 适用题型： 求可行解的最大值或最小值问题 解题步骤： 1.确定答案的最大值和最小值 2.判断二分所得值是否满足条件 3.可行解必须具有单调性（当k可行时，k+1可行或者k-1可行） 注意事项： 每次都要确保right和left有一个被改变。 模板： //可行解最大值问题 int l=minn,r=maxn; int ans=-1; while(l&l...
组合数学：棋盘的完美覆盖问题
hzaukotete的博客
09-21 2014
                                                 完美覆盖问题 结论：m*n棋盘被b格牌完美覆盖的充分必要条件是：b是m或则n的一个因子...
01BFS
hzaukotete的博客
10-17 814
01BFS 先欠着  
二次剩余（模板）
hzaukotete的博客
08-16 696
#include<bits/stdc++.h> using namespace std; typedef long long ll; ll w; struct num{ ll x,y; }; num mul(num a,num b,ll p) { num ans={0,0}; ans.x=((a.x*b.x%p+a.y*b.y%p*w%p)%p+p)%p; ...
莫比乌斯反演详解与应用
“莫比乌斯反演-莫比乌斯反演经典讲稿” 莫比乌斯反演是一种在数论中用于解决计数问题的技巧，它由德国数学家莫比乌斯提出。在数论和组合数学中，它常常用于简化计算，特别是在处理关于数的乘积的问题时。讲稿来自...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值