关闭

BZOJ 1853 [Scoi2010]幸运数字

标签: BZOJ-1853
105人阅读 评论(0) 收藏 举报
分类:

容斥原理。

[l,r]6S18S266S368S4......

21+22+...+210=2112=2046,

66666666943

搜索的时候从数字大的数开始搜索,以减少回溯层数,否则TLE。

求LCM的时候可能爆longlong,要用double- -

时间复杂度O(玄学)…

#include<cstdio>
#include<algorithm>
using namespace std;
const double eps = 1e-3;
int m;
long long l, r, a[3000], ans;
bool vis[3000];
void pre(long long x)
{
    if(x>r)return;
    a[++m]=x;
    pre(x*10+6);
    pre(x*10+8);
}
bool cmp(long long a,long long b){return a>b;}
long long gcd(long long a, long long b){return b==0?a:gcd(b,a%b);}
void dfs(int now, int cnt, long long x)
{
    if(now>m)
    {
        if(cnt&1)ans+=r/x-(l-1)/x;
        else if(cnt)ans-=r/x-(l-1)/x;
        return;
    }
    long long tmp=x/gcd(a[now],x);
    dfs(now+1,cnt,x);
    if((double)a[now]*tmp<=r+eps)
        dfs(now+1,cnt+1,tmp*a[now]);
}
int main()
{
    scanf("%lld%lld",&l,&r);
    pre(6);pre(8);
    sort(a+1,a+1+m,cmp);
    for(int i = m; i >= 1; i--)
        if(!vis[i])
            for(int j = i-1; j >=1; j--)    
                if(a[j]%a[i]==0)
                    vis[j]=1;
    int tmp=0;
    for(int i = 1; i <= m; i++)
        if(!vis[i])
            a[++tmp]=a[i];
    m=tmp;
    dfs(1,0,1);
    printf("%lld\n",ans);
}
1
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:34940次
    • 积分:3291
    • 等级:
    • 排名:第10275名
    • 原创:307篇
    • 转载:0篇
    • 译文:0篇
    • 评论:22条
    公告
    欢迎您 !
    由于博主是一个蒟蒻,
    如有错误之处还望指出!
    o( =•ω•= )o
    鸡汤
    ❤不要轻易否定自己,不要轻易放弃。
    ❤上天撒下一缕光明,我便不应安逸于黑暗之中。
    ❤不到万不得已,为什么要放弃。
    神犇们的博客
    文章分类