知识点 - 容斥原理 鸽巢原理 乘法原理

知识点 - 容斥原理 鸽巢原理 乘法原理

解决问题类型：

容斥：算重了。

1. 首位大于1，末尾小于8的排列数
2. 长度为n的由 (0, 1, 2) 构成，且0，1，2至少出现一次的子序列个数
3. n元一次方程整数解个数。
4. 给定区间互质数对个数
5. 给n个，问给定区间中是某个$a_i$的倍数的个数
6. 符合给定规则的串的个数
7. 从一个格子走到另一个格子的方案数
8. 互质四元组的个数
9. harmonic triplets的个数
10. 错排问题

鸽巢：

题目给了n，考虑n+1的情况。

结论

容斥公式

（1）设$A_1,A_2,\cdot \cdot \cdot ,A_n$是集合 $S$ 的子集，表示以集合 $S$ 代表可能发生的事件中的 $n$ 个子事件，$\left|A_i\right|$表示子事件$A_i$发生的个数$\left(0\le i\le n\right)$，则有
$$\left|A_1\cup A_2\cup A_3\cup \cdot \cdot \cdot \cup A_n\right|=\sum _{i=1}^n\left|A_i\right|-\sum _{1\le i<j\le n}\left|A_i\cap A_j\right|+\cdot \cdot \cdot +{\left(-1\right)}^{n-1}\left|A_1\cap A_2\cap A_3\cap \cdot \cdot \cdot \cap A_n\right|$$

压缩版：

$$\begin{gathered} \left|\bigcup _{i=1}^n{A}_i\right|=\sum _{\varnothing ̸=J\subseteq \\ 1,2,\dots ,n \\ }(-1{)}^{|J|-1}|\bigcap _{j\in J}{A}_j| \end{gathered}$$

概率版：

KaTeX parse error: Expected '}', got '\cal' at position 2: {\̲c̲a̲l̲ ̲P}(A)代表A事件发生的概率。
KaTeX parse error: Expected '}', got '\cal' at position 3: {\̲c̲a̲l̲ ̲P} \left(\bigcu…

补集公式：

（A上面有横线）
$$\left|\overline{A_1}\cap \overline{A_2}\cap \overline{A_3}\cap \cdot \cdot \cdot \cap \overline{A_n}\right|=S-\left|A_1\cup A_2\cup A_3\cup \cdot \cdot \cdot \cup A_n\right|$$

补集版

$$\left|\bigcap _{i=1}^n\overline{A_i}\right|=\sum _{m=0}^n(-1{)}^m\sum _{|X|=m}\left|\bigcap _{i\in X}{A}_i\right|$$

集合大小为r版

$$\left|\bigcup _{|B|=r}\left[\bigcap _{i\in B}{A}_i\cap \bigcap _{j̸\in B}\overline{A_j}\right]\right|=\sum _{m=r}^n(-1{)}^{m-r}(\genfrac{}{}{0ex}{}{m}{r})\sum _{|X|=m}\left|\bigcap _{i\in X}{A}_i\right|$$

容斥问题

上面所有公式证明和问题题解戳这里

错排问题

（1）设$D_n$表示$1,2,3,\cdot \cdot \cdot ,n$这n个数的一个排列的错排个数，有
$$D_n=n!\left[1-\frac{1}{1!}+\frac{1}{2!}-\frac{1}{3!}+\cdot \cdot \cdot +{\left(-1\right)}^n\frac{1}{n!}\right],D_1=0,D_2=1$$

$$D_n=(n-1)[D_{n-1}+D_{n-2}],n>2$$

$$D_n=n{D}_{n-1}+(-1{)}^n$$

带有禁位的错排问题

（1）n个元素$a_1,a_2,a_3,\cdot \cdot \cdot ,a_n$带有禁位$X_1,X_2,X_3,\cdot \cdot \cdot ,X_n$的错排数为
$$P\left(X_1,X_2,X_3,\cdot \cdot \cdot ,X_n\right)=n!-r_1\left(n-1\right)!+r_2\left(n-2\right)!-\cdot \cdot \cdot +{\left(-1\right)}^k{r}_k\left(n-k\right)!+\cdot \cdot \cdot +{\left(-1\right)}^n{r}_n$$
式中$r_k$表示有 $k$ 个元素在禁位上的个数

复杂度：

例题

HDU2204 容斥：给你一个正整数N，确定在1到N之间有多少个可以表示成M^K（K>1)的数。

分析：小于$N$的$M^k$的个数为$\lfloor \sqrt[k]{N}\rfloor$。枚举所有的k的话。显然有东西算重了。(比如1,我们可以每个k都不统计，最后加1)

首先，如果k不是质数$e.g.M^k=M^{p1p2}$，则$M^k$就以$(M^{p1}{)}^{p2}$的形式被$\lfloor \sqrt[p_2]{N}\rfloor$统计过了。于是所有和数都不用统计。

其次，我们发现上面的例子里$M^k$被$\lfloor \sqrt[p_1]{N}\rfloor$$\lfloor \sqrt[p_2]{N}\rfloor$都统计了。于是我们做容斥，K包含1个质数-K包含2个质数+3个…

最后，注意到$2^{60}>10^{18}$所以容斥到三个质数即可，最大的质数取到59。

hdu1205鸽巢：显然的想法题，但不会用鸽巢原理证明。

hdu5525 约数乘积: $n=p_1^{e_1}\cdot p_2^{e_2}\cdots p_k^{e_k}$对于每个质数，它被乘入答案的次数为${\Sigma }_{i=1}^{i\le e}i\cdot d(\frac{n}{p^e})$,d(n)为因数个数。
$$\prod _{i=1}^{i\le k}{p}_i^{({\Sigma }_{j=1}^{j\le e_i}j\cdot d(\frac{n}{p_i^{e_i}}))}$$
注意指数上要模p-1，所以不能用除法。对于/2有$2|(n+1)n$ 直接把2乘到模数上。对于$(1+e_i)$维护前后缀积。

代码

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
int prime[20]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59};
long long ans,n;
int i;
void dfs(int j,int num,int p){
    if(p==0){
        long long t=pow(n,1.0/num);
        if(pow(t,0.0+num)>n)  t--;
        t--;
        if(t>0)
            //i为奇数则加，否则减
            ans+=t*(i&1?1:-1);
        return ;
    }
    if(j>=17)
        return;//超出60以内素数，不予考虑
    if(num*prime[j]<60) //仍在范围内，继续搜
        dfs(j+1,num*prime[j],p-1);
    dfs(j+1,num,p);
}
int main(){
    while(scanf("%I64d",&n)!=EOF){
        ans=0;
        for(i=1;i<=3;i++)
            dfs(0,1,i);
        printf("%I64d\n",ans+1);
    }
    return 0;
 
}

//T2
const int N=2e5+10,M=4e6+10,inf=2147483647;
const ll INF=1e18+10,mod=1e9+7;

///   数组大小
int prime(int n)
{
    if(n<=1)
    return 0;
    if(n==2)
    return 1;
    if(n%2==0)
    return 0;
    int k, upperBound=n/2;
    for(k=3; k<=upperBound; k+=2)
    {
        upperBound=n/k;
        if(n%k==0)
            return 0;
    }
    return 1;
}

vector<int>p;
int si[N];
void init()
{
    for(int i=2;i<=100000;i++)
        if(prime(i))
        p.push_back(i);
}
ll quick(ll a,ll b,ll c)
{
    ll ans=1;
    while(b)
    {
        if(b&1)ans=(ans*a)%c;
        b>>=1;
        a=(a*a)%c;
    }
    return ans;
}
int main()
{
    init();
    int n;
    while(~scanf("%d",&n))
    {
        memset(si,0,sizeof(si));
        for(int i=1;i<=n;i++)
        {
            int z;
            int x=i;
            scanf("%d",&z);
            for(int j=0;j<p.size();j++)
            {
                if(1LL*p[j]*p[j]>x)break;
                if(x==1)break;
                while(x%p[j]==0)
                {
                    si[j]+=z;
                    si[j]=(si[j])%(2*(mod-1));
                    x/=p[j];
                }
            }
            if(x!=1)
            {
                int pos=lower_bound(p.begin(),p.end(),x)-p.begin();
                si[pos]+=z;
            }
        }
        ll sum=1;
        for(int i=0;i<p.size();i++)
        {
            sum*=(si[i]+1);
            sum%=(2*(mod-1));
        }
        //cout<<sum<<endl;
        ll ans=1;
        for(int i=0;i<p.size();i++)
        {
            ans=(ans*quick(p[i],((sum*si[i])%(2*(mod-1)))/2,mod))%(mod);
        }
        printf("%lld\n",ans);
    }
    return 0;
}