When we want to integrate a function f over a manifold M, we may meet some problems, for example, the problem showed in the picture below:

Then people used differential form to integrate. But it confused me:does that really solve the problem of integrating the function f on M? How can f be related directly to a k-form $\omega$? Is there a k-form $\omega$ to represent a specific function f on manifold M?
Differential forms are not introduced to answer the question "How do I integrate functions on manifolds?"
Instead, they are introduced to answer a different question, namely "What are the correct objects to integrate on manifolds?"
The answer to that question is: differential forms! More precisely, the correct objects to integrate on an orientable $m$-dimensional manifold are differential $m$-forms.
The reason they are correct is precisely because the integral of a differential form does transform correctly, without its value changing, under a change of coordinates.
Added: You say you still want to know how to integrate a function on a manifold. As I have tried to explain briefly so far, this is a quixotic quest. One should not attempt to integrate functions on manifolds, because it will not be well-defined when you change coordinates.
Here is an example to better illustrate why instead we integrate differential forms, not functions.
In multivariable calculus you learn that total mass equals the integral of density. What are mass and density? Are they functions? If not, what are they?
Let's take a look at a typical mass integral such as $\int_U f \, dx$. The numerical value of this integral represents the total mass. Now, we cannot measure the mass at a point, it does not make mathematical sense. But what we can measure is the infinitesmal ratio of mass per unit volume at a point: that's called the density function, and that's the integrand $f$ in the density integral. Also, the expression $dx$ represents an infinitesmal amount of volume. The expression $f \, dx$ then represents the product of the infinitesmal ratio of mass per unit volume, multiplied by infinitesmal volume; so you can think of $f \, dx$ as a measurement of a total mass that has been spread around over the whole of the set $U$.
That's one way to think of a differential form such as $f \, dx$ (at least in the special case of a differential $k$-form on a $k$-dimensional manifold). Namely, it is a physical description of a total mass that has been spread over the whole of the manifold.
Note well: the differential form is not the density function $f$; the actual numerical value of density per unit volume depends on the units you choose to measure volume. Instead, the differential form is the density function multiplied by the volume form, namely $f \, dx$. This quantity is well-defined independent of the coordinates you choose.
Let's look at an example to illustrate about what happens under a coordinate change. Start with imagining $U$ as a box $\mathbb{R}^3$ with coordinate functions $x_1,x_2,x_3$, of side length $2$ and volume $8$. Imagine also that we have a total mass $M$ of constant density, the value of that density being $M/8$ at every point. The total mass is $\int_U M/8 \, dx = M/8 \int_U dx = M/8 \times \text{(volume of $U$ in $x$-coordinates)} = M/8 \times 8 = M$.
What happens when you change coordinates? Well, suppose you change coordinates by a function of the form $$(y_1,y_2,y_3) = F(x_1,x_2,....,x_k) = (\frac{3}{2} \, x_1,\, \frac{3}{2} \, x_2,\, \frac{3}{2} \, x_3) $$ This corresponds to getting a new ruler, the $y$-ruler, whose unit length is $2/3$ the size of the unit length on the $x$-ruler you started with. Under this coordinate change function $F$, the box with side length $2$, volume $2^k$, and mass $M$ in $x$-coordinates becomes, in $y$-coordinates, a box with side length $3$, volume $3^k$, and the mass $M$ is unchanged.
Mass does not change, just because you are using a different ruler.
In the new $y$-coordinates, the density function is $M/3^3=M/27$ per unit volume. And the new mass integral is $\int_U M/27 \, dy = M/27 \times \text{(volume of $U$ in $y$-coordinates)} = M/27 \times 27 = M$.
To summarize this example, the expressions $M/8 \, dx$ and $M/27 \, dy$ represent the same differential form on $U$, expressed in two different coordinate systems. The integral of this differential form over $U$ has a well-defined value, namely $M$, the total mass.