NuCalc - C#
Bachelorproject on finding the nucleolus solution to distribution-problems.
code:
using System;
using System.IO;
namespace NC
{
public class Program
{
static void Main()
{
Console.SetWindowSize(80, 50);
Console.ForegroundColor = ConsoleColor.White;
NC.Input Input = new NC.Input();
NC.Calc Calc = new NC.Calc();
NC.Pre Pre = new NC.Pre();
char Again;
do
{
Input.MainIn();
Calc.MainCal(Input);
Pre.PreMain(Calc, Input);
Console.Write("\n\n Vil du lave endnu en beregning?\n (Tast 'J' for ja eller 'N' for nej)\n");
Again = Console.ReadLine()[0];
}
while (Again == 'j' || Again == 'J');
}
}
public class Input
{
public string FilePlace;
public string[] FileDataInit = new string[50];
public string[,] FileData = new string[50, 2];
public char GameType;
public char InputType;
public int Players;
public char GenKoa;
public char SpecKoa;
public string Koa;
public string[] KoaSplit = new string[10];
public string KoaSize;
public double[] GameValues = new double[50];
public double[,] GameMat = new double[50, 100];
public int FileLines;
public int i;
public int j;
public int k;
public int KoaPlrs;
public double KoaValue;
public int MatLine;
public double Tolerance;
public void MainIn()
{
Tolerance = 0.001;
Types();
if (InputType == 'm' || GameType == 'M')
{
ManualSize();
SetupMat();
ManualValuesGen();
ManualValuesSpec();
}
if (InputType == 'f' || InputType == 'F')
{
FileIn();
FileSize();
SetupMat();
FileValues();
}
SetupMat2();
Console.Write("\n\n Spillet er opstillet i en matrix:\n");
UdskrivMat(GameMat, (int)Math.Pow(2, Players), Players + 2);
}
public void Types()
{
Console.WriteLine(" Velkommen til NuCalc, til bestemmelse af\n nucleolus-løsningen for kooperative spil.\n");
do
{
Console.WriteLine(" Er spillet et omkostningsspil eller et værdispil?\n (tast 'C' for omkostning eller 'V' for værdi)");
GameType = Console.ReadLine()[0];
}
while (GameType != 'c' && GameType != 'C' && GameType != 'v' && GameType != 'V');
do
{
Console.WriteLine("\n Du har mulighed for at indtaste værdierne for spillet");
Console.WriteLine("\n manuelt eller indlæse dem fra en fil\n (tast 'M' for manuel eller 'F' for fil)");
InputType = Console.ReadLine()[0];
}
while (InputType != 'm' && InputType != 'M' && InputType != 'f' && InputType != 'F');
}
public void FileIn()
{
Console.WriteLine("\n Indtast placeringen på den fil, som skal bruges til\n indlæsning af værdier: ");
FilePlace = Console.ReadLine();
string[] FileDataInit = File.ReadAllLines(FilePlace);
FileLines = FileDataInit.Length;
for (i = 0; i < FileLines; i++)//Der skal være '\t' i hver linje
{
for (j = 0; j < 2; j++)
{
FileData[i, j] = FileDataInit[i].Split('\t')[j];
}
}
}
public void FileSize()
{
Players = int.Parse(FileData[0, 0]);
}
public void FileValues()
{
for (i = 1; i < FileLines; i++)
{
Koa = FileData[i, 0];
KoaSplit = Koa.Split(new Char[] { ',' });
Console.Write("\n Indsætter koalitionsværdi: {");
for (j = 0; j < (KoaSplit.Length - 1); j++)
{
Console.Write(KoaSplit[j] + ",");
}
Console.Write(KoaSplit[KoaSplit.Length - 1] + "}");
KoaValue = double.Parse(FileData[i, 1]);
MatLine = (int)Math.Pow(2, Players) - 1;
for (j = 0; j < (KoaSplit.Length); j++)
{
MatLine -= (int)Math.Pow(2, Players - int.Parse(KoaSplit[j]));
}
GameMat[MatLine, 0] = KoaValue;
}
}
public void ManualSize()
{
do
{
Console.WriteLine(" Hvor mange spillere er der i spillet?\n (antallet skal være mellem 2 og 10)");
Players = int.Parse(Console.ReadLine());
}
while (Players <= 2 && Players >= 10);
}
public void ManualValuesGen()//Indskrivning af koalitionsværdier efter koalitionskombinationerne er opstillet i SetupMat()
{
Console.WriteLine(" Der skal nu angives værdier for de forskellige koalitioner.\n");
Console.WriteLine(" Der skal ialt angives " + (Math.Pow(2, Players)) + " værdier.");
do
{
Console.WriteLine("\n Vil du angive nogle generelle værdier for forskellige ");
Console.WriteLine("\n koalitionsstørrelser? ('J' for ja eller 'N' for nej)");
GenKoa = Console.ReadLine()[0];
}
while (GenKoa != 'j' && GenKoa != 'J' && GenKoa != 'n' && GenKoa != 'N');
while (GenKoa == 'j' || GenKoa == 'J')
{
Console.WriteLine("\n Angiv den koalitionsstørrelse, som du vil angive en værdi for.");
Console.WriteLine("\n (tast '0' for den tomme koalition, '" + Players + "' for den samlede koalition");
Console.WriteLine("\n eller et tal derimellem for en anden koalitionsstørrelse.)");
Console.WriteLine("\n (tast evt. 'A' for at gennemløbe alle koalitionsstørrelser.)");
KoaSize = Console.ReadLine();
if (KoaSize[0] == 'a' || KoaSize[0] == 'A')
{
Console.WriteLine("\n Hvad er vædien for den tomme koalition?");
GameMat[(int)(Math.Pow(2, Players) - 1), 0] = double.Parse(Console.ReadLine());
Console.WriteLine(" \n Hvad er værdien for koalitioner der kun har én spiller?");
KoaValue = double.Parse(Console.ReadLine());
for (j = 0; j < Math.Pow(2, Players); j++)
{
KoaPlrs = 0;
for (k = 1; k <= Players; k++)
{
KoaPlrs += (int)GameMat[j, k];
}
if (KoaPlrs == 1)
{
GameMat[j, 0] = KoaValue;
}
}
for (i = 2; i < Players; i++)
{
Console.WriteLine("\n Hvad er værdien for koalitioner med " + i + " spillere?");
KoaValue = double.Parse(Console.ReadLine());
for (j = 0; j < Math.Pow(2, Players); j++)
{
KoaPlrs = 0;
for (k = 1; k <= Players; k++)
{
KoaPlrs += (int)GameMat[j, k];
}
if (KoaPlrs == i)
{
GameMat[j, 0] = KoaValue;
}
}
}
Console.WriteLine("\n Hvad er værdien for den samlede koalition?");
GameMat[0, 0] = double.Parse(Console.ReadLine());
}
if (KoaSize[0] != 'a' && KoaSize[0] != 'A')
{
if (int.Parse(KoaSize) >= 0 && int.Parse(KoaSize) <= Players)
{
Console.WriteLine("\n Hvad er værdien for koalitioner af den afgivne størrelse?");
KoaValue = double.Parse(Console.ReadLine());
for (j = 0; j < Math.Pow(2, Players); j++)
{
KoaPlrs = 0;
for (k = 1; k <= Players; k++)
{
KoaPlrs += (int)GameMat[j, k];
}
if (KoaPlrs == int.Parse(KoaSize))
{
GameMat[j, 0] = KoaValue;
}
}
}
}
do
{
Console.WriteLine("\n Vil du angive flere generelle koalitionsværdier?");
GenKoa = Console.ReadLine()[0];
}
while (GenKoa != 'j' && GenKoa != 'J' && GenKoa != 'n' && GenKoa != 'N');
}
}
public void ManualValuesSpec()
{
do
{
Console.WriteLine("\n Vil du angive værdier for nogle specifikke koalitioner?\n ('J' for ja eller 'N' for nej)");
SpecKoa = Console.ReadLine()[0];
}
while (SpecKoa != 'j' && SpecKoa != 'J' && SpecKoa != 'n' && SpecKoa != 'N');
while (SpecKoa == 'j' || SpecKoa == 'J')
{
Console.WriteLine("\n Hvilken koaltion vil du angive en værdi for?");
Console.WriteLine("\n (indtast numrene på de spillere der indgår i\n koalitionen adskilt af kommaer)");
Koa = Console.ReadLine();
KoaSplit = Koa.Split(new Char[] { ',' });
Console.Write("\n Hvad er værdien for koalitionen der indeholder spillerne {");
for (i = 0; i < (KoaSplit.Length - 1); i++)
{
Console.Write(KoaSplit[i] + ",");
}
Console.WriteLine(KoaSplit[(KoaSplit.Length - 1)] + "}?");
KoaValue = double.Parse(Console.ReadLine());
MatLine = (int)Math.Pow(2, Players) - 1;
for (i = 0; i < (KoaSplit.Length); i++)
{
MatLine -= (int)Math.Pow(2, Players - int.Parse(KoaSplit[i]));
}
GameMat[MatLine, 0] = KoaValue;
do
{
Console.WriteLine("\n Vil du angive flere specifikke koalitionsværdier?");
Console.WriteLine("\n ('J' for ja eller 'N' for nej)");
SpecKoa = Console.ReadLine()[0];
}
while (SpecKoa != 'j' && SpecKoa != 'J' && SpecKoa != 'n' && SpecKoa != 'N');
}
}
public void SetupMat()//Opstiller alle koalitionskombinationer inkl. fuld og tom
{
for (i = 0; i < Math.Pow(2, Players); i++)
{
for (j = 0; j < (2 * (Math.Pow(2, Players)) + Players); j++)
{
GameMat[i, j] = 0;
}
}
for (i = 1; i <= Players; i++)
{
for (j = 0; j < Math.Pow(2, (i - 1)); j++)
{
for (k = 0; k < Math.Pow(2, (Players - i)); k++)
{
GameMat[(int)Math.Pow(2, (Players - (i - 1))) * j + k, i] = 1;
}
}
}
for (i = 1; i < Math.Pow(2, Players) - 1; i++)
{
if (GameType == 'v' || GameType == 'V')
{
GameMat[i, (Players + 1)] = 1;
GameMat[i, (Players + 1 + i)] = -1;
GameMat[i, (Players + (int)Math.Pow(2, Players) + i - 1)] = 1;
}
if (GameType == 'c' || GameType == 'C')
{
GameMat[i, (Players + 1)] = -1;
GameMat[i, (Players + 1 + i)] = 1;
}
}
GameMat[(int)Math.Pow(2, Players) - 1, (Players + 2 * (int)Math.Pow(2, Players) - 2)] = 1;
}
public void SetupMat2()//Sletter den tomme koalition og flytter den fulde koalition
{
for (i = 0; i < Math.Pow(2, Players) - 1; i++)
{
GameMat[i, 0] -= GameMat[(int)Math.Pow(2, Players) - 1, 0];//Trækker værdien af den tomme koalition fra for at standardisere
}
for (i = 0; i <= Players; i++)
{
GameMat[(int)Math.Pow(2, Players) - 1, i] = GameMat[0, i];
GameMat[0, i] = 0;
}
}
public void UdskrivMat(double[,] Mat, int Height, int Width)
{
Console.ForegroundColor = ConsoleColor.Yellow;
for (i = 0; i < Height; i++)
{
for (j = 0; j < Width; j++)
{
Console.Write(" " + Mat[i, j] + " ");
}
Console.Write("\n");
}
Console.ReadLine();
Console.ForegroundColor = ConsoleColor.White;
}
}
public class Calc
{
public double[] Shapley = new double[10];
public double Alpha;
public bool CoreEmpty;
public bool AlphaMax;
public int i;
public int j;
public int Rows;
public int Columns;
public int Row;
public int Column;
int HasAlpha;
public double[,] GameMatHlp = new double[50, 100];
public double[,] GameMatMain = new double[50, 100];
public double Tolerance;
public double Max;
public int[] RowNumbers = new int[50];
public int[] ColumnNumbers = new int[50];
public int[] SlacksRows = new int[50];
public int Slacks;
public int NewSlacks;
public double[] Nucleolus = new double[11];
public void MainCal(Input Input)
{
FindShapley(Input);
AlphaMax = false;
CoreEmpty = false;
Slacks = 0;
Rows = (int)Math.Pow(2, Input.Players) - 1;
Columns = (int)(2 * (Math.Pow(2, Input.Players)) + Input.Players - 2);
if (Input.GameType == 'v' || Input.GameType == 'V')
{
SimpHlpSetupV(Input);
}
if (Input.GameType == 'c' || Input.GameType == 'C')
{
SimpHlpSetupC(Input);
}
Tolerance = Input.Tolerance;
if (AlphaMax == false)
{
do
{
Simplex(GameMatHlp);
FindBasis(GameMatHlp);
//FindNucleolus(Input, GameMatHlp);
HasAlpha = 0;
for (i = 0; i < Rows; i++)
{
if (ColumnNumbers[i] == Input.Players + 1)
{ HasAlpha = RowNumbers[i]; }
}
NewSlacks = 0;
Columns = (int)(Math.Pow(2, Input.Players) + Input.Players - 1);
if (HasAlpha != 0)
{
SimpMainSetup(Input);
Simplex(GameMatMain);
FindBasis(GameMatMain);
HasAlpha = 0;
Alpha = 0;
for (i = 0; i < Rows; i++)
{
if (ColumnNumbers[i] == Input.Players + 1)
{ HasAlpha = RowNumbers[i]; }
}
if (HasAlpha != 0)
{
Alpha = GameMatMain[HasAlpha, 0];
}
if (Alpha != 0)
{
FindSlacks(Input, GameMatMain);
UpdateMain(Input, GameMatMain);
FindNucleolus(Input, GameMatMain);
}
}
if (NewSlacks != 0)
{
Columns = (int)(2 * (Math.Pow(2, Input.Players)) + Input.Players - 2);
if (Input.GameType == 'v' || Input.GameType == 'V')
{
SimpHlpSetupV(Input);
}
if (Input.GameType == 'c' || Input.GameType == 'C')
{
SimpHlpSetupC(Input);
}
SimpHlpSetupSlack(Input);
CoreEmpty = true;
}
if (NewSlacks == 0 && Math.Abs(Nucleolus[Input.Players]) <= Tolerance)
{
AlphaMax = true;//Ændre til maksimering; skift fortegn for alpha og nyt hovedproblem
}
}
while (NewSlacks != 0);
}
if (AlphaMax == true)
{
for (i = 1; i < Rows; i++)
{
Input.GameMat[i, (Input.Players + 1)] *= -1;
}
Rows = (int)Math.Pow(2, Input.Players) - 1;
Columns = (int)(2 * (Math.Pow(2, Input.Players)) + Input.Players - 2);
if (Input.GameType == 'v' || Input.GameType == 'V')
{
SimpHlpSetupV(Input);
}
if (Input.GameType == 'c' || Input.GameType == 'C')
{
SimpHlpSetupC(Input);
}
do
{
Simplex(GameMatHlp);
//FindNucleolus(Input, GameMatHlp);
Columns = (int)(Math.Pow(2, Input.Players) + Input.Players - 1);
SimpMainSetup2(Input);//NY FUNKTION - simpel: udskift kriterie med nuller og et 1-tal
Simplex(GameMatMain);
FindBasis(GameMatMain);
NewSlacks = 0;
FindSlacks(Input, GameMatMain);
UpdateMain(Input, GameMatMain);
FindNucleolus(Input, GameMatMain);
if (NewSlacks != 0)
{
Columns = (int)(2 * (Math.Pow(2, Input.Players)) + Input.Players - 2);
if (Input.GameType == 'v' || Input.GameType == 'V')
{
SimpHlpSetupV(Input);
}
if (Input.GameType == 'c' || Input.GameType == 'C')
{
SimpHlpSetupC(Input);
}
SimpHlpSetupSlack(Input);
}
}
while (NewSlacks != 0);
}
}
public static Int64 Fact(int t)
{
long factorial = t;
while (t > 1)
factorial *= --t;
return factorial;
}
public void FindShapley(Input Input)//Skal køres først, da matricen nedbrydes
{
for (int j = 1; j <= Input.Players; j++)
{
Shapley[j-1] = 0;
for (int i = 1; i < Math.Pow(2,Input.Players); i++)
{
if (Input.GameMat[i, j] == 1)
{
int ks = 0;
for (int k = 1; k <= Input.Players; k++)
{
ks += (int)Input.GameMat[i, k];
}
if(ks==1)
{
Shapley[j-1]+=(Input.GameMat[i,0]*Fact(Input.Players-1));
}
if (ks > 1 && ks < Input.Players)
{
Shapley[j-1]+=((Input.GameMat[i,0]-Input.GameMat[i + (int)Math.Pow(2, (Input.Players-j)),0])*Fact(ks-1)*Fact(Input.Players-ks));
}
if (ks == Input.Players)
{
Shapley[j-1] += ((Input.GameMat[i, 0] - Input.GameMat[(int)Math.Pow(2, (Input.Players-j)), 0]) * Fact(Input.Players - 1));
}
}
}
Shapley[j-1] /= Fact(Input.Players);
}
}
public void UpdateMain(Input Input, double[,] Mat)
{
for (i = 0; i < NewSlacks; i++)
{
if (Input.GameType == 'v' || Input.GameType == 'V')
{
Input.GameMat[SlacksRows[Slacks - NewSlacks + i], 0] -= Mat[0, 0];
}
if (Input.GameType == 'c' || Input.GameType == 'C')
{
Input.GameMat[SlacksRows[Slacks - NewSlacks + i], 0] += Mat[0, 0];
}
Input.GameMat[SlacksRows[Slacks - NewSlacks + i], Input.Players + 1] = 0;
}
}
public void SimpMainSetup(Input Input)
{
for (i = 0; i <= Rows; i++)
{
for (j = 0; j <= Columns; j++)
{
GameMatMain[i, j] = GameMatHlp[i, j];
}
}
for (i = 0; i < Columns; i++)
{
GameMatMain[0, i] = GameMatHlp[HasAlpha, i];
}
GameMatMain[0, Input.Players + 1] = 0;
}
public void SimpMainSetup2(Input Input)
{
for (i = 0; i <= Rows; i++)
{
for (j = 0; j <= Columns; j++)
{
GameMatMain[i, j] = GameMatHlp[i, j];
}
}
for (i = 0; i < Columns; i++)
{
GameMatMain[0, i] = 0;
}
GameMatMain[0, Input.Players + 1] = 1;
}
public void SimpHlpSetupV(Input Input)
{
for (i = 0; i <= Rows; i++)
{
for (j = 0; j <= Columns; j++)
{
GameMatHlp[i, j] = Input.GameMat[i, j];
}
}
for (i = 0; i <= (Input.Players + 1); i++)
{
for (j = 1; j < Math.Pow(2, Input.Players); j++)
{
GameMatHlp[0, i] += Input.GameMat[j, i];
}
}
for (i = Input.Players + 2; i < (Math.Pow(2, Input.Players) + Input.Players); i++)
{
{ GameMatHlp[0, i] = -1; }
}
}
public void SimpHlpSetupC(Input Input)
{
for (i = 0; i <= Rows; i++)
{
for (j = 0; j <= Columns; j++)
{
GameMatHlp[i, j] = Input.GameMat[i, j];
}
}
for (i = 0; i <= (Input.Players + 1); i++)
{
GameMatHlp[0, i] += Input.GameMat[(int)Math.Pow(2, Input.Players) - 1, i];
}
}
public void SimpHlpSetupSlack(Input Input)
{
for (i = 0; i < Slacks; i++)
{
GameMatHlp[SlacksRows[i], SlacksRows[i] + Input.Players + 1] = 0;//Fjerner slack
GameMatHlp[0, SlacksRows[i] + Input.Players + 1] = 0;//Fjerner slack fra kriterie
GameMatHlp[SlacksRows[i], Input.Players + (int)Math.Pow(2, Input.Players) - 1 + SlacksRows[i]] = 1;//Tilføjer kunstig slack
if (Input.GameType == 'c' || Input.GameType == 'C')
{
for (j = 0; j <= Input.Players; j++)
{
GameMatHlp[0, j] += GameMatHlp[SlacksRows[i], j];//tilføjer for omkostningsspil rækker uden slack/alpha
}
}
}
}
public void Simplex(double[,] Mat) //Udfører de sædvanelige simplex udregninger således at vi en ny basis
{
int l = 0;
do //Løkken køre sålænge den største værdi i c-vektoren er større end en værdi lidt større end 0
{
FindBig(Mat); //Finder største værdi i c-vektoren og gemmer hvilken søjle den er placeret i
if (Max > Tolerance) //Tjekker om største værdi er større end 0
{
l++;
FindPivotElement(Mat); //Finder Pivotelementet og gemmer rækken den er placeret i
if (Max > Tolerance)
{
NewTable(Row, Column, Mat); //Laver rækker operationer således at der står nuller over og under pivotelementet, som sættes til 1
}
}
}
while (Max > Tolerance);
}
public void FindBig(double[,] Mat) //Finder største værdi i c-vektoren og gemmer hvilken søjle den er placeret i
{
Column = 1;
Max = Mat[0, 1]; //Start værdien
for (int j = 2; j <= Columns; j++) //Løkken køre igennem hvert tal i rækken, på nær det første
{
if (Max < Mat[0, j]) //Tjekker om værdien på j'te plads i rækken er større end den tidligere gemte største værdi
{
Max = Mat[0, j]; //Gemmer værdien på j'te plads i rækken som "største"
Column = j; //Gemmer søjlenummeret, som bruges i "FindPivotelement"
}
}
}
public void FindPivotElement(double[,] Mat) //Finder Pivotelementet og gemmer rækken den er placeret i
{
double Frac, PivotFrac;
int a = 0;
PivotFrac = 0;
Row = -1; //Start værdi
if (Mat[Rows, Column] > Tolerance)//Tester nederste række før alle andre
{
PivotFrac = Mat[Rows, 0] / Mat[Rows, Column];
Row = Rows;
a = 1;
}
for (int i = 1; i <= Rows; i++) //Løkken køre igennem hvert tal i søjlen, på nær det første
{//Tilføj evt Mat[i,0]/Mat[i,Column]>Tolerance for test
if (Mat[i, Column] > Tolerance) //Tjekker om værdien i søjlen, fundet i "FindStørste", er større end en værdi lidt større end 0
{
Frac = Mat[i, 0] / Mat[i, Column]; //Finder forholdet mellem i'te række i søjle og i'te række i søjlen vi ser på
if (a == 0) //Første gange gemmes brøken i pivotbrøken
{
PivotFrac = Frac;
Row = i;
}
if (Frac < PivotFrac) //Hvis brøken er mindre end den tidligere gemt pivotbrøk, gemmes den i pivotbrøk
{
PivotFrac = Frac;
Row = i;
}
a++;
}
}
if (Row == -1) //Hvis række nummeret ikke er blevet ændret, så skriver programmet ud at der ikke er fundet et pivotelemenent og sætter "største" = -1 således at løkken stopper i "SimplexUdregninger"
{
//Console.ForegroundColor = ConsoleColor.Magenta;
//Console.WriteLine("\n\n Programmet har desværre ikke kunnet finde\n et pivotelement");
Row = 1;
Max = -1;
}
}
public void NewTable(int Row, int Column, double[,] Mat) //Laver rækker operationer således at der står nuller over og under pivotelementet, som sættes til 1
{ //Som en af de to eneste funktion bruger den ikke kun "public" variabler, eftersom den bruges mere end én gang
DivPivotRow(Row, Column, Mat); //Dividerer rækken med pivotelementet igennem med pivotelementet
for (int i = 0; i <= Rows; i++) //De to for-løkker køre igennem alle pladser i tableauet, på nær alle pladserne i den række hvor pivotelementet er placeret
{
if (i != Row)
{
double PivotColumnElement = Mat[i, Column]; //Gemmer i'te plads i søjlen i "pivotsøjleelement"
for (int j = 0; j <= Columns; j++)
{
Mat[i, j] -= PivotColumnElement * Mat[Row, j]; //Laver rækker operationer således at der opstår nuller over og under pivotelementet
}
}
}
}
public void DivPivotRow(int Row, int Column, double[,] Mat) //Dividerer rækken med pivotelementet igennem med pivotelementet
{ //Som en af de to eneste funktion bruger den ikke kun "public" variabler, eftersom den bruges mere end én gang
double PivotElement;
PivotElement = Mat[Row, Column]; //Gemmer pivotelementets værdi i "pivotelement"
for (int j = 0; j <= Columns; j++) //Løkken dividere igennem alle pladser i rækken hvor pivotelementet er placeret med pivotelementet selv
{
Mat[Row, j] /= PivotElement;
}
}
public void FindBasis(double[,] Mat) //Finder basissøjler og gemmer dem i TableauMatricen
{
int l = 0;
int m = 0;
for (int i = 0; i < Rows; i++)//Starter med at slette tidligere fundne basissøjler
{
RowNumbers[i] = 0;
ColumnNumbers[i] = 0;
}
for (int i = 1; i <= Rows; i++) //Løkkerne køre igennem alle rækker og søjler på nær den første
{
for (int j = 1; j <= Columns; j++)
{
if (Math.Abs(Mat[i, j] - 1) < Tolerance) //Tjekker om værdien på i,j'te plads er nær 1
{
m = 0;
for (int n = 0; n <= Rows; n++) //Løkken køre igennem alle pladser for søjle j
{
if (Math.Abs(Mat[n, j]) < Tolerance) //Tjekker om værdien på i,j'te plads er nær 0
{
m++; //Tæller antallet af nuller i søjle j
}
}
if (m == Rows) //Tjekker om der er det passende antal nuller i søjle j
{
RowNumbers[l] = i; //Gemmer placeringen af 1-tallet
ColumnNumbers[l] = j;
l++; //Tæller basis søjler
i++; //Går til næste række
j = -1; //Starter forfra i rækken
}
}
}
}
/*if (l != Rows) //Hvis antallet af basis søjler er forskellig fra det antal der forventes, skriver programmet at vi ikke er i basis
{
Console.ForegroundColor = ConsoleColor.Magenta;
Console.WriteLine("\n\n Der er ikke tilstrækkelige basissøjler til\n at udføre Simplex");
}*/
}
public void FindSlacks(Input Input, double[,] Mat)
{
NewSlacks = 0;
for (i = Input.Players + 2; i <= Columns; i++)
{
if (Math.Abs(Mat[0, i]) > Tolerance)
{
SlacksRows[Slacks] = i - (Input.Players + 1);
Slacks++;
NewSlacks++;
}
}
}
public void FindNucleolus(Input Input, double[,] Mat)
{
for (i = 0; i <= Input.Players; i++)
{
for (j = 0; j < Rows; j++)
{
if (ColumnNumbers[j] == (i + 1))
{
Nucleolus[i] = Mat[RowNumbers[j], 0];
}
}
} //IGNORER RESULTAT OG GÅ TIL MAX VED ALPHA=0!!!!!
}
}
public class Pre
{
public void PreMain(Calc Calc, Input Input)
{
PreCore(Calc);
PreNuc(Calc, Input);
PreShapley(Calc, Input);
Console.ReadLine();
}
public void PreNuc(Calc Calc, Input Input)
{
Console.Write("\n\n Nucleolus blev beregnet til: ");
Console.ForegroundColor=ConsoleColor.Green;
Console.Write("\n { ");
for (int i = 0; i < (Input.Players - 1); i++)
{
Console.Write(Math.Round(Calc.Nucleolus[i],3) + " ; ");
}
Console.Write(Math.Round(Calc.Nucleolus[(Input.Players - 1)],3) + " }");
Console.ForegroundColor = ConsoleColor.White;
}
public void PreCore(Calc Calc)
{
Console.Write("\n\n Kernen for spillet er ");
Console.ForegroundColor = ConsoleColor.Magenta;
if (Calc.CoreEmpty == false)
{ Console.Write("IKKE "); }
Console.Write("TOM.");
Console.ForegroundColor = ConsoleColor.White;
}
public void PreShapley(Calc Calc, Input Input)
{
Console.Write("\n\n Shapley-værdien er beregnet til: ");
Console.ForegroundColor = ConsoleColor.Green;
Console.Write("\n { ");
for (int i = 0; i < (Input.Players - 1); i++)
{
Console.Write(Math.Round(Calc.Shapley[i], 3) + " ; ");
}
Console.Write(Math.Round(Calc.Shapley[(Input.Players - 1)],3) + " }");
Console.ForegroundColor = ConsoleColor.White;
}
}
}
