#include<cstdio>
#include<cmath>
using namespace std;
#define N 1059000
#define M 23
int n,mod,su=1,l[M],w[M],p[M],sp,as=1,sv;
int pw(int a,int b){int as=1;while(b){if(b&1)as=1ll*as*a%mod;a=1ll*a*a%mod;b>>=1;}return as;}
struct comp{double a,b;};
comp operator +(comp a,comp b){return (comp){a.a+b.a,a.b+b.b};}
comp operator -(comp a,comp b){return (comp){a.a-b.a,a.b-b.b};}
comp operator *(comp a,comp b){return (comp){a.a*b.a-a.b*b.b,a.a*b.b+a.b*b.a};}
comp gr[2][N*2],c1[N],c2[N],c3[N],c4[N],fft[N];
int rev[N*2];
void init(int d=20)
{
	double pi=acos(-1);
	for(int l=2;l<=1<<d;l<<=1)for(int i=0;i<l;i++)rev[i+l]=(rev[(i>>1)+l]>>1)|((i&1)*(l>>1));
	for(int t=0;t<2;t++)
	for(int l=2;l<=1<<d;l<<=1)for(int i=0;i<l>>1;i++)
	gr[t][l+i]=(comp){cos(2*pi/l*i),sin(2*pi/l*i*(t?1:-1))};
}
void dft(int s,comp *a,int t)
{
	for(int i=0;i<s;i++)fft[rev[i+s]]=a[i];
	for(int l=2;l<=s;l<<=1)
	for(int i=0;i<s;i+=l)
	for(int j=0;j<l>>1;j++)
	{
		comp v1=fft[i+j],v2=fft[i+j+(l>>1)]*gr[t][l+j];
		fft[i+j]=v1+v2;fft[i+j+(l>>1)]=v1-v2;
	}
	double tp=t?1:1.0/s;
	for(int i=0;i<s;i++)a[i]=(comp){fft[i].a*tp,fft[i].b*tp};
}
void mtt(int l,int *a,int *b,int *c)
{
	for(int i=0;i<l;i++)
	{
		c1[i]=(comp){a[i]>>15,0};c2[i]=(comp){a[i]&32767,0};
		c3[i]=(comp){b[i]>>15,0};c4[i]=(comp){b[i]&32767,0};
	}
	dft(l,c1,1);dft(l,c2,1);dft(l,c3,1);dft(l,c4,1);
	for(int i=0;i<l;i++)
	{
		comp v1=c1[i]*c3[i],v2=c1[i]*c4[i]+c2[i]*c3[i],v3=c2[i]*c4[i];
		c1[i]=v1,c2[i]=v2,c3[i]=v3;
	}
	dft(l,c1,0);dft(l,c2,0);dft(l,c3,0);
	for(int i=0;i<l;i++)
	{
		long long f1=c1[i].a+0.4,f2=c2[i].a+0.4,f3=c3[i].a+0.4;
		c[i]=(((f1%mod)<<30)+((f2%mod)<<15)+f3)%mod;
	}
}
int a[N],b[N],c[N];
int v1[N],v2[N],f[N],g[N],rw[N];
void bluestein(int l,int *f,int *g,int w)
{
	if(l<=64)
	{
		for(int i=0;i<l;i++)g[i]=0;
		for(int i=0,a=1;i<l;i++,a=1ll*a*w%mod)
		for(int j=0,b=1;j<l;j++,b=1ll*b*a%mod)
		g[j]=(g[j]+1ll*f[i]*b)%mod;
		return;
	}
	for(int i=0,a=1,b=1;i<=l*2;i++,a=1ll*a*b%mod,b=1ll*b*w%mod)rw[i]=a;
	int le=1;while(le<l*2)le<<=1;
	for(int i=0;i<le;i++)a[i]=b[i]=0;
	for(int i=0;i<l*2;i++)a[i]=rw[i];
	for(int i=1;i<=l;i++)b[i]=1ll*pw(rw[l-i],mod-2)*f[l-i]%mod;
	mtt(le,a,b,c);
	for(int i=0;i<l;i++)g[i]=1ll*c[i+l]*pw(rw[i],mod-2)%mod;
}
void trs(int *v)
{
	for(int i=1,le=1;i<=n;le*=l[i],i++)
	for(int j=0;j<su;j+=le*l[i])
	for(int k=0;k<le;k++)
	{
		for(int p=0;p<l[i];p++)f[p]=v[j+k+p*le];
		bluestein(l[i],f,g,w[i]);
		for(int p=0;p<l[i];p++)v[j+k+p*le]=g[p];
	}
}
void dfs(int d,int x,int v)
{
	if(d==0){v1[x]=v;return;}
	x*=l[d];
	for(int i=0,q=1;i<l[d];i++,q=1ll*q*w[d]%mod)dfs(d-1,x,(v+mod-1ll*q*p[d]%mod)%mod),x++;
}
void dfs2(int d,int x,int y)
{
	if(d==0){v2[x]=(2ll*(sv+mod)-v1[x]-v1[y])%mod;return;}
	x*=l[d],y*=l[d];
	for(int i=0;i<l[d];i++)dfs2(d-1,x+i,y+(l[d]-i)%l[d]);
}
int main()
{
	int mx=0,li=1;
	scanf("%d%d%*d",&n,&mod);
	for(int i=1;i<=n;i++)scanf("%d%d%d",&l[i],&p[i],&w[i]),su*=l[i],sp+=p[i],mx=mx<l[i]?l[i]:mx;
	for(int i=1;i<=n;i++)p[i]=1ll*p[i]*pw(sp,mod-2)%mod;
	while((1<<li)<mx*2)li++;
	init(li);
	dfs(n,0,1);
	for(int i=0;i<su;i++)v1[i]=pw(v1[i],mod-2),sv=(sv+v1[i])%mod;
	trs(v1);
	dfs2(n,0,0);
	for(int i=0;i<su;i++)as=(as+1ll*su*pw(v2[i],mod-2))%mod;
	printf("%d\n",as);
}