Bzoj 1853: [Scoi2010]幸运数字 容斥原理,深搜

1853: [Scoi2010]幸运数字

Time Limit: 2 Sec  Memory Limit: 64 MB
Submit: 1774  Solved: 644
[Submit][Status][Discuss]

Description

在中国，很多人都把6和8视为是幸运数字！lxhgww也这样认为，于是他定义自己的“幸运号码”是十进制表示中只包含数字6和8的那些号码，比如68，666，888都是“幸运号码”！但是这种“幸运号码”总是太少了，比如在[1,100]的区间内就只有6个（6，8，66，68，86，88），于是他又定义了一种“近似幸运号码”。lxhgww规定，凡是“幸运号码”的倍数都是“近似幸运号码”，当然，任何的“幸运号码”也都是“近似幸运号码”，比如12，16，666都是“近似幸运号码”。 现在lxhgww想知道在一段闭区间[a, b]内，“近似幸运号码”的个数。

Input

输入数据是一行，包括2个数字a和b

Output

输出数据是一行，包括1个数字，表示在闭区间[a, b]内“近似幸运号码”的个数

Sample Input

【样例输入1】
1 10
【样例输入2】
1234 4321

Sample Output

【样例输出1】
2
【样例输出2】
809

HINT

【数据范围】
对于30%的数据，保证1 < =a < =b < =1000000
对于100%的数据，保证1 < =a < =b < =10000000000

Source

Day1

题解：

容斥原理＋爆搜。

记得开double或unsigned long long。

 #include<bits/stdc++.h>
 using namespace std;
 #define LL unsigned long long
 LL sum,cc[],lc,a,b;
 bool vis[];
 LL Gcd(LL aa,LL bb){if(bb==)return aa;else return Gcd(bb,aa%bb);}
 void DFS(LL x,LL y,LL gs,LL lcm, LL gcd,LL ii)
 {
     LL LCM,GCD,i;
     //if(lcm>x||lcm>y)return;
     if(lcm>x&&lcm>y)return;
     if(gs>)
     {
         if(gs%!=)sum+=((LL)(y/lcm)-(LL)(x/lcm));
         else sum-=((LL)(y/lcm)-(LL)(x/lcm));
     }
     if(gs==lc+||ii+>lc)return;
     for(i=ii+;i<=lc;i++)
     {
         if(vis[i]==false)
         {
             vis[i]=true;
             LCM=lcm;GCD=gcd;
             gcd=Gcd(lcm,cc[i]);lcm=(lcm*cc[i])/gcd;
             DFS(x,y,gs+,lcm,gcd,i);
             lcm=LCM;gcd=GCD;
             vis[i]=false;
         }
     }
 }
 LL calc(LL x,LL y)
 {
     /*for(i=1;i<=lc;i++)
     {
         tot=0;
         DFS(x);
     }*/
     memset(vis,false,sizeof(vis));sum=;
     DFS(x-,y,,,,);
     return sum;
 }
 void dfs(LL k)
 {
     if(k>b)return;
     if(k!=)cc[++lc]=k;
     dfs(k*+);
     dfs(k*+);
 }
 int main()
 {
     LL len,i;
     scanf("%lld %lld",&a,&b);
     memset(cc,,sizeof(cc));lc=;
     dfs();
     sort(cc+,cc+lc+);
     len=unique(cc+,cc+lc+)-(cc+);
     for(i=;i<=len/;i++)swap(cc[i],cc[len-i+]);
     lc=len;
     printf("%lld",calc(a,b));
     return ;
 }