-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathGRPH.go
More file actions
69 lines (62 loc) · 1.34 KB
/
Copy pathGRPH.go
File metadata and controls
69 lines (62 loc) · 1.34 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
package main
import (
"bytes"
"fmt"
"io/ioutil"
"strings"
)
type Dna struct {
id string
code string
}
func main() {
rawBytes, _ := ioutil.ReadFile(`data\rosalind_grph.txt`)
fasta := parse(strings.Fields(string(rawBytes)))
//fmt.Println(fasta)
idDnaList := makeList(fasta)
//fmt.Println(idDnaList)
adjList_O3 := adjList(idDnaList, 3)
fmt.Println(stringifyTable(adjList_O3))
}
func adjList(list []Dna, lng int) (table map[string][]string) {
table = make(map[string][]string)
for _, dna1 := range list {
suffix := dna1.code[len(dna1.code)-lng:]
for _, dna2 := range list {
if dna1 == dna2 {
continue
}
prefix := dna2.code[:lng]
if suffix == prefix {
table[dna1.id] = append(table[dna1.id], dna2.id)
}
}
}
return
}
func makeList(fasta []string) (list []Dna) {
for i := 0; i < len(fasta); i += 2 {
record := Dna{id: fasta[i][1:], code: fasta[i+1]}
list = append(list, record)
}
return list
}
func stringifyTable(table map[string][]string) (out string) {
for key, val := range table {
for _, v := range val {
out += fmt.Sprintf("%v %v\n", key, v)
}
}
return
}
func parse(fasta []string) []string {
var buf bytes.Buffer
for _, line := range fasta {
if line[0] == '>' || line == "" {
buf.WriteString("\n" + line + "\n")
} else {
buf.WriteString(line)
}
}
return strings.Fields(buf.String())
}