1
#include<cstdio>
2
#include<cstring>
3
#include<cmath>
4
#include<algorithm>
5
using namespace std;
6
struct edge{
7
    int END,next;
8
    bool b[32];
9
}v[30005];
10
int first[105],p,t[105],du[105];
11
double f[105][35][2];
12
double a[105][105],x[105];
13
void add(int a,int b,int c){
14
    int i=1;
15
    du[a]++;
16
    while(c){
17
        v[p].b[i]=c&1;
18
        i++;c>>=1;
19
    }
20
    v[p].END=b;v[p].next=first[a];
21
    first[a]=p++;
22
}
23
int n,m;
24
void gauss(){
25
    int im,num=1;
26
    for(int k=1;k<=n;k++,num++){
27
        im=k;
28
        for(int i=k+1;i<=n;i++){
29
            if(fabs(a[i][k])>fabs(a[im][k]))
30
                im=i;
31
        }
32
        if(im!=k){
33
            for(int i=k;i<=n+1;i++){
34
                swap(a[num][i],a[im][i]);
35
            }
36
        }
37
        if(!a[num][k]){
38
            num--;continue;
39
        }
40
        for(int i=num+1;i<=n;i++){
41
            if(!a[num][k])continue;
42
            long double t=a[i][k]/a[num][k];
43
            for(int j=k;j<=n+1;j++){
44
                a[i][j]-=t*a[k][j];
45
            }
46
        }
47
    }
48
    for(int i=n;i>=1;i--){
49
        for(int j=n;j>=i+1;j--){
50
            a[i][n+1]-=a[i][j]*x[j];
51
        }
52
        x[i]=a[i][n+1]/a[i][i];
53
    }
54
}
55
int main()
56
{
57
    //freopen("xorpath.in", "r", stdin);
58
    //freopen("xorpath.out", "w", stdout);
59
    memset(first,-1,sizeof(first));
60
    scanf("%d%d",&n,&m);
61
    for(int a1,b,c,i=1;i<=m;i++){
62
        scanf("%d%d%d",&a1,&b,&c);
63
        add(a1,b,c);
64
        if(a1!=b) add(b,a1,c);
65
    }
66
    //for(int i=first[n];i!=-1;i=v[i].next)
67
     //   rdu[v[i].END]--;
68
    long double ans=0;
69
    for(int k=1;k<=31;k++){
70
        memset(a,0,sizeof(a));
71
        for(int i=1;i<n;i++){
72
            a[i][i]=1;
73
            memset(t,0,sizeof(t));
74
            for(int j=first[i];j!=-1;j=v[j].next){
75
                //if(v[j].END==n)continue;
76
                if(v[j].b[k]) a[i][v[j].END]+=double(1)/du[i],a[i][n+1]+=(double)1/du[i];
77
                else a[i][v[j].END]-=double(1)/du[i];
78
            }
79
        }
80
        a[n][n]=1;
81
        gauss();
82
        ans+=x[1]*(1<<(k-1));
83
    }
84
    printf("%.3lf",(double)ans);
85
}