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
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
#include<stdio.h>
#include<stdlib.h>
#include<message.h>
#include<vector>
#include"poszukiwania.h"

//#define DEBUG

#ifdef DEBUG
#define log(format, ...)     fprintf(stderr, format "\n", ##__VA_ARGS__)
#else
#define log(...)
#endif


using namespace std;

int from=1;               //!< Zakres początkowy danych (+1)
int sLen;               //!< długość sygnału
int mLen;               //!< długość ciągu matematycznego
int chckLen;            //!< Długość danych do sprawdzania (mLen-sLen)

inline int x(int pos) { return SignalAt(pos+1); }
inline int y(int pos) { return SeqAt(pos+from); }

inline bool check(int pos) {
    pos+=from-1;
    for(int j=1;j<=sLen;++j) {
        if(SignalAt(j)!=SeqAt(j+pos)) return false;   // jeżeli jest miejsce niepazujące, to zwracamy false
    }
    return true;    // jeżeli wszystkie elementy się zgadzają
}

/* Computing of the maximal suffix for <= */
int maxSuf(int m, int *p) {
   int ms, j, k;
   char a, b;

   ms = -1;
   j = 0;
   k = *p = 1;
   while (j + k < m) {
      a = x(j + k);
      b = x(ms + k);
      if (a < b) {
         j += k;
         k = 1;
         *p = j - ms;
      }
      else
         if (a == b)
            if (k != *p)
               ++k;
            else {
               j += *p;
               k = 1;
            }
         else { /* a > b */
            ms = j;
            j = ms + 1;
            k = *p = 1;
         }
   }
   return(ms);
}

/* Computing of the maximal suffix for >= */
int maxSufTilde(int m, int *p) {
   int ms, j, k;
   char a, b;

   ms = -1;
   j = 0;
   k = *p = 1;
   while (j + k < m) {
      a = x(j + k);
      b = x(ms + k);
      if (a > b) {
         j += k;
         k = 1;
         *p = j - ms;
      }
      else
         if (a == b)
            if (k != *p)
               ++k;
            else {
               j += *p;
               k = 1;
            }
         else { /* a < b */
            ms = j;
            j = ms + 1;
            k = *p = 1;
         }
   }
   return(ms);
}

int cmp(int k, int len) {
    for(int i=0;i<len;++i) if(x(i)!=x(k+i)) return false;
    return true;
}

int MAX(int a, int b) { return a>b?a:b; }

/* Two Way string matching algorithm. */
int TW() {
   int i, j, ell, memory, p, per, q;

   int m=sLen;
   int n=mLen;
   int res=0;

   /* Preprocessing */
   i = maxSuf(m, &p);
   j = maxSufTilde(m, &q);
   if (i > j) {
      ell = i;
      per = p;
   }
   else {
      ell = j;
      per = q;
   }

   /* Searching */
   if ( cmp(per, ell + 1) ) {
      j = 0;
      memory = -1;
      while (j <= n - m) {
         i = MAX(ell, memory) + 1;
         while (i < m && x(i) == y(i + j) )
            ++i;
         if (i >= m) {
            i = ell;
            while (i > memory && x(i) == y(i + j) )
               --i;
            if (i <= memory) ++res;
            j += per;
            memory = m - per - 1;
         }
         else {
            j += (i - ell);
            memory = -1;
         }
      }
   }
   else {
      per = MAX(ell + 1, m - ell - 1) + 1;
      j = 0;
      while (j <= n - m) {
         i = ell + 1;
         while (i < m && x(i) == y(i + j))
            ++i;
         if (i >= m) {
            i = ell;
            while (i >= 0 && x(i) == y(i + j) )
               --i;
            if (i < 0) ++res;
            j += per;
         }
         else
            j += (i - ell);
      }
   }
   return res;
}

int main(int argc, char** argv) {
	sLen=SignalLength();
	mLen=SeqLength();

    // Każy wykonuje 1/n pracy i zrzuca wyniki
    if(mLen<NumberOfNodes()*100)
    //if(false)
    {  // jeżeli jest mało zadań, to wykonujemy na jednym komputerze
        if(MyNodeId()!=0) return 0;
        log("Small data node");
        // prostą metodą
        int chckLen=mLen-sLen;
        int res=0;
        for(int i=0;i<chckLen;++i) {
            if(check(i)) ++res;
        }
        printf("%d\n",res); // wypisujemy wsumaryczny wynik
        return 0;
    }
    // tryb pracy na wielu komputerach
    int part=(mLen-SignalLength()+NumberOfNodes()-1)/NumberOfNodes();  // jedna paczka danych
    from=MyNodeId()*part;
    mLen=part+SignalLength()-1;
    if(mLen+from>SeqLength()) mLen=SeqLength()-from;
    log("Multinode from %d to %d (%d) of %d", from, mLen+from-1, mLen, (int)SeqLength());
    from=from+1;

    int res;
    if(sLen<10)
    {   // to wykonujemy zwykłą metodą, bo narzut będzie zbyt duży
        log("Small pattern mode");
        int chckLen=mLen-sLen;
        res=0;
        for(int i=0;i<chckLen;++i) {
            if(check(i)) ++res;
        }
    } else {
        // uruchomienie kompletnego algorytmu "Two Way"
        res=TW();
    }
    if(MyNodeId()>0) {  // to odsyłamy wyniki to 0-rowego
        log("Node %d sending result to core", MyNodeId());
        PutInt(0,res);
        Send(0);
        return 0;           // i kończym prace
    }
    // jak jesteśmy zerowym, to zbieramy wyniki od pozostałych
    for(int i=1;i<NumberOfNodes();++i) {
        Receive(i);
        int cres=GetInt(i);
        res+=cres;
        log("Received summary from node: %d = %d",i,cres);
    }

	printf("%d\n",res); // wypisujemy wsumaryczny wynik

	return 0;
}