#include #define rep(i, a, n) for(int i = a; i < (n); i++) using namespace std; using ll = long long; using P = pair; const int INF = 1001001001; const ll LINF = 1001002003004005006ll; //const int mod = 1000000007; //const int mod = 998244353; //FFT namespace FFT{ struct num{ double x,y; num(){x=y=0;} num(double x,double y):x(x),y(y){} }; inline num operator+(num a,num b){ return num(a.x+b.x,a.y+b.y);} inline num operator-(num a,num b){ return num(a.x-b.x,a.y-b.y);} inline num operator*(num a,num b){ return num(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);} inline num conj(num a){ return num(a.x,-a.y);} int base=1; vector rts={{0,0},{1,0}}; vector rev={0,1}; const double PI=acosl(-1.0); void ensure_base(int nbase){ if(nbase<=base) return; rev.resize(1<>1]>>1)+((i&1)<<(nbase-1)); rts.resize(1< &a,int n=-1){ if(n == -1) n = (int)a.size(); assert((n&(n-1))==0); int zeros=__builtin_ctz(n); ensure_base(zeros); int shift=base-zeros; for(int i=0;i>shift)) swap(a[i],a[rev[i]>>shift]); for(int k=1;k fa; template vector multiply(const vector &a,const vector &b){ int need=(int)a.size()+b.size()-1; int nbase=0; while((1<(int)fa.size()) fa.resize(sz); for(int i=0;i>1);i++){ int j=(sz-i)&(sz-1); num z=(fa[j]*fa[j]-conj(fa[i]*fa[i]))*r; if(i!=j) fa[j]=(fa[i]*fa[i]-conj(fa[j]*fa[j]))*r; fa[i]=z; } fft(fa,sz); vector res(need); for(int i=0;i> n >> x; vector a(n); rep(i, 0, n) cin >> a[i]; vector cnt(100005); rep(i, 0, n) cnt[a[i]]++; auto c = FFT::multiply(cnt, cnt); cout << c[x] << endl; return 0; }