#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int mod=200003,G=5,N=8e5+50;
const double pi=acos(-1);
int n,id[N],r[N],lim=1,val[N];LL ans;
struct node{
    double x,y;
    node friend operator +(node a,node b){return node{a.x+b.x,a.y+b.y};}
    node friend operator -(node a,node b){return node{a.x-b.x,a.y-b.y};}
    node friend operator *(node a,node b){return node{a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x};}
}c[N];
void FFT(node *a,int n,int op){
    for(int i=0;i<n;i++)r[i]=r[i>>1]>>1^(i&1?n>>1:0);
    for(int i=0;i<n;i++)if(i<r[i])swap(a[i],a[r[i]]);
    for(int i=1;i<n;i<<=1){
        node Wn=node{cos(pi/i),sin(pi/i)};
        for(int j=0;j<n;j+=i<<1){
            node w=node{1,0};
            for(int k=0;k<i;k++,w=w*Wn){
                node x=a[j+k],y=a[i+j+k]*w;
                a[j+k]=x+y;a[i+j+k]=x-y;
            }
        }
    }
    if(op==-1){
        reverse(a+1,a+n);
        for(int i=0;i<n;i++)a[i].x/=n;
    }
}
int power(int x,int y){
    int z=1;
    for(;y;y>>=1,x=1ll*x*x%mod)if(y&1)z=1ll*z*x%mod;
    return z;
}
int main(){
    for(int i=1,j=0;j<mod-1;i=1ll*i*G%mod,j++)id[i]=j,val[j]=i;
    scanf("%d",&n);
    for(int i=1,x;i<=n;i++){
        scanf("%d",&x);
        ans-=1ll*x*x%mod;
        if(x)c[id[x]].x++;
    }
    while(lim<mod*2)lim<<=1;
    FFT(c,lim,1);
    for(int i=0;i<lim;i++)c[i]=c[i]*c[i];
    FFT(c,lim,-1);
    for(int i=0;i<mod*2;i++)
        ans+=(LL)(c[i].x+0.5)*val[i%(mod-1)];
    cout<<ans/2;
    return 0;
}
### Prompt

In Cpp, your task is to solve the following problem:

Let’s take a prime P = 200\,003. You are given N integers A_1, A_2, \ldots, A_N. Find the sum of ((A_i \cdot A_j) \bmod P) over all N \cdot (N-1) / 2 unordered pairs of elements (i < j).

Please note that the sum isn't computed modulo P.

Constraints

* 2 \leq N \leq 200\,000
* 0 \leq A_i < P = 200\,003
* All values in input are integers.

Input

Input is given from Standard Input in the following format.


N
A_1 A_2 \cdots A_N


Output

Print one integer — the sum over ((A_i \cdot A_j) \bmod P).

Examples

Input

4
2019 0 2020 200002


Output

474287


Input

5
1 1 2 2 100000


Output

600013

### Response

```cpp
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int mod=200003,G=5,N=8e5+50;
const double pi=acos(-1);
int n,id[N],r[N],lim=1,val[N];LL ans;
struct node{
    double x,y;
    node friend operator +(node a,node b){return node{a.x+b.x,a.y+b.y};}
    node friend operator -(node a,node b){return node{a.x-b.x,a.y-b.y};}
    node friend operator *(node a,node b){return node{a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x};}
}c[N];
void FFT(node *a,int n,int op){
    for(int i=0;i<n;i++)r[i]=r[i>>1]>>1^(i&1?n>>1:0);
    for(int i=0;i<n;i++)if(i<r[i])swap(a[i],a[r[i]]);
    for(int i=1;i<n;i<<=1){
        node Wn=node{cos(pi/i),sin(pi/i)};
        for(int j=0;j<n;j+=i<<1){
            node w=node{1,0};
            for(int k=0;k<i;k++,w=w*Wn){
                node x=a[j+k],y=a[i+j+k]*w;
                a[j+k]=x+y;a[i+j+k]=x-y;
            }
        }
    }
    if(op==-1){
        reverse(a+1,a+n);
        for(int i=0;i<n;i++)a[i].x/=n;
    }
}
int power(int x,int y){
    int z=1;
    for(;y;y>>=1,x=1ll*x*x%mod)if(y&1)z=1ll*z*x%mod;
    return z;
}
int main(){
    for(int i=1,j=0;j<mod-1;i=1ll*i*G%mod,j++)id[i]=j,val[j]=i;
    scanf("%d",&n);
    for(int i=1,x;i<=n;i++){
        scanf("%d",&x);
        ans-=1ll*x*x%mod;
        if(x)c[id[x]].x++;
    }
    while(lim<mod*2)lim<<=1;
    FFT(c,lim,1);
    for(int i=0;i<lim;i++)c[i]=c[i]*c[i];
    FFT(c,lim,-1);
    for(int i=0;i<mod*2;i++)
        ans+=(LL)(c[i].x+0.5)*val[i%(mod-1)];
    cout<<ans/2;
    return 0;
}
```