|
1 |
| package org.sat4j.reader; |
|
2 |
| |
|
3 |
| import java.io.BufferedReader; |
|
4 |
| import java.io.FileInputStream; |
|
5 |
| import java.io.FileNotFoundException; |
|
6 |
| import java.io.FileReader; |
|
7 |
| import java.io.IOException; |
|
8 |
| import java.io.InputStreamReader; |
|
9 |
| import java.io.LineNumberReader; |
|
10 |
| import java.io.Serializable; |
|
11 |
| import java.math.BigInteger; |
|
12 |
| import java.util.HashMap; |
|
13 |
| import java.util.Map; |
|
14 |
| import java.util.zip.GZIPInputStream; |
|
15 |
| |
|
16 |
| import org.sat4j.core.Vec; |
|
17 |
| import org.sat4j.core.VecInt; |
|
18 |
| import org.sat4j.specs.ContradictionException; |
|
19 |
| import org.sat4j.specs.IProblem; |
|
20 |
| import org.sat4j.specs.ISolver; |
|
21 |
| import org.sat4j.specs.IVec; |
|
22 |
| import org.sat4j.specs.IVecInt; |
|
23 |
| |
|
24 |
| |
|
25 |
| |
|
26 |
| |
|
27 |
| |
|
28 |
| |
|
29 |
| |
|
30 |
| |
|
31 |
| public class OPBReader2005 implements Reader, Serializable { |
|
32 |
| |
|
33 |
| |
|
34 |
| |
|
35 |
| |
|
36 |
| private static final long serialVersionUID = 1L; |
|
37 |
| |
|
38 |
| private ISolver solver; |
|
39 |
| |
|
40 |
| private final IVec<String> decode = new Vec<String>(); |
|
41 |
| |
|
42 |
| private IVecInt lits; |
|
43 |
| |
|
44 |
| private IVec<BigInteger> coeffs; |
|
45 |
| |
|
46 |
| private BigInteger d; |
|
47 |
| |
|
48 |
| private String operator; |
|
49 |
| |
|
50 |
| private final Map<String, Integer> map = new HashMap<String, Integer>(); |
|
51 |
| |
|
52 |
| private final IVecInt objectiveVars = new VecInt(); |
|
53 |
| |
|
54 |
| private final IVec<BigInteger> objectiveCoeffs = new Vec<BigInteger>(); |
|
55 |
| |
|
56 |
| |
|
57 |
| private boolean hasObjFunc = false; |
|
58 |
| |
|
59 |
| |
|
60 |
| |
|
61 |
| |
|
62 |
| |
|
63 |
| |
|
64 |
| |
|
65 |
| |
|
66 |
| |
|
67 |
| |
|
68 |
0
| protected void metaData(int nbvar, int nbconstr) {
|
|
69 |
0
| solver.newVar(nbvar);
|
|
70 |
| } |
|
71 |
| |
|
72 |
| |
|
73 |
| |
|
74 |
| |
|
75 |
0
| protected void beginObjective() {
|
|
76 |
| } |
|
77 |
| |
|
78 |
| |
|
79 |
| |
|
80 |
| |
|
81 |
0
| protected void endObjective() {
|
|
82 |
| assert lits.size() == coeffs.size(); |
|
83 |
| |
|
84 |
| |
|
85 |
| |
|
86 |
| |
|
87 |
| |
|
88 |
| |
|
89 |
| assert lits.size() == coeffs.size(); |
|
90 |
0
| for (int i = 0; i < lits.size(); i++) {
|
|
91 |
0
| objectiveVars.push(lits.get(i));
|
|
92 |
0
| objectiveCoeffs.push(coeffs.get(i));
|
|
93 |
| } |
|
94 |
| } |
|
95 |
| |
|
96 |
| |
|
97 |
| |
|
98 |
| |
|
99 |
0
| protected void beginConstraint() {
|
|
100 |
0
| lits.clear();
|
|
101 |
0
| coeffs.clear();
|
|
102 |
| assert lits.size() == 0; |
|
103 |
| assert coeffs.size() == 0; |
|
104 |
| } |
|
105 |
| |
|
106 |
| |
|
107 |
| |
|
108 |
| |
|
109 |
| |
|
110 |
| |
|
111 |
| |
|
112 |
0
| protected void endConstraint() throws ContradictionException {
|
|
113 |
| |
|
114 |
| assert !(lits.size() == 0); |
|
115 |
| assert !(coeffs.size() == 0); |
|
116 |
| assert lits.size() == coeffs.size(); |
|
117 |
| |
|
118 |
0
| if (operator.equals("<=") || operator.equals("="))
|
|
119 |
0
| solver.addPseudoBoolean(lits, coeffs, false, d);
|
|
120 |
0
| if (operator.equals(">=") || operator.equals("="))
|
|
121 |
0
| solver.addPseudoBoolean(lits, coeffs, true, d);
|
|
122 |
| } |
|
123 |
| |
|
124 |
| |
|
125 |
| |
|
126 |
| |
|
127 |
| |
|
128 |
| |
|
129 |
| |
|
130 |
| |
|
131 |
| |
|
132 |
0
| protected void constraintTerm(BigInteger coeff, String var) {
|
|
133 |
0
| coeffs.push(coeff);
|
|
134 |
0
| Integer id = map.get(var);
|
|
135 |
0
| if (id == null) {
|
|
136 |
| assert map.size()<solver.nVars(); |
|
137 |
0
| map.put(var, id = map.size()+1);
|
|
138 |
| assert id > 0; |
|
139 |
| assert id <= solver.nVars(); |
|
140 |
0
| decode.push(var);
|
|
141 |
| assert decode.size() == id; |
|
142 |
| } |
|
143 |
0
| int lid = ((savedChar == '-') ? -1 : 1) * id.intValue();
|
|
144 |
0
| lits.push(lid);
|
|
145 |
| } |
|
146 |
| |
|
147 |
| |
|
148 |
| |
|
149 |
| |
|
150 |
| |
|
151 |
| |
|
152 |
| |
|
153 |
0
| protected void constraintRelOp(String relop) {
|
|
154 |
0
| operator = relop;
|
|
155 |
| } |
|
156 |
| |
|
157 |
| |
|
158 |
| |
|
159 |
| |
|
160 |
| |
|
161 |
| |
|
162 |
| |
|
163 |
| |
|
164 |
0
| protected void constraintRightTerm(BigInteger val) {
|
|
165 |
0
| d = val;
|
|
166 |
| } |
|
167 |
| |
|
168 |
| transient BufferedReader in; |
|
169 |
| |
|
170 |
| char savedChar; |
|
171 |
| |
|
172 |
| boolean charAvailable = false; |
|
173 |
| |
|
174 |
| boolean eofReached = false; |
|
175 |
| |
|
176 |
| |
|
177 |
| |
|
178 |
| |
|
179 |
| |
|
180 |
| |
|
181 |
0
| private char get() throws IOException {
|
|
182 |
0
| int c;
|
|
183 |
| |
|
184 |
0
| if (charAvailable) {
|
|
185 |
0
| charAvailable = false;
|
|
186 |
0
| return savedChar;
|
|
187 |
| } |
|
188 |
| |
|
189 |
0
| c = in.read();
|
|
190 |
0
| if (c == -1)
|
|
191 |
0
| eofReached = true;
|
|
192 |
| |
|
193 |
0
| return (char) c;
|
|
194 |
| } |
|
195 |
| |
|
196 |
0
| public IVecInt getVars() {
|
|
197 |
0
| return objectiveVars;
|
|
198 |
| } |
|
199 |
| |
|
200 |
0
| public IVec<BigInteger> getCoeffs() {
|
|
201 |
0
| return objectiveCoeffs;
|
|
202 |
| } |
|
203 |
| |
|
204 |
| |
|
205 |
| |
|
206 |
| |
|
207 |
0
| private void putback(char c) {
|
|
208 |
0
| savedChar = c;
|
|
209 |
0
| charAvailable = true;
|
|
210 |
| } |
|
211 |
| |
|
212 |
| |
|
213 |
| |
|
214 |
| |
|
215 |
0
| private boolean eof() {
|
|
216 |
0
| return eofReached;
|
|
217 |
| } |
|
218 |
| |
|
219 |
| |
|
220 |
| |
|
221 |
| |
|
222 |
| |
|
223 |
| |
|
224 |
0
| private void skipSpaces() throws IOException {
|
|
225 |
0
| char c;
|
|
226 |
| |
|
227 |
0
| while (Character.isWhitespace(c = get()))
|
|
228 |
| ; |
|
229 |
| |
|
230 |
0
| putback(c);
|
|
231 |
| } |
|
232 |
| |
|
233 |
| |
|
234 |
| |
|
235 |
| |
|
236 |
| |
|
237 |
| |
|
238 |
| |
|
239 |
0
| public String readWord() throws IOException {
|
|
240 |
0
| StringBuffer s = new StringBuffer();
|
|
241 |
0
| char c;
|
|
242 |
| |
|
243 |
0
| skipSpaces();
|
|
244 |
| |
|
245 |
0
| while (!Character.isWhitespace(c = get()) && !eof())
|
|
246 |
0
| s.append(c);
|
|
247 |
| |
|
248 |
0
| return s.toString();
|
|
249 |
| } |
|
250 |
| |
|
251 |
| |
|
252 |
| |
|
253 |
| |
|
254 |
| |
|
255 |
| |
|
256 |
| |
|
257 |
| |
|
258 |
0
| public void readInteger(StringBuffer s) throws IOException {
|
|
259 |
0
| char c;
|
|
260 |
| |
|
261 |
0
| skipSpaces();
|
|
262 |
0
| s.setLength(0);
|
|
263 |
| |
|
264 |
0
| c = get();
|
|
265 |
0
| if (c == '-' || Character.isDigit(c))
|
|
266 |
0
| s.append(c);
|
|
267 |
| |
|
268 |
| |
|
269 |
0
| while (Character.isDigit(c = get()) && !eof())
|
|
270 |
0
| s.append(c);
|
|
271 |
| |
|
272 |
0
| putback(c);
|
|
273 |
| } |
|
274 |
| |
|
275 |
| |
|
276 |
| |
|
277 |
| |
|
278 |
| |
|
279 |
| |
|
280 |
| |
|
281 |
0
| private boolean readIdentifier(StringBuffer s) throws IOException {
|
|
282 |
0
| char c;
|
|
283 |
| |
|
284 |
0
| s.setLength(0);
|
|
285 |
| |
|
286 |
0
| skipSpaces();
|
|
287 |
| |
|
288 |
| |
|
289 |
0
| c = get();
|
|
290 |
0
| if (eof())
|
|
291 |
0
| return false;
|
|
292 |
| |
|
293 |
0
| if (!Character.isLetter(c) && c != '_') {
|
|
294 |
0
| putback(c);
|
|
295 |
0
| return false;
|
|
296 |
| } |
|
297 |
| |
|
298 |
0
| s.append(c);
|
|
299 |
| |
|
300 |
| |
|
301 |
0
| while (true) {
|
|
302 |
0
| c = get();
|
|
303 |
0
| if (eof())
|
|
304 |
0
| break;
|
|
305 |
| |
|
306 |
0
| if (Character.isLetter(c) || Character.isDigit(c) || c == '_')
|
|
307 |
0
| s.append(c);
|
|
308 |
| else { |
|
309 |
0
| putback(c);
|
|
310 |
0
| break;
|
|
311 |
| } |
|
312 |
| } |
|
313 |
| |
|
314 |
0
| return true;
|
|
315 |
| } |
|
316 |
| |
|
317 |
| |
|
318 |
| |
|
319 |
| |
|
320 |
| |
|
321 |
| |
|
322 |
| |
|
323 |
0
| private String readRelOp() throws IOException {
|
|
324 |
0
| char c;
|
|
325 |
| |
|
326 |
0
| skipSpaces();
|
|
327 |
| |
|
328 |
0
| c = get();
|
|
329 |
0
| if (eof())
|
|
330 |
0
| return null;
|
|
331 |
| |
|
332 |
0
| if (c == '=')
|
|
333 |
0
| return "=";
|
|
334 |
| |
|
335 |
0
| if (c == '>' && get() == '=')
|
|
336 |
0
| return ">=";
|
|
337 |
| |
|
338 |
0
| return null;
|
|
339 |
| } |
|
340 |
| |
|
341 |
| |
|
342 |
| |
|
343 |
| |
|
344 |
| |
|
345 |
| |
|
346 |
| |
|
347 |
| |
|
348 |
0
| private void readMetaData() throws IOException, ParseFormatException {
|
|
349 |
0
| char c;
|
|
350 |
0
| String s;
|
|
351 |
0
| int nbVars, nbConstr;
|
|
352 |
| |
|
353 |
| |
|
354 |
0
| c = get();
|
|
355 |
0
| if (c != '*')
|
|
356 |
0
| throw new ParseFormatException(
|
|
357 |
| "First line of input file should be a comment"); |
|
358 |
| |
|
359 |
0
| s = readWord();
|
|
360 |
0
| if (eof() || !s.equals("#variable="))
|
|
361 |
0
| throw new ParseFormatException(
|
|
362 |
| "First line should contain #variable= as first keyword"); |
|
363 |
| |
|
364 |
0
| nbVars = Integer.parseInt(readWord());
|
|
365 |
| |
|
366 |
0
| s = readWord();
|
|
367 |
0
| if (eof() || !s.equals("#constraint="))
|
|
368 |
0
| throw new ParseFormatException(
|
|
369 |
| "First line should contain #constraint= as second keyword"); |
|
370 |
| |
|
371 |
0
| nbConstr = Integer.parseInt(readWord());
|
|
372 |
| |
|
373 |
| |
|
374 |
0
| in.readLine();
|
|
375 |
| |
|
376 |
| |
|
377 |
0
| metaData(nbVars, nbConstr);
|
|
378 |
| } |
|
379 |
| |
|
380 |
| |
|
381 |
| |
|
382 |
| |
|
383 |
| |
|
384 |
| |
|
385 |
0
| private void skipComments() throws IOException {
|
|
386 |
0
| char c = ' ';
|
|
387 |
| |
|
388 |
| |
|
389 |
| |
|
390 |
0
| while (!eof() && (c = get()) == '*') {
|
|
391 |
0
| in.readLine();
|
|
392 |
| } |
|
393 |
| |
|
394 |
0
| putback(c);
|
|
395 |
| } |
|
396 |
| |
|
397 |
| |
|
398 |
| |
|
399 |
| |
|
400 |
| |
|
401 |
| |
|
402 |
| |
|
403 |
| |
|
404 |
| |
|
405 |
| |
|
406 |
| |
|
407 |
0
| private void readTerm(StringBuffer coeff, StringBuffer var)
|
|
408 |
| throws IOException, ParseFormatException { |
|
409 |
0
| char c;
|
|
410 |
| |
|
411 |
0
| readInteger(coeff);
|
|
412 |
| |
|
413 |
0
| skipSpaces();
|
|
414 |
0
| c = get();
|
|
415 |
0
| if (c != '*')
|
|
416 |
0
| throw new ParseFormatException(
|
|
417 |
| "'*' expected between a coefficient and a variable"); |
|
418 |
| |
|
419 |
0
| if (!readIdentifier(var))
|
|
420 |
0
| throw new ParseFormatException("identifier expected");
|
|
421 |
| } |
|
422 |
| |
|
423 |
| |
|
424 |
| |
|
425 |
| |
|
426 |
| |
|
427 |
| |
|
428 |
| |
|
429 |
| |
|
430 |
0
| private void readObjective() throws IOException, ParseFormatException {
|
|
431 |
0
| char c;
|
|
432 |
0
| StringBuffer var = new StringBuffer();
|
|
433 |
0
| StringBuffer coeff = new StringBuffer();
|
|
434 |
| |
|
435 |
| |
|
436 |
| |
|
437 |
0
| skipSpaces();
|
|
438 |
0
| c = get();
|
|
439 |
0
| if (c != 'm') {
|
|
440 |
| |
|
441 |
0
| putback(c);
|
|
442 |
0
| return;
|
|
443 |
| } |
|
444 |
| |
|
445 |
0
| hasObjFunc = true;
|
|
446 |
0
| if (get() == 'i' && get() == 'n' && get() == ':') {
|
|
447 |
0
| beginObjective();
|
|
448 |
| |
|
449 |
0
| while (!eof()) {
|
|
450 |
0
| readTerm(coeff, var);
|
|
451 |
0
| constraintTerm(new BigInteger(coeff.toString()), var.toString());
|
|
452 |
| |
|
453 |
0
| skipSpaces();
|
|
454 |
0
| c = get();
|
|
455 |
0
| if (c == ';')
|
|
456 |
0
| break;
|
|
457 |
| |
|
458 |
0
| else if (c == '-' || c == '+' || Character.isDigit(c))
|
|
459 |
0
| putback(c);
|
|
460 |
| else |
|
461 |
0
| throw new ParseFormatException(
|
|
462 |
| "unexpected character in objective function"); |
|
463 |
| } |
|
464 |
| |
|
465 |
0
| endObjective();
|
|
466 |
| } else |
|
467 |
0
| throw new ParseFormatException(
|
|
468 |
| "input format error: 'min:' expected"); |
|
469 |
| } |
|
470 |
| |
|
471 |
| |
|
472 |
| |
|
473 |
| |
|
474 |
| |
|
475 |
| |
|
476 |
| |
|
477 |
| |
|
478 |
0
| private void readConstraint() throws IOException, ParseFormatException,
|
|
479 |
| ContradictionException { |
|
480 |
0
| StringBuffer var = new StringBuffer();
|
|
481 |
0
| StringBuffer coeff = new StringBuffer();
|
|
482 |
0
| char c;
|
|
483 |
| |
|
484 |
0
| beginConstraint();
|
|
485 |
| |
|
486 |
0
| while (!eof()) {
|
|
487 |
0
| readTerm(coeff, var);
|
|
488 |
0
| constraintTerm(new BigInteger(coeff.toString()), var.toString());
|
|
489 |
| |
|
490 |
0
| skipSpaces();
|
|
491 |
0
| c = get();
|
|
492 |
0
| if (c == '>' || c == '=') {
|
|
493 |
| |
|
494 |
0
| putback(c);
|
|
495 |
0
| break;
|
|
496 |
0
| } else if (c == '-' || c == '+' || Character.isDigit(c))
|
|
497 |
0
| putback(c);
|
|
498 |
| else { |
|
499 |
0
| throw new ParseFormatException(
|
|
500 |
| "unexpected character in constraint"); |
|
501 |
| } |
|
502 |
| } |
|
503 |
| |
|
504 |
0
| if (eof())
|
|
505 |
0
| throw new ParseFormatException(
|
|
506 |
| "unexpected EOF before end of constraint"); |
|
507 |
| |
|
508 |
0
| String relop;
|
|
509 |
0
| if ((relop = readRelOp()) != null)
|
|
510 |
0
| constraintRelOp(relop);
|
|
511 |
| else |
|
512 |
0
| throw new ParseFormatException(
|
|
513 |
| "unexpected relational operator in constraint"); |
|
514 |
| |
|
515 |
0
| readInteger(coeff);
|
|
516 |
0
| constraintRightTerm(new BigInteger(coeff.toString()));
|
|
517 |
| |
|
518 |
0
| skipSpaces();
|
|
519 |
0
| c = get();
|
|
520 |
0
| if (eof() || c != ';')
|
|
521 |
0
| throw new ParseFormatException(
|
|
522 |
| "semicolon expected at end of constraint"); |
|
523 |
| |
|
524 |
0
| endConstraint();
|
|
525 |
| } |
|
526 |
| |
|
527 |
0
| public OPBReader2005(ISolver solver) {
|
|
528 |
0
| this.solver = solver;
|
|
529 |
0
| lits = new VecInt();
|
|
530 |
0
| coeffs = new Vec<BigInteger>();
|
|
531 |
| } |
|
532 |
| |
|
533 |
| |
|
534 |
| |
|
535 |
| |
|
536 |
| |
|
537 |
| |
|
538 |
| |
|
539 |
| |
|
540 |
| |
|
541 |
0
| public void parse() throws IOException, ParseFormatException,
|
|
542 |
| ContradictionException { |
|
543 |
0
| readMetaData();
|
|
544 |
| |
|
545 |
0
| skipComments();
|
|
546 |
| |
|
547 |
0
| readObjective();
|
|
548 |
| |
|
549 |
| |
|
550 |
0
| while (!eof()) {
|
|
551 |
0
| skipSpaces();
|
|
552 |
0
| if (eof())
|
|
553 |
0
| break;
|
|
554 |
| |
|
555 |
0
| readConstraint();
|
|
556 |
| } |
|
557 |
| } |
|
558 |
| |
|
559 |
0
| public IProblem parseInstance(String filename)
|
|
560 |
| throws FileNotFoundException, ParseFormatException, IOException, |
|
561 |
| ContradictionException { |
|
562 |
| |
|
563 |
0
| if (filename.endsWith(".gz")) {
|
|
564 |
0
| parseInstance(new LineNumberReader(new InputStreamReader(
|
|
565 |
| new GZIPInputStream(new FileInputStream(filename))))); |
|
566 |
| } else { |
|
567 |
0
| parseInstance(new LineNumberReader(new FileReader(filename)));
|
|
568 |
| } |
|
569 |
| |
|
570 |
0
| return solver;
|
|
571 |
| } |
|
572 |
| |
|
573 |
0
| public void parseInstance(LineNumberReader in) throws ParseFormatException,
|
|
574 |
| ContradictionException { |
|
575 |
0
| solver.reset();
|
|
576 |
0
| this.in = in;
|
|
577 |
0
| try {
|
|
578 |
0
| parse();
|
|
579 |
| } catch (IOException e) { |
|
580 |
| |
|
581 |
0
| e.printStackTrace();
|
|
582 |
| } |
|
583 |
| } |
|
584 |
| |
|
585 |
0
| public String decode(int[] model) {
|
|
586 |
0
| StringBuffer stb = new StringBuffer();
|
|
587 |
| |
|
588 |
0
| for (int i = 0; i < model.length; i++) {
|
|
589 |
0
| if (model[i] < 0) {
|
|
590 |
0
| stb.append("-");
|
|
591 |
0
| stb.append(decode.get(-model[i] - 1));
|
|
592 |
| } else { |
|
593 |
0
| stb.append(decode.get(model[i] - 1));
|
|
594 |
| } |
|
595 |
0
| stb.append(" ");
|
|
596 |
| } |
|
597 |
0
| return stb.toString();
|
|
598 |
| } |
|
599 |
| |
|
600 |
0
| public ObjectiveFunction getObjectiveFunction() {
|
|
601 |
0
| if (hasObjFunc)
|
|
602 |
0
| return new ObjectiveFunction(getVars(), getCoeffs());
|
|
603 |
0
| return null;
|
|
604 |
| } |
|
605 |
| |
|
606 |
| } |