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
|
#include "parser.h"
#include "printer.h"
#include "polynomial.h"
#include "division.h"
#include "buchberger.h"
#include "wallideal.h"
#include "lp.h"
#include "reversesearch.h"
#include "termorder.h"
#include "ep_standard.h"
#include "ep_xfig.h"
#include "gfanapplication.h"
#include "symmetry.h"
class PolynomialSetUnionApplication : public GFanApplication
{
SimpleOption optionSpecialMode;
class PolynomialCmp
{
public:
bool operator()(const Polynomial &a, const Polynomial &b)const
{
if(a.getMarked().m.exponent.sum()<b.getMarked().m.exponent.sum())return true;
if(a.getMarked().m.exponent.sum()>b.getMarked().m.exponent.sum())return false;
return LexicographicTermOrder()(a.getMarked().m.exponent,b.getMarked().m.exponent);
}
};
public:
const char *helpText()
{
return "This program computes the union of a list of polynomial sets given as input. The polynomials must all belong to the same ring. The ring is specified an the input. After this follows the list of polynomial sets.\n";
}
PolynomialSetUnionApplication():
optionSpecialMode("-s","Sort output by degree.\n")
{
registerOptions();
}
char *name()
{
return "_polynomialsetunion";
}
int main()
{
/* SymmetryGroup g(16);
// g.computeClosure(StringParser("(1,2,3,0,5,6,7,4,9,10,11,8,13,14,15,12)").parseIntegerVector());
// g.computeClosure(StringParser("(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)").parseIntegerVector());
g.computeClosure(StringParser("(1,0,2,3,5,4,6,7,9,8,10,11,13,12,14,15)").parseIntegerVector());
g.computeClosure(StringParser("(0,2,1,3,4,6,5,7,8,10,9,11,12,14,13,15)").parseIntegerVector());
g.computeClosure(StringParser("(0,1,3,2,4,5,7,6,8,9,11,10,12,13,15,14)").parseIntegerVector());
g.computeClosure(StringParser("(0,2,1,3,4,6,5,7,8,10,9,11,12,14,13,15)").parseIntegerVector());
g.computeClosure(StringParser("(4,5,6,7,0,1,2,3,8,9,10,11,12,13,14,15)").parseIntegerVector());
g.computeClosure(StringParser("(8,9,10,11,4,5,6,7,0,1,2,3,12,13,14,15)").parseIntegerVector());
g.computeClosure(StringParser("(12,13,14,15,4,5,6,7,8,9,10,11,0,1,2,3)").parseIntegerVector());
g.computeClosure(StringParser("(8,9,10,11,4,5,6,7,0,1,2,3,12,13,14,15)").parseIntegerVector());
g.print(Stderr);
*/
FileParser P(Stdin);
PolynomialRing r=P.parsePolynomialRing();
int c=P.nextNonBlank();
bool first=true;
assert(P.isLeftBracket(c));
PolynomialSet s=P.parsePolynomialSet(r);
while(','==c || first)
{
PolynomialSet temp=P.parsePolynomialSet(r);
// if(optionSpecialMode.getValue())temp.markAndScale(LexicographicTermOrder());
s.unionSet(temp);
c=P.nextNonBlank();
first=false;
}
assert(P.isRightBracket(c));
/* PolynomialSet S;
for(PolynomialSet::const_iterator i=s.begin();i!=s.end();i++)
{
S.unionPolynomial((*i));
// S.unionPolynomial(g.computeUniqueRepresentative(*i));
}
*/
if(optionSpecialMode.getValue())s.sort(PolynomialCmp());
AsciiPrinter(Stdout).printPolynomialRing(r);
AsciiPrinter(Stdout).printNewLine();
AsciiPrinter(Stdout).printPolynomialSet(s);
/*
PolynomialSet total;
for(PolynomialSet::const_iterator i=S.begin();i!=S.end();i++)
{
AsciiPrinter(Stdout).printPolynomial(*i);
PolynomialSet a;
for(SymmetryGroup::ElementContainer::const_iterator j=g.elements.begin();j!=g.elements.end();j++)
{
a.unionPolynomial(g.permutePolynomial(*i,*j));
total.unionPolynomial(g.permutePolynomial(*i,*j));
}
fprintf(Stdout,"\nelments in class: %i\n",a.size());
fprintf(Stdout,"Total so far: %i\n\n",total.size());
}
*/
return 0;
}
};
static PolynomialSetUnionApplication theApplication;
|
