Variational operator learning: A unified paradigm marrying training neural operators and solving partial differential equations