2020河南省ccpc省赛太阳轰炸

#include<bits/stdc++.h>

#define IOS ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);

using namespace std;


const int maxn = 5e6+5;
const int mod = 1e9 + 7;
typedef long long ll;
int n,m,k;
ll r1,r2,r,a,h;
int pm[maxn],pni[maxn];
int jc[maxn],inv[maxn];
ll p,q,p2,q2;
ll qmi(ll a,ll k,ll mod){
	ll res = 1;
	while(k){
		if(k & 1) res = 1ll*res * a % mod;
		a = 1ll*a * a % mod;
		k >>= 1; 
	}