2393: Cirno的完美算数教室容斥原理

最新推荐文章于 2018-05-02 22:49:30 发布

ws_fqk

最新推荐文章于 2018-05-02 22:49:30 发布

阅读量724

点赞数

CC 4.0 BY-SA版权

分类专栏： My Code

本文链接：https://blog.youkuaiyun.com/Phenix_2015/article/details/50630110

My Code 专栏收录该内容

349 篇文章

订阅专栏

本文介绍了一个使用容斥原理解决特定数学问题的C++程序实例。通过递归生成特定数字序列，并运用最大公约数计算及深度优先搜索来求解问题答案。该程序实现了对指定区间内满足条件的数字进行计数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

这题好像做过一次。。想复习一下容斥原理。。
qwq其实自己随便yy一下就好了。。没什么好复习的。。

#include<iostream>
#include<cstdio>
#include<algorithm>
#define ll long long 
using namespace std;
int l,r,cnt,tot;
ll ans;
ll a[2005],b[2005];
bool vis[2005];
inline void pre(ll x)
{
    if (x>r) return;
    a[++cnt]=x;
    pre(x*10+2); pre(x*10+9);
}
ll gcd(ll a,ll b)
{
    return b==0?a:gcd(b,a%b);
}
void dfs(int x,int y,ll z)
{
    if (x>tot)
    {
        if (y&1) ans+=r/z-(l-1)/z;
        else if (y) ans-=r/z-(l-1)/z;
        return;
    }
    dfs(x+1,y,z);
    ll tmp=b[x]*z/gcd(b[x],z);
    if (tmp<=r) dfs(x+1,y+1,tmp);
}
int main()
{
    scanf("%d%d",&l,&r);
    pre(2); pre(9); 
    sort(a+1,a+cnt+1);
    for (int i=1;i<=cnt;i++)
        if (!vis[i])
        {
            b[++tot]=a[i];
            for (int j=i+1;j<=cnt;j++)
                if (!(a[j]%a[i])) vis[j]=1;
        }
    dfs(1,0,1);
    cout << ans << endl;
    return 0;
}