Let's derive a generalization of the Inclusion-Exclusion Principle based on some combinatoric identities.
Cancellation Lemma
Suppose $n$ and $k$ are non-negative integers. Then
$$
\sum_{j=k}^n(-1)^{j-k}\binom{n}{j}\binom{j}{k}
=\left\{\begin{array}{}
1&\text{if }n=k\\
0&\text{otherwise}
\end{array}\right.\tag{1}
$$
Note that if $n\lt k$, then $\binom{n}{j}\binom{j}{k}=0$ for all $j$, so we can assume that $n\ge k$.
Proof 1: Using $\displaystyle\binom{n}{k}=\frac{n!}{k!\,(n-k)!}$, we get
$$
\begin{align}
\sum_{j=k}^n(-1)^{j-k}\binom{n}{j}\binom{j}{k}
&=\sum_{j=k}^n(-1)^{j-k}\binom{n}{k}\binom{n-k}{j-k}\\
&=\binom{n}{k}(1-1)^{n-k}\tag{2}
\end{align}
$$
which is $1$ if $n=k$ and $0$ otherwise. $\quad\square$
Proof 2: Using Vandermonde's Identity, we get
$$
\begin{align}
\sum_{j=k}^n(-1)^{j-k}\binom{n}{j}\binom{j}{k}
&=\sum_{j=k}^n(-1)^{j-k}\binom{n}{n-j}\binom{j}{j-k}\\
&=\sum_{j=k}^n\binom{n}{n-j}\binom{-k-1}{j-k}\\
&=\binom{n-k-1}{n-k}\tag{3}
\end{align}
$$
which is $1$ if $n=k$ and $0$ otherwise. $\quad\square$
Proof 3: Consider the generating function for $(1)$ as a function of $k$:
$$
\begin{align}
\sum_{k=0}^n\sum_{j=k}^n(-1)^{j-k}\binom{n}{j}\binom{j}{k}x^k
&=\sum_{j=0}^n\sum_{k=0}^j(-1)^{j-k}\binom{n}{j}\binom{j}{k}x^k\\
&=\sum_{j=0}^n(-1)^j\binom{n}{j}(1-x)^j\\
&=(1-(1-x))^n\\[12pt]
&=x^n\tag{4}
\end{align}
$$
Equating the coefficients of $x^k$ in $(4)$ proves the result. $\quad\square$
Theorem (Generalized Inclusion-Exclusion Principle)
Let $\{S(i)\}_{i=1}^m$ be a finite collection of sets from a finite universe.
Let $N(j)$ be the sum of the sizes of all intersections of $j$ of the
$S(i)$:
$$
N(j)=\sum_{|A|=j}\left|\,\bigcap_{i\in A} S(i)\,\right|\tag{5}
$$
Thus, $N(0)$ is the size of the universe.
Then, the number of elements in exactly $k$ of the $S(i)$ is
$$
\sum_{j=0}^m(-1)^{j-k}\binom{j}{k}N(j)\tag{6}
$$
Proof:
An element that is in $n$ of the $S(i)$ is counted $\binom{n}{j}$ times in $N(j)$.
That is, there are $\binom{n}{j}$ ways to choose the $j$ sets in the intersection
from the $n$ sets that contain the element. The Cancellation Lemma says
that this element is counted once in $(6)$ if $n = k$ and is cancelled out
otherwise. $\quad\square$
Using the identity
$$
\sum_{k=0}^n(-1)^k\binom{j}{k}=(-1)^n\binom{j-1}{n}\tag{7}
$$
where $\binom{-1}{n}=$$(-1)^n\binom{n}{n}$$=(-1)^n$$[n\ge0]$, we get the following
Corollary 1:
The number of elements in at most $k$ of the $S(i)$ is
$$
\sum_{j=0}^m(-1)^{j-k}\binom{j-1}{k}N(j)\tag{8}
$$
Using the identity
$$
\sum_{k=n}^j(-1)^k\binom{j}{k}=(-1)^n\binom{j-1}{j-n}\tag{9}
$$
again where $\binom{-1}{n}=(-1)^n\binom{n}{n}=(-1)^n[n\ge0]$, we get the following
Corollary 2:
The number of elements in at least $k$ of the $S(i)$ is
$$
\sum_{j=k}^m(-1)^{j-k}\binom{j-1}{j-k}N(j)\tag{10}
$$
The Inclusion-Exclusion Principle
The usual Inclusion-Exclusion Principle is the case $k=1$ of Corollary 2.
However, the usual Inclusion-Exclusion Principle can also be derived by subtracting the number of objects that are in none of the $S(i)$ from the size of the universe. Using $(6)$ yields
$$
\underbrace{\ \ \ N(0)\ \ \ \vphantom{\sum_j\binom{j}{0}}}_{\substack{\text{size of the}\\\text{universe}}}-\underbrace{\sum_{j=0}^m(-1)^{j}\binom{j}{0}N(j)}_{\substack{\text{number of elements}\\\text{in none of the $S(i)$}}}=\underbrace{\sum_{j=1}^m(-1)^{j-1}N(j)\vphantom{\binom{j}{0}}}_{\substack{\text{number of elements in}\\\text{at least one of the $S(i)$}}}
$$
Best Answer
As Asaf said in the comments, this is bookkeeping and induction. I’m going to write out the induction step in full, but I’m going to do it in a very compact fashion that may not be easy to follow at first. I’m doing this for two reasons. First, you really ought to try to work out the induction step on your own, so I don’t want to make it too easy to read. (If you have trouble doing it in general, I suggest writing out the argument for $n=3$ and perhaps $n=4$ first; at that point you should be able to see fairly clearly what’s going on, and the main problem will likely be expressing it clearly in the general case.) Secondly, at some point you’ll want to get used to reading this sort of compact argument.
First let’s write out the sum in full:
$$\left|\bigoplus_{k=1}^nA_k\right|=\sum_{k=1}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\left|\bigcap_{i\in S}A_i\right|\;,\tag{1}$$
where I’ve used $\oplus$ for symmetric difference, and $[n]=\{1,\dots,n\}$. To avoid having too many subscripts, let $\chi_k$ be the characteristic function of $A_k$, and let $\varphi_n$ be the characteristic function of $\bigoplus_{k=1}^nA_k$. Then $(1)$ will certainly be true if
$$\varphi_n=\sum_{k=1}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\prod_{i\in S}\chi_i\;,\tag{2}$$
since $\prod_{i\in S}\chi_i$ is the characteristic function of $\bigcap_{i\in S}A_i$.
You can prove $(2)$ by induction on $n$. You need to use the fact that for any characteristic function $\chi$, $\chi^2=\chi$. You already know that $\chi_{A\oplus B}=(\chi_A-\chi_B)^2$. Now for the induction step you have
$$\begin{align*} \varphi_{n+1}&=\left(\varphi_n-\chi_{n+1}\right)^2\\ &=\varphi_n^2+\chi_{n+1}^2-2\varphi_n\chi_{n+1}\\ &=\varphi_n+\chi_{n+1}-2\varphi_n\chi_{n+1}\\ &=\sum_{k=1}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\prod_{i\in S}\chi_i+\chi_{n+1}-2\chi_{n+1}\sum_{k=1}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\prod_{i\in S}\chi_i\tag{3} \end{align*}$$
for starters. The first term in $(3)$ can be split as
$$\sum_{i\in[n]}\chi_i+\sum_{k=2}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\prod_{i\in S}\chi_i\tag{4}$$
and the first term of $(4)$ can be combined with the $\chi_{n+1}$ term in $(3)$ to yield
$$\begin{align*} \varphi_{n+1}&=\sum_{i\in[n+1]}\chi_i+\sum_{k=2}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\prod_{i\in S}\chi_i-2\chi_{n+1}\sum_{k=1}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\prod_{i\in S}\chi_i\\ &=\sum_{i\in[n+1]}\chi_i+\sum_{k=2}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\prod_{i\in S}\chi_i-2\sum_{k=1}^n(-1)^{k-1}2^{k-1}\sum_{S\subseteq[n]\atop|S|=k}\prod_{i\in S}\chi_i\chi_{n+1}\tag{5}\;. \end{align*}$$
Now let $S$ be a subset of $[n+1]$ such that $|S|=k$. If $k=1$, then $\prod_{i\in S}$ is one of the $\chi_i$, and it appears in the first term of $(5)$ with coefficient $1=(-1)^02^0=(-1)^{k-1}2^{k-1}$.
Suppose now that $1<k\le n$. If $n+1\notin S$, $\prod_{i\in S}\chi_i$ is counted in the second term of $(5)$, where it appears with coefficient $(-1)^{k-1}2^{k-1}$. If $n+1\in S$, then $\prod_{i\in S}\chi_i$ appears in the last term of $(5)$ as $\prod_{i\in S'}\chi_i\chi_{n+1}$, where $S'=S\setminus\{n+1\}$, and $|S'|=k-1$. In this case the coefficient of $\prod_{i\in S}\chi_i$ is $-2(-1)^{k-2}2^{k-2}=(-1)^{k-1}2^{k-1}$.
Finally, suppose that $k=n+1$. Then $S=[n+1]$, and $\prod_{i\in S}\chi_i=\prod_{i\in[n]}\chi_i\chi_{n+1}$ appears in the last term of $(5)$ with coefficient $-2(-1)^{n-1}2^{n-1}=(-1)^n2^n=(-1)^{k-1}2^{k-1}$.
It follows that $(5)$ (and hence $\varphi_{n+1}$) is equal to
$$\sum_{k=1}^{n+1}(-1)^{k-1}2^{k-1}\sum_{S\subset[n+1]\atop|S|=k}\prod_{i\in S}\chi_i\;,$$
which completes the induction step.
By the way, a simpler approach is to show that $a\in\bigoplus_{k=1}^nA_k$ iff $|\{k\in[n]:a\in A_k\}|$ is odd, which is an easy induction, and then to show that that the righthand side of $(1)$ counts the elements of $\bigcup_{k=1}^nA_k$ that are in an odd number of the $A_k$. This is a fairly easy consequence of the binomial theorem. If an element $a$ is in exactly $m$ of the $n$ sets, for each $k$ from $1$ through $m$ it’s in $\binom{m}k$ $k$-fold intersections, so it contributes
$$\begin{align*} \sum_{k=1}^m\binom{m}k(-1)^{k-1}2^{k-1}&=-\frac12\sum_{k=1}^m\binom{m}k(-2)^k\\ &=-\frac12\left(\sum_{k=0}^m\binom{m}k(-2)^k-1\right)\\ &=\frac12\Big(1-(1-2)^m\Big)\\ &=\frac{1-(-1)^m}2\\ &=\begin{cases} 1,&\text{if }m\text{ is odd}\\ 0,&\text{if }m\text{ is even} \end{cases} \end{align*}$$
to the righthand side of $(1)$.