PuzzleSolver/Source/functions/solve/structure.cpp

310 lines
12 KiB
C++
Raw Normal View History

#include "../../header.h"
bool SetBestOrMoreLayersArithmetical(vector<LogEntry>& log, qualityVector& cqVector);
2017-12-22 22:55:55 +01:00
void calculateTrueDestructionPower(vector<LogEntry>& log, Puzzle& puzzleMat, float Layerworth);
void cut(vector<LogEntry>& log, int& cutID);
float capLogElements(vector<LogEntry>& log);
void CalculateNewCombinedQuality(vector<LogEntry>& log, qualityVector& qVector, qualityVector& cqVector);
bool next(vector<LogEntry>& log,Puzzle& puzzleMat)
{
//last log element is set, create new log element or log not yet started
2017-12-03 20:12:32 +01:00
if(!(log.size()) || log.back().isSet())
if((puzzleMat.allSet())) return false; //puzzle solved
else createNextLogElement(log,puzzleMat);
2018-01-07 20:08:50 +01:00
//last log element is empty, backtrack
2018-01-07 20:08:50 +01:00
else if(!(log.back().PieceCollector.size()))
{
if(!(backtrack(log,puzzleMat)))
return false;
}
//case last log element has multiple entries
else if(log.back().PieceCollector.size() > 1)
2017-12-22 22:55:55 +01:00
//moreLayers is 0, setbest is 1
if (SetBestOrMoreLayersArithmetical(log, puzzleMat.combinedQualityVector)) setsolution(log, puzzleMat);
else solve(log, puzzleMat);
2018-01-07 20:08:50 +01:00
//case last log exactly one solution
else if(log.back().PieceCollector.size() == 1)
if(log.back().hasRandomed())
if(log.back().abstractionLevel < 2)//do 2 at least two best abstractions to check if part is okay
solve(log,puzzleMat);
else
setsolution(log,puzzleMat);
else
setsolution(log,puzzleMat);
return true;
}
void createNextLogElement(vector<LogEntry>& log, Puzzle& puzzleMat)
{
log.emplace_back(LogEntry(coor(0, 0)));
log.back().myCoor = calculateNextCoor(log, puzzleMat);
2017-12-22 23:09:23 +01:00
puzzleMat.dp.DestructionOfSurrounding(log.back().myCoor);//calculate dp from surrounding
2018-01-07 20:08:50 +01:00
//get all not set pieces
for(auto it:puzzleMat.p_myBox)
2018-01-07 20:08:50 +01:00
if(!it->set)
log.back().PieceCollector.emplace_back(pair<float,Part*>(0,it));
solve(log,puzzleMat);
}
coor calculateNextCoor(vector<LogEntry>& log, Puzzle& puzzleMat)
{
//level 1:
//go left to right, then increase current row
2017-12-03 20:12:32 +01:00
if (log.size() == 1)
return {0,0};
2017-12-03 20:12:32 +01:00
2018-01-07 20:08:50 +01:00
unsigned int col= log.rbegin()[1].myCoor.col;
unsigned int row= log.rbegin()[1].myCoor.row;
2017-12-03 20:12:32 +01:00
2018-01-07 20:08:50 +01:00
if(row<puzzleMat.getSizeAsCoor().row-1) row++;
else if(col<puzzleMat.getSizeAsCoor().col-1){ row=0; col++;}
return {col,row};
}
2018-01-07 20:08:50 +01:00
void solve(vector<LogEntry>& log,Puzzle& puzzleMat)
{
log.back().abstractionLevel = puzzleMat.dp.getNextAbstractionLayer(log.back().myCoor,log.back().abstractionLevel); //sets in abstractionLevel
//status(log,p_Box,puzzleMat);
switch(log.back().abstractionLevel)
{
2018-01-07 20:08:50 +01:00
case 0://pömpel
2017-12-22 22:47:14 +01:00
puzzleMat.a1.EvaluateQuality(log.back().myCoor, log.back().PieceCollector);
break;
2018-01-07 20:08:50 +01:00
case -1://random
setsolution(log,puzzleMat);
return;
default:
break;
}
float worth = capLogElements(log);
calculateTrueDestructionPower(log,puzzleMat, worth);
CalculateNewCombinedQuality(log, log.back().PieceCollector, puzzleMat.combinedQualityVector);
}
//removes from box and makes log "set"
void setsolution(vector<LogEntry>& log, Puzzle& puzzleMat)
{
//advance number of randomed part count
if(log.back().PieceCollector.size()>1) log.back().advanceRandomed();
2018-01-07 20:08:50 +01:00
//'set=true' all 4 rotations of pieces in puzzleBox
for(int i=0;i<puzzleMat.p_myBox.size();i++)
if(puzzleMat.p_myBox[i]->GetPartID()==log.back().PieceCollector.begin()->second->GetPartID())
puzzleMat.p_myBox[i]->set=true;
2017-12-03 20:12:32 +01:00
2018-01-07 20:08:50 +01:00
puzzleMat.combinedQualityVector.clear(); //clear data from temp variable
//tell log entry that it is set
log.back().Set();
2018-01-07 20:08:50 +01:00
puzzleMat.setConstraints(log.back().myCoor,log.back().PieceCollector.begin()->second);
cout << "set:" << log.back().myCoor.col << "," << log.back().myCoor.row << endl;
2017-12-03 20:12:32 +01:00
}
bool backtrack(vector<LogEntry>& log, Puzzle& puzzleMat)
{
2018-01-07 20:08:50 +01:00
if(log.empty())
{
cout << "Puzzle not solveable!" << endl;
return false;
}
puzzleMat.combinedQualityVector.clear(); //remove all data from temp quality save
//if more pieces possible, take next piece
if((log.back().PieceCollector.size())>1)
{
puzzleMat.p_myBox.push_back(log.back().PieceCollector.begin()->second);
2017-12-03 20:12:32 +01:00
log.back().PieceCollector.erase(log.back().PieceCollector.begin());
if(log.back().PieceCollector.size()==1)
log.back().decreaseRandomed();
2017-12-03 20:12:32 +01:00
setsolution(log,puzzleMat);
return true;
}
//else remove log element and backtrack once more
else
{
puzzleMat.removeConstrains(log.back().myCoor); //this should remove constraints from all layers
if((log.back().PieceCollector.size()))
2018-01-07 20:08:50 +01:00
for(int i=0;i<puzzleMat.p_myBox.size();i++)
if(puzzleMat.p_myBox[i]->GetPartID()==log.back().PieceCollector.begin()->second->GetPartID())//sets all with partid
puzzleMat.p_myBox[i]->set=false;
log.pop_back();
2018-01-07 20:08:50 +01:00
if(!backtrack(log,puzzleMat))
return false;
return true;
}
}
//this is addon stuff that should later all be extracted into a sererate cpp as it is not core dispatcher functionality
2017-12-22 22:55:55 +01:00
void calculateTrueDestructionPower(vector<LogEntry>& log, Puzzle& puzzleMat, float Layerworth) {
float destructionPower = sqrt(
2018-01-07 20:08:50 +01:00
Layerworth * puzzleMat.dp.m_constraintMatrix[0][0].SpeedTable[log.back().abstractionLevel+1]);
//save destructionArray for when coor is done
if(puzzleMat.tmp_destructionArray.empty())
for(auto it:puzzleMat.dp.m_constraintMatrix[log.back().myCoor.col][log.back().myCoor.row].DestructionArray)
puzzleMat.tmp_destructionArray.emplace_back(it);
puzzleMat.tmp_destructionArray[log.back().abstractionLevel]=destructionPower;
2017-12-22 22:55:55 +01:00
}
// PART RAUER_WEIDINGER
float capLogElements(vector<LogEntry>& log)
2017-12-22 22:55:55 +01:00
{
// Till Now only ground structure -> incorrect variable ans vector names
double limit = 0.6;
double diff = 0;
int id=0;
2017-12-22 22:55:55 +01:00
double maxdiff = 0;
int vectorsizeBefore = 0;
int vectorsizeAfter = 0;
double destroyed = 0; // destroyed parts in %
vectorsizeBefore = log.back().PieceCollector.size();
sort(log.back().PieceCollector.begin(),log.back().PieceCollector.end()); // Sort the vector after probabilities
reverse(log.back().PieceCollector.begin(),log.back().PieceCollector.end());
for(;id<log.back().PieceCollector.size();id++)
{
if(log.back().PieceCollector[id].first < limit)
2017-12-22 22:55:55 +01:00
break;
}
2018-01-07 20:08:50 +01:00
int newid=0;
if(id>0)
newid = --id; //set to the one just over limit
2017-12-22 22:55:55 +01:00
2018-01-07 20:08:50 +01:00
while(id<(log.back().PieceCollector.size()-1)) //find maximum difference in function
2017-12-22 22:55:55 +01:00
{
2018-01-07 20:08:50 +01:00
if(!log.back().PieceCollector[id].first)
break;
diff = log.back().PieceCollector[id].first - log.back().PieceCollector[++id].first;
2017-12-22 22:55:55 +01:00
if(diff > maxdiff)
{
maxdiff = diff;
newid = id;
2017-12-22 22:55:55 +01:00
}
}
cut(log,newid);
2017-12-22 22:55:55 +01:00
vectorsizeAfter = log.back().PieceCollector.size();
destroyed = ((double)vectorsizeBefore - (double)vectorsizeAfter) / (double)vectorsizeBefore;
return (float)sqrt(destroyed*maxdiff);
2017-12-22 22:55:55 +01:00
}
void cut(vector<LogEntry>& log, int& cutID)
{
while(cutID<log.back().PieceCollector.size())
log.back().PieceCollector.erase(log.back().PieceCollector.begin()+cutID);
}
2017-12-22 22:55:55 +01:00
// -------------------- Part David: SetBest and CalculateCombinedQuality --------------------
// pruefen, ob mehr als X combinedQualities ueber dem Grenzwert sind. Wenn nur noch Y Pieces ueber bleiben, dann setBest!
// geeignete Threshold values muessen noch getestet werden
bool SetBestOrMoreLayersArithmetical(vector<LogEntry>& log, qualityVector& cqVector)
2017-12-22 22:55:55 +01:00
{
2018-01-07 20:08:50 +01:00
float threshold, tempBest = 0.0;
unsigned int countHigherThreshold = 0;
2017-12-22 22:55:55 +01:00
if(cqVector.empty())
2017-12-22 22:55:55 +01:00
{
2018-01-07 20:08:50 +01:00
//cerr << "combinedQualityVector is empty." << endl; // should not be empty => backtrack?
return false; // Warning: can only return true or false. What return for error?
2017-12-22 22:55:55 +01:00
}
else
{
2018-01-07 20:08:50 +01:00
switch(log.back().abstractionLevel+1)
{
case 1: threshold = 0.90; break;
case 2: threshold = 0.80; break;
case 3: threshold = 0.75; break;
case 4: threshold = 0.66; break;
case 5: threshold = 0.60; break;
default: threshold = 0.5; break;
}
2018-01-07 20:08:50 +01:00
// check Quality of current Puzzle Piece in combinedQualityVector with Threshold value
for (qualityVector::iterator it = cqVector.begin(); it != cqVector.end(); it++)
if ((cqVector.back().first / log.back().abstractionLevel) >= threshold) // const threshold values
// count how many Pieces are greater than the threshold value
countHigherThreshold++;
else
if ((cqVector.back().first / log.back().abstractionLevel) > tempBest)
tempBest = cqVector.back().first; // could be used, for additional constraints
2017-12-22 22:55:55 +01:00
// return true if only one piece is left
if (1 == countHigherThreshold)
{
return true;
}
else
{
return false;
}
}
2017-12-22 22:55:55 +01:00
}
// jede Quality vom Piececollector zu einer combinedQuality aufsummieren (von jedem bereits verwendetem Layer)
// Achtung: Es muss noch der Mittelwert gebildet werden => SetBestOrMoreLayersArithmetical
void CalculateNewCombinedQuality(vector<LogEntry>& log, qualityVector& qVector, qualityVector& cqVector)
2017-12-22 22:55:55 +01:00
{
bool summarizedVectors = false;
int countSummarizedVectors = 0;
2017-12-22 22:55:55 +01:00
// check if both qualityVectors are not empty
if(qVector.empty())
{
2018-01-07 20:08:50 +01:00
//cerr << "qualityVector is empty." << endl; // should not be empty => backtrack?
return;
}
2018-01-07 20:08:50 +01:00
if(cqVector.empty())
{
2018-01-07 20:08:50 +01:00
//cout << "combinedQualityVector was initialized." << endl; //first layer stuff eh
for(auto it:qVector)
cqVector.emplace_back(it);
2017-12-22 22:55:55 +01:00
}
2018-01-07 20:08:50 +01:00
else
{
for (unsigned int i = 0; i < cqVector.size(); i++) {
for (unsigned int j = 0; j < qVector.size(); j++) {
// search same PuzzlePart of qualityVector and combinedQualityVector
if (&cqVector.at(i).second == &qVector.at(j).second) {
// sum Quality of PieceCollector (qualityVector) to combinedQualityVector
cqVector.at(j).first += qVector.at(i).first;
countSummarizedVectors++;
break; // skip remaining for loop => save time!
}
// remove element at poisition X in combinedQualityVector, because it was not summarized
// inefficient way to delete element X
//cqVector->erase(cqVector->begin()+i);
// efficient way, but no sorted cqVector => wayne //echt? lol
swap(cqVector.at(i), cqVector.back());
cqVector.pop_back();
}
}
// cqVector should have the same size now as newest qVector
if (cqVector.size() != qVector.size()) {
cerr << "Size of combinedQualityVector doenst match with size of qualityVector!" << endl;
cout << "Size of combinedQualityVector: " << cqVector.size() << endl;
cout << "Size of qualityVector: " << qVector.size() << endl;
cout << "Size of countSummarizedVectors: " << countSummarizedVectors << endl;
}
}
}