176
|
1 using Implab;
|
|
2 using System;
|
|
3 using System.Collections.Generic;
|
|
4 using System.Linq;
|
|
5
|
|
6 namespace Implab.Automaton {
|
|
7 public class DFATable : IDFATableBuilder {
|
|
8 int m_stateCount;
|
|
9 int m_symbolCount;
|
|
10 int m_initialState;
|
|
11
|
|
12 readonly HashSet<int> m_finalStates = new HashSet<int>();
|
|
13 readonly HashSet<AutomatonTransition> m_transitions = new HashSet<AutomatonTransition>();
|
|
14
|
|
15
|
|
16 #region IDFADefinition implementation
|
|
17
|
|
18 public bool IsFinalState(int s) {
|
|
19 Safe.ArgumentInRange(s, 0, m_stateCount, "s");
|
|
20
|
|
21 return m_finalStates.Contains(s);
|
|
22 }
|
|
23
|
|
24 public IEnumerable<int> FinalStates {
|
|
25 get {
|
|
26 return m_finalStates;
|
|
27 }
|
|
28 }
|
|
29
|
|
30 public int StateCount {
|
|
31 get { return m_stateCount; }
|
|
32 }
|
|
33
|
|
34 public int AlphabetSize {
|
|
35 get { return m_symbolCount; }
|
|
36 }
|
|
37
|
|
38 public int InitialState {
|
|
39 get { return m_initialState; }
|
|
40 }
|
|
41
|
|
42 #endregion
|
|
43
|
|
44 public void SetInitialState(int s) {
|
|
45 Safe.ArgumentAssert(s >= 0, "s");
|
|
46 m_initialState = s;
|
|
47 }
|
|
48
|
|
49 public void MarkFinalState(int state) {
|
|
50 m_finalStates.Add(state);
|
|
51 }
|
|
52
|
|
53 public void Add(AutomatonTransition item) {
|
|
54 Safe.ArgumentAssert(item.s1 >= 0, "item");
|
|
55 Safe.ArgumentAssert(item.s2 >= 0, "item");
|
|
56 Safe.ArgumentAssert(item.edge >= 0, "item");
|
|
57
|
|
58 m_stateCount = Math.Max(m_stateCount, Math.Max(item.s1, item.s2) + 1);
|
|
59 m_symbolCount = Math.Max(m_symbolCount, item.edge);
|
|
60
|
|
61 m_transitions.Add(item);
|
|
62 }
|
|
63
|
|
64 public void Clear() {
|
|
65 m_stateCount = 0;
|
|
66 m_symbolCount = 0;
|
|
67 m_finalStates.Clear();
|
|
68 m_transitions.Clear();
|
|
69 }
|
|
70
|
|
71 public bool Contains(AutomatonTransition item) {
|
|
72 return m_transitions.Contains(item);
|
|
73 }
|
|
74
|
|
75 public void CopyTo(AutomatonTransition[] array, int arrayIndex) {
|
|
76 m_transitions.CopyTo(array, arrayIndex);
|
|
77 }
|
|
78
|
|
79 public bool Remove(AutomatonTransition item) {
|
180
|
80 return m_transitions.Remove(item);
|
176
|
81 }
|
|
82
|
|
83 public int Count {
|
|
84 get {
|
|
85 return m_transitions.Count;
|
|
86 }
|
|
87 }
|
|
88
|
|
89 public bool IsReadOnly {
|
|
90 get {
|
|
91 return false;
|
|
92 }
|
|
93 }
|
|
94
|
|
95 public IEnumerator<AutomatonTransition> GetEnumerator() {
|
|
96 return m_transitions.GetEnumerator();
|
|
97 }
|
|
98
|
|
99 System.Collections.IEnumerator System.Collections.IEnumerable.GetEnumerator() {
|
|
100 return GetEnumerator();
|
|
101 }
|
|
102
|
|
103 public int[,] CreateTransitionTable() {
|
|
104 var table = new int[StateCount,AlphabetSize];
|
|
105
|
|
106 for (int i = 0; i < StateCount; i++)
|
|
107 for (int j = 0; i < AlphabetSize; j++)
|
178
|
108 table[i, j] = AutomatonConst.UNREACHABLE_STATE;
|
176
|
109
|
|
110 foreach (var t in this)
|
|
111 table[t.s1,t.edge] = t.s2;
|
|
112
|
|
113 return table;
|
|
114 }
|
|
115
|
|
116 public bool[] CreateFinalStateTable() {
|
|
117 var table = new bool[StateCount];
|
|
118
|
|
119 foreach (var s in FinalStates)
|
|
120 table[s] = true;
|
|
121
|
|
122 return table;
|
|
123 }
|
|
124
|
|
125 /// <summary>Формирует множества конечных состояний перед началом работы алгоритма минимизации.</summary>
|
|
126 /// <remarks>
|
|
127 /// В процессе построения минимального автомата требуется разделить множество состояний,
|
|
128 /// на два подмножества - конечные состояния и все остальные, после чего эти подмножества
|
|
129 /// будут резделены на более мелкие. Иногда требуется гарантировать различия конечных сосотяний,
|
|
130 /// для этого необходимо переопределить даннцю фукнцию, для получения множеств конечных состояний.
|
|
131 /// </remarks>
|
|
132 /// <returns>The final states.</returns>
|
|
133 protected virtual IEnumerable<HashSet<int>> GroupFinalStates() {
|
|
134 return new HashSet<int>[] { m_finalStates };
|
|
135 }
|
|
136
|
|
137 protected void Optimize(
|
|
138 IDFATableBuilder optimalDFA,
|
|
139 IDictionary<int,int> alphabetMap,
|
|
140 IDictionary<int,int> stateMap
|
|
141 ) {
|
|
142 Safe.ArgumentNotNull(optimalDFA, "dfa");
|
|
143 Safe.ArgumentNotNull(alphabetMap, "alphabetMap");
|
|
144 Safe.ArgumentNotNull(stateMap, "stateMap");
|
|
145
|
|
146
|
|
147 var setComparer = new CustomEqualityComparer<HashSet<int>>(
|
|
148 (x, y) => x.SetEquals(y),
|
|
149 s => s.Sum(x => x.GetHashCode())
|
|
150 );
|
|
151
|
|
152 var optimalStates = new HashSet<HashSet<int>>(setComparer);
|
|
153 var queue = new HashSet<HashSet<int>>(setComparer);
|
|
154
|
|
155 // получаем конечные состояния, сгруппированные по маркерам
|
|
156 optimalStates.UnionWith(
|
|
157 GroupFinalStates()
|
|
158 );
|
|
159
|
|
160 var state = new HashSet<int>(
|
|
161 Enumerable
|
|
162 .Range(0, m_stateCount - 1)
|
|
163 .Where(i => !m_finalStates.Contains(i))
|
|
164 );
|
|
165
|
|
166 optimalStates.Add(state);
|
|
167 queue.Add(state);
|
|
168
|
|
169 var rmap = m_transitions
|
|
170 .GroupBy(t => t.s2)
|
180
|
171 .ToDictionary(
|
176
|
172 g => g.Key, // s2
|
180
|
173 g => g.GroupBy(t => t.edge, t => t.s1).ToDictionary(p => p.Key)
|
176
|
174 );
|
|
175
|
|
176 while (queue.Count > 0) {
|
|
177 var stateA = queue.First();
|
|
178 queue.Remove(stateA);
|
|
179
|
|
180 for (int c = 0; c < m_symbolCount; c++) {
|
|
181 var stateX = new HashSet<int>();
|
|
182 foreach(var a in stateA)
|
180
|
183 stateX.UnionWith(rmap[a][c]); // all states from wich the symbol 'c' leads to the state 'a'
|
176
|
184
|
|
185 foreach (var stateY in optimalStates.ToArray()) {
|
|
186 if (stateX.Overlaps(stateY) && !stateY.IsSubsetOf(stateX)) {
|
|
187 var stateR1 = new HashSet<int>(stateY);
|
|
188 var stateR2 = new HashSet<int>(stateY);
|
|
189
|
|
190 stateR1.IntersectWith(stateX);
|
|
191 stateR2.ExceptWith(stateX);
|
|
192
|
|
193 optimalStates.Remove(stateY);
|
|
194 optimalStates.Add(stateR1);
|
|
195 optimalStates.Add(stateR2);
|
|
196
|
|
197 if (queue.Contains(stateY)) {
|
|
198 queue.Remove(stateY);
|
|
199 queue.Add(stateR1);
|
|
200 queue.Add(stateR2);
|
|
201 } else {
|
|
202 queue.Add(stateR1.Count <= stateR2.Count ? stateR1 : stateR2);
|
|
203 }
|
|
204 }
|
|
205 }
|
|
206 }
|
|
207 }
|
|
208
|
|
209 // карта получения оптимального состояния по соотвествующему ему простому состоянию
|
|
210 var nextState = 0;
|
|
211 foreach (var item in optimalStates) {
|
|
212 var id = nextState++;
|
|
213 foreach (var s in item)
|
|
214 stateMap[s] = id;
|
|
215 }
|
|
216
|
|
217 // получаем минимальный алфавит
|
|
218 // входные символы не различимы, если Move(s,a1) == Move(s,a2), для любого s
|
|
219 // для этого используем алгоритм кластеризации, сначала
|
|
220 // считаем, что все символы не различимы
|
|
221
|
|
222 var minClasses = new HashSet<HashSet<int>>(setComparer);
|
|
223 var alphaQueue = new Queue<HashSet<int>>();
|
|
224 alphaQueue.Enqueue(new HashSet<int>(Enumerable.Range(0,AlphabetSize)));
|
|
225
|
|
226 // для всех состояний, будем проверять каждый класс на различимость,
|
|
227 // т.е. символы различимы, если они приводят к разным состояниям
|
|
228 for (int s = 0 ; s < optimalStates.Count; s++) {
|
|
229 var newQueue = new Queue<HashSet<int>>();
|
|
230
|
|
231 foreach (var A in alphaQueue) {
|
|
232 // классы из одного символа делить бесполезно, переводим их сразу в
|
|
233 // результирующий алфавит
|
|
234 if (A.Count == 1) {
|
|
235 minClasses.Add(A);
|
|
236 continue;
|
|
237 }
|
|
238
|
|
239 // различаем классы символов, которые переводят в различные оптимальные состояния
|
|
240 // optimalState -> alphaClass
|
|
241 var classes = new Dictionary<int, HashSet<int>>();
|
|
242
|
|
243 foreach (var term in A) {
|
|
244 // ищем все переходы класса по символу term
|
|
245 var res = m_transitions.Where(t => stateMap[t.s1] == s && t.edge == term).Select(t => stateMap[t.s2]).ToArray();
|
|
246
|
|
247 var s2 = res.Length > 0 ? res[0] : -1;
|
|
248
|
|
249 HashSet<int> a2;
|
|
250 if (!classes.TryGetValue(s2, out a2)) {
|
|
251 a2 = new HashSet<int>();
|
|
252 newQueue.Enqueue(a2);
|
|
253 classes[s2] = a2;
|
|
254 }
|
|
255 a2.Add(term);
|
|
256 }
|
|
257 }
|
|
258
|
|
259 if (newQueue.Count == 0)
|
|
260 break;
|
|
261 alphaQueue = newQueue;
|
|
262 }
|
|
263
|
|
264 // после окончания работы алгоритма в очереди останутся минимальные различимые классы
|
|
265 // входных символов
|
|
266 foreach (var A in alphaQueue)
|
|
267 minClasses.Add(A);
|
|
268
|
|
269 // построение отображения алфавитов входных символов.
|
|
270 // поскольку символ DFAConst.UNCLASSIFIED_INPUT может иметь
|
|
271 // специальное значение, тогда сохраним минимальный класс,
|
|
272 // содержащий этот символ на томже месте.
|
|
273
|
|
274 var nextCls = 0;
|
|
275 foreach (var item in minClasses) {
|
178
|
276 if (nextCls == AutomatonConst.UNCLASSIFIED_INPUT)
|
176
|
277 nextCls++;
|
|
278
|
|
279 // сохраняем DFAConst.UNCLASSIFIED_INPUT
|
178
|
280 var cls = item.Contains(AutomatonConst.UNCLASSIFIED_INPUT) ? AutomatonConst.UNCLASSIFIED_INPUT : nextCls;
|
176
|
281
|
|
282 foreach (var a in item)
|
|
283 alphabetMap[a] = cls;
|
|
284
|
|
285 nextCls++;
|
|
286 }
|
|
287
|
|
288 // построение автомата
|
|
289 optimalDFA.SetInitialState(stateMap[m_initialState]);
|
|
290
|
|
291 foreach (var sf in m_finalStates.Select(s => stateMap[s]).Distinct())
|
|
292 optimalDFA.MarkFinalState(sf);
|
|
293
|
|
294 foreach (var t in m_transitions.Select(t => new AutomatonTransition(stateMap[t.s1],stateMap[t.s2],alphabetMap[t.edge])).Distinct())
|
|
295 optimalDFA.Add(t);
|
|
296 }
|
|
297
|
|
298 protected void PrintDFA<TInput, TState>(IAlphabet<TInput> inputAlphabet, IAlphabet<TState> stateAlphabet) {
|
|
299 Safe.ArgumentNotNull(inputAlphabet, "inputAlphabet");
|
|
300 Safe.ArgumentNotNull(stateAlphabet, "stateAlphabet");
|
|
301
|
|
302 foreach(var t in m_transitions)
|
|
303 Console.WriteLine(
|
|
304 "[{0}] -{{{1}}}-> [{2}]{3}",
|
|
305 String.Join(",", stateAlphabet.GetSymbols(t.s1)),
|
|
306 String.Join("", inputAlphabet.GetSymbols(t.edge)),
|
|
307 String.Join(",", stateAlphabet.GetSymbols(t.s2)),
|
|
308 m_finalStates.Contains(t.s2) ? "$" : ""
|
|
309 );
|
|
310 }
|
|
311
|
|
312 }
|
|
313 }
|