diff options
| -rw-r--r-- | main.cpp | 10 |
1 files changed, 5 insertions, 5 deletions
@@ -44,7 +44,7 @@ bitset<k> test(string & probe) { bitset<k> bs; double avg[7]; for (int i = 0; i < 7; ++i) { - for (int j = 0; j < 2; ++j) { + for (int j = 0; j < 10; ++j) { int r = rand()%20; avg[i] += sequenceCmp(probe, vrt[i][r]); } @@ -71,7 +71,7 @@ void analyze_spike_sequences(){ } int cnt = 0; for (int i = 0; i < 7; ++i) { - for (int j = 0; j <= vrt[i][0].size()-150; ++j) { + for (int j = 0; j <= vrt[i][0].size()-150; j += rand()%100) { string probe(begin(vrt[i][0])+j, begin(vrt[i][0])+j+150); bitset<k> bs = test(probe); cout << probe << ' ' << cnt++ << endl; @@ -95,7 +95,7 @@ void analyze_spike_sequences(){ int retx = 0, rety = 0; //find best dp coverage - for(int i = 1; i <= 3; i++) + for(int i = 3; i <= 3; i++) for(int j = 0; j < (1 << k); j++){ if(!dp[i][j].empty() && __builtin_popcount(j) > __builtin_popcount(rety)){ retx = i, rety = j; @@ -162,10 +162,10 @@ int main(){ ios::sync_with_stdio(0); cin.tie(0); if (fopen("spike_sequences.txt", "r")) { - analyze_spike_sequences(); if (P[0] != "") mkscores(); + else analyze_spike_sequences(); return 0; } answer(); return 0; -}
\ No newline at end of file +} |