幸运数字(容斥原理+高精度与低精度的加 减 除)

幸运数字

题目描述：
LYK 最近运气很差，例如在 NOIP 初赛中仅仅考了 90 分，刚刚卡进复赛，于是它决定使用一些方法来增加自己的运气值。
它觉得，通过收集幸运数字可以快速的增加它的 RP 值。
它给幸运数字下了一个定义：如果一个数 x 能被 3 整除或被 5 整除或被 7 整除，则这个数为幸运数字。
于是它想让你帮帮它在 L~R 中存在多少幸运数字。
输入格式：
第一行两个数 L,R。
输出格式：
一个数表示答案。
输入样例
10 15
输出样例
4
数据范围：
对于 50%的数据 1<=L<=R<=10^5。
对于 60%的数据 1<=L<=R<=10^9。
对于 80%的数据 1<=L<=R<=10^18。
对于 90%的数据 1<=L<=R<=10^100。
对于另外 10%的数据 L=1， 1<=R<=10^100。
对于 100%的数据 L， R 没有前导 0。
思路：
容斥原理
设num(n,i)表示1~n中是i的倍数
则求1~n中3或5或7的倍数为
num(n,3)+num(n,5)+num(n,7)-num(n,3*5)-num(n,3*7)-num(5*7)+num(n,3*5*7)
这样求出1~r的减去1~l-1的即可
高精度！！！

#include<iostream>#include<cstdio>#include<cstring>using namespace std;const int maxn=1010;int len,a[maxn],b[maxn],t[maxn],ans[maxn];int sa[maxn],sb[maxn],c[2]={1,1};char s1[maxn],s2[maxn];void jian2(int *a,int *b){    memset(sa,0,sizeof(sa));    memset(sb,0,sizeof(sb));    sa[0]=a[0];    for(int i=1;i<=a[0];i++)    sa[i]=a[a[0]-i+1];    for(int i=1;i<=b[0];i++)    sb[i]=b[b[0]-i+1];    int len=1;    while(len<=sa[0])    {        sa[len]=sa[len]-sb[len];        if(sa[len]<0)        {            sa[len]=sa[len]+10;            sa[len+1]--;        }        len++;    }    while(sa[0]>1&&!sa[a[0]])    sa[0]--;a[0]=sa[0];    for(int i=1;i<=sa[0];i++)    a[i]=sa[sa[0]-i+1];}void div(int *r,int *a,int x){    int s=0;r[0]=a[0];    for(int i=1;i<=a[0];i++)    {        s=s*10+a[i];        r[i]=s/x;        s=s%x;    }}void add(int *a,int *b){    memset(sa,0,sizeof(sa));    memset(sb,0,sizeof(sb));    int len1=a[0],len2=b[0],len=1;    for(int i=1;i<=len2;i++)    sb[i]=b[len2-i+1];    while(len<=len1||len<=len2)    {        sa[len]+=a[len]+sb[len];        sa[len+1]=sa[len]/10;        sa[len]=sa[len]%10;        len++;    }    while(len>1&&!sa[len])    len--;sa[0]=len;    for(int i=0;i<=sa[0];i++) a[i]=sa[i];}void jian(int *a,int *b){    memset(sb,0,sizeof(sb));    for(int i=1;i<=b[0];i++)    sb[i]=b[b[0]-i+1];    int len=1;    while(len<=a[0])    {        a[len]=a[len]-sb[len];        if(a[len]<0)        {            a[len]=a[len]+10;            a[len+1]--;        }        len++;    }    while(a[0]>1&&!a[a[0]])    a[0]--;}int main(){    freopen("number.in","r",stdin);    freopen("number.out","w",stdout);    cin>>s1>>s2;    int l1=strlen(s1),l2=strlen(s2);    a[0]=l1,b[0]=l2;    for(int i=1;i<=l1;i++) a[i]=s1[i-1]-'0';     for(int i=1;i<=l2;i++) b[i]=s2[i-1]-'0';    jian2(a,c);    div(t,b,3);add(ans,t);    div(t,b,5);add(ans,t);    div(t,b,7);add(ans,t);    div(t,a,3*5);add(ans,t);    div(t,a,3*7);add(ans,t);    div(t,a,5*7);add(ans,t);    div(t,b,3*5*7);add(ans,t);    div(t,b,3*5);jian(ans,t);    div(t,b,3*7);jian(ans,t);    div(t,b,5*7);jian(ans,t);    div(t,a,3);jian(ans,t);    div(t,a,5);jian(ans,t);    div(t,a,7);jian(ans,t);    div(t,a,3*5*7);jian(ans,t);    for(int i=ans[0];i>=1;i--)    cout<<ans[i];    fclose(stdin);fclose(stdout);    return 0;}