#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define X first #define Y second #define pb push_back #define rep(X,Y) for (int (X) = 0;(X) < (Y);++(X)) #define rrep(X,Y) for (int (X) = (Y-1);(X) >=0;--(X)) #define repe(X,Y) for ((X) = 0;(X) < (Y);++(X)) #define peat(X,Y) for (;(X) < (Y);++(X)) #define all(X) (X).begin(),(X).end() #define rall(X) (X).rbegin(),(X).rend() using namespace std; typedef long long ll; typedef pair pii; typedef pair pll; template using vv=vector>; template ostream& operator<<(ostream &os, const vector &t) { os<<"{"; rep(i,t.size()) {os< ostream& operator<<(ostream &os, const pair &t) { return os<<"("<1)n=(n%(m-1)+m-1)%(m-1); if(n==0)return 1ll; ll tmp=modpow(r,n/2,m); return tmp*tmp%m*(n%2?r:1ll)%m; } const int MAX_P=2123; int prime[MAX_P]; map Fact(int n){ map re; while(prime[n]){ ++re[prime[n]]; n/=prime[n]; } if(n)++re[n]; return re; } void pre(){ int i,j,k; for(i=2;i fact=Fact(n); int re=n; for(auto pq:fact){ re/=pq.X; re*=(pq.X-1); } return re; } ll fun(ll a,ll n,ll m){ //cout<>a>>n>>m; //if(m==1){cout<<0<