Skip to Content

3. Sieve methods

2019-09-04·Combinatorics

PIE

  • PIE(Priciple of Inclusion-Exclusion)
    • AI=⋂i∈IAiA_I=\bigcap_{i\in I}A_i
    • ∣⋃i=1nAi∣=∑I⊆{1,⋯ ,n},∣I∣≥1(−1)∣I∣−1∣AI∣|\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)∣I∣∣AI∣|\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:N→Mf:N\rightarrow M
  • U={f:[n]→[m]},∣U∣=mnU=\{f:[n]\rightarrow[m]\}, |U|=m^n
  • Ai={f:[n]→[m]\{i}},∣Ai∣=(m−1)nA_i=\{f:[n]\rightarrow[m]\backslash\{i\}\},|A_i|=(m-1)^n
  • AI={f:[n]→[m]\I},∣AI∣=(m−∣I∣)nA_I=\{f:[n]\rightarrow[m]\backslash I\},|A_I|=(m-|I|)^n
  • ∣⋂i=1nAi‾∣=∑I⊆{1,⋯ ,n}(−1)∣I∣∣AI∣=∑j=0m(−1)j(mj)(m−j)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∣=(n−1)!|A_i|=(n-1)!
  • ∣AI∣=(n−∣I∣)!|A_I|=(n-|I|)!
  • ∑I⊆[n](−1)∣I∣∣AI∣=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=∣{π∣B∩Gπ=∅}N_0=|\{\pi|B\cap G_\pi=\emptyset\} (Non-attacking rooks)
  • rk:∣{S∈(Bk)∣∀(i1,j1),(i2,j2)∈S,i1≠i2,j1≠j2}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(n−k)!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),⋯ ,(n−1,n−1),(0,1),(1,2),⋯ ,(n−2,n−1),(n−1,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: (m−k+1k){m-k+1\choose k}
      • circle: mm−1(m−kk)\frac{m}{m-1}{m-k\choose k}
  • N0=∑k=0n(−1)k2n2n−k(2n−kk)(n−k)!N_0=\sum_{k=0}^n(-1)^k\frac{2n}{2n-k}{2n-k\choose k}(n-k)!

Euler totient function

  • ϕ(n)=n∏i=1r(1−1pi)\phi(n)=n\prod_{i=1}^r(1-\frac{1}{p_i})
  • ∑d∣nϕ(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: x≤yx\leq y or y≤xy\leq x
  • incidence algebra: algebra I(P)I(P)
    • I(P)={α:P×P→R∣α(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)=∑z∈Pα(x,z)β(z,y)=∑x≤z≤yα(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)=[x≤Py]∈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]
    • ∑x≤z≤yμ(x,z)=[x=y]\sum_{x\leq z\leq y}\mu(x,z)=[x=y]
μ(x,y)={−∑x≤z<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 x≤x′,y≤y′x\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)=∑y∈Pf(y)ζ(y,x)=∑y≤xf(y)(f\zeta)(x)=\sum_{y\in P}f(y)\zeta(y,x)=\sum_{y\leq x}f(y)
  • (gμ)(x)=∑y∈Pg(y)μ(y,x)=∑y≤xg(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:P→Rf,g:P\rightarrow \mathbb{R}
  • Mobius inversion formula
    • ∀x∈P,g(x)=∑y≤xf(y)  ⟺  ∀x∈P,f(x)=∑y≤xg(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=j−1i+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)∣T∣−∣S∣S⊆T0o.w.\mu(S,T)= \begin{cases} (-1)^{|T|-|S|} & S\subseteq T\newline 0 & o.w. \end{cases}
  • g(J)=∣⋂i∈JAi∣=∑i⊇Jf(I)g(J)=|\bigcap_{i\in J}A_i|=\sum_{i\supseteq J}f(I)
  • belongs to exactly the sets Ai,i∈JA_i,i\in J: f(J)=∣(⋂i∈JAi)\(⋃i∉JAi)∣=∑I⊇Jg(I)μ(J,I)=∑I⊇J(−1)∣I∣−∣J∣∣⋂i∈IAi∣f(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)∣I∣∣⋂i∈IAi∣f(\emptyset)=|U\backslash\bigcup_iA_i|=\sum_{I\subseteq[n]}(-1)^{|I|}|\bigcap_{i\in I}A_i|
    • ∑T⊆I⊆S(−1)∣S∣−∣I∣=0\sum_{T\subseteq I\subseteq S}(-1)^{|S|-|I|}=0 if S≠TS\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) = ∑π∈Sn∏i∈[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: ∑π∈Sn∏i∈[n]Ai,π(i)=∑I⊆[n](−1)n−∣I∣∏i∈[n]∑j∈IAi,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>0∣a∣n},≤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)
    • ∑d∣nμ(d)=[n=1]\sum_{d|n}\mu(d)=[n=1]
    • ϕζ=id\phi\zeta=\text{id}
    • ϕ=idμ=∑d∣nμ(d)nd\phi=\text{id}\mu=\sum_{d|n}\mu(d)\frac{n}{d}
  • g(n)=∑d∣nf(d)g(n)=\sum_{d|n}f(d)
  • f(n)=∑d∣ng(d)μ(nd)f(n)=\sum_{d|n}g(d)\mu(\frac{n}{d})