#include<cstdio>
using namespace std;
const int o=2010;
int n,m,p,h[o],cnt,w[o],s,fa[o],ans,Ans[o][o],tot,d[o],uu,vv,st[o],tp,St[o],Tp,ban;bool vis[o];
inline int qp(int b,int f){int res=1;for(;f;f>>=1,b=b*1ll*b%p) if(f&1) res=res*1ll*b%p;return res;}
struct Edge{int v,p,id;}e[o];
inline void ad(int U,int V,int ID){e[++cnt].v=V;e[cnt].p=h[U];e[h[U]=cnt].id=ID;}
void Dfs(int nw){
	vis[nw]=1;
	for(int i=h[nw];i;i=e[i].p) if(!vis[e[i].v])
		Ans[1][++tot]=fa[e[i].v]=e[i].id,w[e[i].id]=(w[e[i].id]+p-s)%p,d[e[i].v]=d[nw]+1,Dfs(e[i].v);
}
void ddfs(int nw){
	for(int i=h[nw];i;i=e[i].p) if(e[i].id-ban&&!vis[e[i].v]) vis[e[i].v]=1,Ans[ans][++tot]=e[i].id;
}
inline void slv(int x,int y,int w){
	Ans[++ans][0]=w;tot=0;
	for(int i=1;i<=n;++i) vis[i]=0;
	for(int i=1;i<=Tp;++i) vis[e[St[i]*2-1].v]=1,vis[e[St[i]*2].v]=1;
	for(int i=1;i<=Tp;++i) if(St[i]^y) Ans[ans][++tot]=St[i];
	ban=y;
	for(int i=1;i<=n;++i) if(!vis[i]) ddfs(i);
}
void dfs(int nw){
	st[++tp]=nw;
	for(int i=h[nw];i;i=e[i].p) if(e[i].id^fa[nw]){
		if(fa[e[i].v]==e[i].id){
			dfs(e[i].v);
			if(uu==nw) uu=0,vv=0;
		}
		if(d[e[i].v]<d[nw]&&(!uu||d[e[i].v]*2+1<d[uu]+d[vv])){
			uu=st[d[e[i].v]];vv=st[d[e[i].v]+1];St[Tp=1]=e[i].id;
			for(int j=d[e[i].v]+1;j<=d[nw];++j) St[++Tp]=fa[st[j]];
		}
		if(w[e[i].id]&&uu)
			slv(e[i].id,fa[vv],w[e[i].id]),slv(fa[vv],e[i].id,p-w[e[i].id]),
			w[fa[vv]]=(w[fa[vv]]+w[e[i].id])%p,w[e[i].id]=0;
	}
	--tp;
}
int main(){
	scanf("%d%d%d",&n,&m,&p);
	for(int i=1,u,v;i<=m;++i) scanf("%d%d%d",&u,&v,&w[i]),s=(s+w[i])%p,ad(u,v,i),ad(v,u,i);
	Ans[ans=1][0]=s=s*1ll*qp(n-1,p-2)%p;
	Dfs(d[1]=1);dfs(1);
	for(int i=1;i<=m;++i) if(w[i]){printf("-1\n");return 0;}
	printf("%d\n",ans);
	for(int i=1,j;i<=ans;++i,putchar('\n')) for(printf("%d ",Ans[i][0]),j=1;j<n;++j) printf("%d ",Ans[i][j]);
	return 0;
}