牛客容斥原理专题-无关(relationship)

目录

原题

 解题历程

最开始的思路

经过许久的思考之后的分析

引入

AC代码

---

原题

原题链接如下：点这里呀呀呀

下面这个是这个专题的链接：(●ˇ∀ˇ●)

牛客竞赛_ACM/NOI/CSP/CCPC/ICPC算法编程高难度练习赛_牛客竞赛OJ (nowcoder.com)

 

 

 解题历程

最开始的思路

 最开始当让想到的是，呀，这个数据开个long long就够了，这样子没什么问题，之后用数组把这个大小为k的数组存起来，然后写两层循环这样统计起来，看有多少的数是无关的，这些繁琐的事情让电脑去做，，但是！！！！！！！！L和R的数据范围是10^18,  k<=20,这样一看，呀！这个循环好大，必定超时，看，我开始就是这样写的。/(ㄒoㄒ)/~~

#include<bits/stdc++.h>
typedef long long ll;
ll l,r,k;
const int N=25;
ll a[N];
using namespace std;
int main()
{
    cin>>l>>r>>k;
    ll sum=0,ans=0;
    for(int i=0;i<k;i++)
        cin>>a[i];
    for(ll i=l;i<=r;i++)//看这里超大的两层循环，这就必定超时了/(ㄒoㄒ)/~~
      {  
          sum=0;
         for(int j=0;j<k;j++)
          {
            if(i%a[j]==0)
                sum++;
          }
            if(sum==0)
                ans++;
    }
    cout<<ans<<endl;
    return 0;
}

 看,wuwuuwuwu,太惨了

经过许久的思考之后的分析

这里说了，k小于20，而且这个大小为k的数组里面的元素都是素数，就是要把n以内的那些有素数因子的数给去掉。

就比如说，说就是要把1-n,每个数，看他是不是2，3，5，7的倍数，但是他这个数可以同时是2和5的倍数，这样计算的时候就算了两次，这里就想到容斥原理。

公式如下图。

这里就相当于是n-n/a[i]+n/a[i]/a[j]......// 0<=i<k;

引入

这里引入这个专题的第一道题

 下面是这个简单题的AC代码，，你好好看看，找找感觉，能理解这个题目的时候，你大概再看看就可以理解这个题目了

#include<bits/stdc++.h>
typedef long long ll;
using namespace std;
int main()
{
    ll n;
    ll a1,a2,a3,a4;
    ll b1,b2,b3,b4,b5,b6;
    ll c1,c2,c3,c4;
    ll d;
    while(cin>>n){
    a1=n/2,a2=n/5,a3=n/11,a4=n/13;
    b1=n/10,b2=n/22,b3=n/26,b4=n/55,b5=n/65,b6=n/143;
    c1=n/110,c2=n/130,c3=n/715,c4=n/286;
    d=n/1430;
    long long sum=n-(a2+a5+a11+a13)+(b1+b2+b3+b4+b5+b6)-(c1+c2+c3+c4)+d;
    cout<<sum<<endl;
    }
    return 0;
}

这里直接写循环有很麻烦，就想到了写dfs ,毕竟这个k的值这么小，这不就是天造地设，就是得这么写吗哈哈哈哈哈哈哈哈哈哈哈哈！！！

AC代码

相当于是这样（dfs）

 

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll l,r,k,a[25];
ll dfs(int i,ll b)
{
    if(i>k) return b;
    if(!b) return 0;
    return dfs(i+1,b)-dfs(i+1,b/a[i]);
}
int main()
{
    cin>>l>>r>>k;
    for(int i=1;i<=k;i++) 
        cin>>a[i];
    cout<<dfs(1,r)-dfs(1,l-1)<<endl;
    return 0;
}

这里就是他输入l,r  所以用dfs(1，r)-dfs(1,l-1);

这样子就可以得知道结果了

dfs里面就需要你动动小手算算了~~里面是一系列的套娃操作，x=x/a[i] 又把x带进去继续套娃。