Exact matrix computation by multiple P-adic arithmetic