Skip to Content
Combinatorics

3. Sieve methods

2019-09-04Original-language archivelegacy assets may be incomplete

PIE

  • PIE(Priciple of Inclusion-Exclusion)
    • AI=iIAiA_I=\bigcap_{i\in I}A_i
    • i=1nAi=I{1,,n},I1(1)I1AI|\bigcup_{i=1}^nA_i|=\sum_{I\subseteq\{1,\cdots,n\},|I|\geq 1}(-1)^{|I|-1}|A_I|
    • i=1nAi=I{1,,n}(1)IAI|\bigcap_{i=1}^n\overline{A_i}|=\sum_{I\subseteq\{1,\cdots,n\}}(-1)^{|I|}|A_I|
  • Sieve methods: 计算恰好具备某种性质的 \Rightarrow 计算不具有一系列坏性质的
    • Define U|U|
    • Define 'bad' properties AiA_i
    • Use PIE count the number

Surjections

  • Count Surjections f:NMf:N\rightarrow M
  • U={f:[n][m]},U=mnU=\{f:[n]\rightarrow[m]\}, |U|=m^n
  • Ai={f:[n][m]\{i}},Ai=(m1)nA_i=\{f:[n]\rightarrow[m]\backslash\{i\}\},|A_i|=(m-1)^n
  • AI={f:[n][m]\I},AI=(mI)nA_I=\{f:[n]\rightarrow[m]\backslash I\},|A_I|=(m-|I|)^n
  • i=1nAi=I{1,,n}(1)IAI=j=0m(1)j(mj)(mj)n=m!{nm}|\bigcap_{i=1}^n\overline{A_i}|=\sum_{I\subseteq\{1,\cdots,n\}}(-1)^{|I|}|A_I|=\sum_{j=0}^m(-1)^j{m\choose j}(m-j)^n=m!{n\brace m}

Derangements

  • fixed point ii: i{1,2,,n},π(i)=ii\in\{1,2,\cdots,n\},\pi(i)=i
  • derangement of {1,2,,n}\{1,2,\cdots,n\} is π\pi that has no fixed points
  • U=Sn,U=n!U=S_n,|U|=n!
  • Ai={ππ(i)=i}A_i=\{\pi|\pi(i)=i\}, Ai=(n1)!|A_i|=(n-1)!
  • AI=(nI)!|A_I|=(n-|I|)!
  • I[n](1)IAI=n!k=0n(1)kk!n!e\sum_{I\subseteq[n]}(-1)^{|I|}|A_I|=n!\sum_{k=0}^n\frac{(-1)^k}{k!}\approx\frac{n!}{e}

Permutations with restricted positions

  • B[n]×[n]B\subseteq [n]\times[n]
  • Gπ(V,E)={(i,π(i))i[n]}G_\pi(V,E)=\{(i,\pi(i))|i\in [n]\}
  • N0={πBGπ=}N_0=|\{\pi|B\cap G_\pi=\emptyset\} (Non-attacking rooks)
  • rk:{S(Bk)(i1,j1),(i2,j2)S,i1i2,j1j2}r_k:|\{S\in{B\choose k}|\forall(i_1,j_1),(i_2,j_2)\in S,i_1\not= i_2,j_1\not=j_2\}
  • N0=k=0n(1)krk(nk)!N_0=\sum_{k=0}^n(-1)^kr_k(n-k)!

Derangement problem

B={(i,i)i[n]}B=\{(i,i)|i\in [n]\}

Probleme des menages

  • nn couples at a circular table, M&F in alternate and no one sits next to his/her spouse
  • Lady: 2(n!)2(n!)
  • Men: B={(0,0),(1,1),,(n1,n1),(0,1),(1,2),,(n2,n1),(n1,0)}B=\{(0,0),(1,1),\cdots,(n-1,n-1),(0,1),(1,2),\cdots,(n-2,n-1),(n-1,0)\}
    • rkr_k: choosingn kk points, no two consecutive, from a collection of 2n2n points arranged in a circle
      • line: (mk+1k){m-k+1\choose k}
      • circle: mm1(mkk)\frac{m}{m-1}{m-k\choose k}
  • N0=k=0n(1)k2n2nk(2nkk)(nk)!N_0=\sum_{k=0}^n(-1)^k\frac{2n}{2n-k}{2n-k\choose k}(n-k)!

Euler totient function

  • ϕ(n)=ni=1r(11pi)\phi(n)=n\prod_{i=1}^r(1-\frac{1}{p_i})
  • dnϕ(d)=n\sum_{d|n}\phi(d)=n

Inversion

  • Poset PP: partially ordered set
    • Set PP with P\leq_P satisfying reflexivity, antisymmetry and transitivity
    • comparable: xyx\leq y or yxy\leq x
  • incidence algebra: algebra I(P)I(P)
    • I(P)={α:P×PRα(x,y)=0,x̸Py}I(P)=\{\alpha:P\times P\rightarrow\mathbb{R}|\alpha(x,y)=0,\forall x\not\leq_P y\} (α\alpha can be seen as matrix)
    • if α,βI(P)\alpha,\beta\in I(P) then α+βI(P)\alpha+\beta\in I(P)
    • if αI(P)\alpha\in I(P) then cαI(P)c\alpha\in I(P)
    • (αβ)(x,y)=zPα(x,z)β(z,y)=xzyα(x,z)β(z,y)(\alpha\beta)(x,y)=\sum_{z\in P}\alpha(x,z)\beta(z,y)=\sum_{x\leq z\leq y}\alpha(x,z)\beta(z,y)
  • zeta function ζ(x,y)=[xPy]I(P)\zeta(x,y)=[x\leq_P y]\in I(P)
  • Möbius function: μζ=I\mu\zeta=I
    • I(x,y)=[x=y]I(x,y)=[x=y]
    • xzyμ(x,z)=[x=y]\sum_{x\leq z\leq y}\mu(x,z)=[x=y]
μ(x,y)={xz<yμ(x,z)x<y1x=y0x≰y\mu(x,y)= \begin{cases} -\sum_{x\leq z<y}\mu(x,z) & x<y \newline 1 & x=y \newline 0 & x\not\leq y \end{cases}
  • P×QP\times Q: (x,y),(x,y)P×Q,(x,y)(x,y)\forall (x,y),(x',y')\in P\times Q,(x,y)\leq(x',y') iff xx,yyx\leq x',y\leq y'
    • product rule: μP×Q((x,y),(x,y))=μP(x,x)μQ(y,y)\mu_{P\times Q}((x,y),(x',y'))=\mu_P(x,x')\mu_Q(y,y')
  • (fζ)(x)=yPf(y)ζ(y,x)=yxf(y)(f\zeta)(x)=\sum_{y\in P}f(y)\zeta(y,x)=\sum_{y\leq x}f(y)
  • (gμ)(x)=yPg(y)μ(y,x)=yxg(y)μ(y,x)(g\mu)(x)=\sum_{y\in P}g(y)\mu(y,x)=\sum_{y\leq x}g(y)\mu(y,x)
    • f,g:PRf,g:P\rightarrow \mathbb{R}
  • Mobius inversion formula
    • xP,g(x)=yxf(y)    xP,f(x)=yxg(y)μ(y,x)\forall x\in P,g(x)=\sum_{y\leq x}f(y)\iff\forall x\in P,f(x)=\sum_{y\leq x}g(y)\mu(y,x)
    • fζ=g    f=gμf\zeta=g\iff f=g\mu
    • dual form: ζf    f=μg\zeta f\iff f=\mu g

剩余系

P=[n],PP=[n],\leq_P if i=j+1i=j+1

μ(i,j)={1i=j1i+1=j0o.w.\mu(i,j)= \begin{cases} 1 & i=j \newline -1 & i+1=j \newline 0 & o.w. \end{cases}

Boolean algebra

P=2U,PP=2^U,\leq_P is \subseteq

μ(S,T)={(1)TSST0o.w.\mu(S,T)= \begin{cases} (-1)^{|T|-|S|} & S\subseteq T\newline 0 & o.w. \end{cases}
  • g(J)=iJAi=iJf(I)g(J)=|\bigcap_{i\in J}A_i|=\sum_{i\supseteq J}f(I)
  • belongs to exactly the sets Ai,iJA_i,i\in J: f(J)=(iJAi)\(i∉JAi)=IJg(I)μ(J,I)=IJ(1)IJiIAif(J)=|(\bigcap_{i\in J}A_i)\backslash(\bigcup_{i\not\in J}A_i)|=\sum_{I\supseteq J}g(I)\mu(J,I)=\sum_{I\supseteq J}(-1)^{|I|-|J|}|\bigcap_{i\in I}A_i|
  • PIE:
    • f()=U\iAi=I[n](1)IiIAif(\emptyset)=|U\backslash\bigcup_iA_i|=\sum_{I\subseteq[n]}(-1)^{|I|}|\bigcap_{i\in I}A_i|
    • TIS(1)SI=0\sum_{T\subseteq I\subseteq S}(-1)^{|S|-|I|}=0 if STS\not=T else 11, A1=A2==An={1}A_1=A_2=\cdots=A_n=\{1\}
  • Bipartite Perfect Matching
    • Ai,j=1A_{i,j}=1 if (i,j)E(i,j)\in E
    • perm(A) = πSni[n]Ai,π(i)\sum_{\pi\in S_n}\prod_{i\in[n]}A_{i,\pi(i)} (#P-hard)
    • det(A) = πSn(1)r(π)i[n]Ai,π(i)\sum_{\pi\in S_n}(-1)^{r(\pi)}\prod_{i\in[n]}A_{i,\pi(i)} (poly-time)
    • Ryser's formula: πSni[n]Ai,π(i)=I[n](1)nIi[n]jIAi,j\sum_{\pi\in S_n}\prod_{i\in[n]}A_{i,\pi(i)}=\sum_{I\subseteq[n]}(-1)^{n-|I|}\prod_{i\in[n]}\sum_{j\in I}A_{i,j}
      • O(n!)O(n2n)O(n!)\rightarrow O(n2^n)

Number Theory

P={a>0an},pP=\{a>0|a|n\},\leq_p is |

μ(a,b)={(1)rba is the product of r distinct primes0o.w.\mu(a,b)= \begin{cases} (-1)^r & \frac{b}{a} \text{ is the product of }r\text{ distinct primes} \newline \newline 0 & o.w. \end{cases}
  • number-theoretical Mobius function: μ(nd)=μ(n,d)\mu(\frac{n}{d})=\mu(n,d)
    • dnμ(d)=[n=1]\sum_{d|n}\mu(d)=[n=1]
    • ϕζ=id\phi\zeta=\text{id}
    • ϕ=idμ=dnμ(d)nd\phi=\text{id}\mu=\sum_{d|n}\mu(d)\frac{n}{d}
  • g(n)=dnf(d)g(n)=\sum_{d|n}f(d)
  • f(n)=dng(d)μ(nd)f(n)=\sum_{d|n}g(d)\mu(\frac{n}{d})