//Inverse la matrice avec l'affichage en sortie.
public Matrice Inverse(out List <String> display)
{
if (IsInversible())
{
Matrice inverse;
display = new List <string>();
Matrice LPrime = InverseL();
Matrice UPrime = InverseU();
inverse = UPrime.Product(LPrime);
//Permute les colonnes si il y a eu des permutation de lignes.
for (int i = Swaps.GetLength(0) - 1; i >= 0; i--)
{
if (Swaps[i, 0] != -1)
{
inverse.SwapCol(Swaps[i, 0], Swaps[i, 1]);
display.Add("Permutation de colonne : " + (Swaps[i, 0] + 1) + " <->" + (Swaps[i, 1] + 1));
}
}
return(inverse);
}
else
{
throw new Exception("La matrice n'est pas inversible : déterminant = 0");
}
}
public void TimeFromIntervalTest()
{
const string period = "2W";
const double expected = 2 / 52d;
double actual = new Swaps().TimeFromInterval(period);
Assert.AreEqual(expected, actual);
}
public AABB2D(float minX, float minY, float maxX, float maxY)
{
Swaps.minMax(ref minX, ref maxX);
Swaps.minMax(ref minY, ref maxY);
center_.x = (maxX + minX) * 0.5f;
center_.y = (maxY + minY) * 0.5f;
len_.x = (maxX - minX) * 0.5f;
len_.y = (maxY - minY) * 0.5f;
}
protected override void Read(TextReader reader)
{
var lines = ReadLines(reader);
// TODO: could trim the string, should also check for caps, and/or just ToUpper
_workbench = new Workbench(int.Parse(lines[0]), lines[1]);
_swaps = new Swaps(int.Parse(lines[2]),
(from x in lines[3].Split(' ') select int.Parse(x)).ToArray());
}
public void Evaluate(Swaps swaps)
{
foreach (var swap in swaps)
{
var crystals = ElementalCrystals;
if (TrySwap(crystals, swap, out crystals))
ElementalCrystals = crystals;
}
}
//Inverse la matrice.
public void Inverse()
{
if (IsInversible())
{
Matrice LPrime = InverseL();
Matrice UPrime = InverseU();
this.values = (UPrime.Product(LPrime)).values;
//Permute les colonnes si il y a eu des permutation de lignes.
for (int i = 0; i < Swaps.GetLength(0); i++)
{
if (Swaps[i, 0] != -1)
{
SwapCol(Swaps[i, 0], Swaps[i, 1]);
}
}
}
else
{
throw new Exception("La matrice n'est pas inversible : déterminant = 0");
}
}
/// <summary>
/// Сканирует игровое поле на возможные смены камешков
/// </summary>
public void DetectPossibleSwaps()
{
// очищаем список доступных свопов
Swaps.Clear();
for (int row = 0; row < RowsNumber; row++)
{
for (int column = 0; column < ColumnsNumber; column++)
{
// Запоминаем камешек в текущей позиции
Gem gem = GemArray[row, column];
if (gem != null)
{
// Проверка обмена вправо
// Последний ряд не сканируем
if (column < ColumnsNumber - 1)
{
// Запоминаем камешек в следующем столбце
Gem otherGem = GemArray[row, column + 1];
if (otherGem != null)
{
// Меняем местами камешки в массиве
GemArray[row, column] = otherGem;
GemArray[row, column + 1] = gem;
// Если получилась цепь - запоминаем свопы в одну и в другую сторону
if (GetChainAt(row, column).Count >= 3 || GetChainAt(row, column + 1).Count >= 3)
{
Swaps.Add(new Swap(gem, otherGem));
Swaps.Add(new Swap(otherGem, gem));
}
// Возвращаем массив в исходное состояние
GemArray[row, column] = gem;
GemArray[row, column + 1] = otherGem;
}
}
// Проверка обмена вверх, принцип тот же
if (row < RowsNumber - 1)
{
Gem otherGem = GemArray[row + 1, column];
if (otherGem != null)
{
GemArray[row, column] = otherGem;
GemArray[row + 1, column] = gem;
if (GetChainAt(row, column).Count >= 3 || GetChainAt(row + 1, column).Count >= 3)
{
Swaps.Add(new Swap(gem, otherGem));
Swaps.Add(new Swap(otherGem, gem));
}
GemArray[row, column] = gem;
GemArray[row + 1, column] = otherGem;
}
}
}
}
}
// ________ДЛЯ ДЕБАГА_______
foreach (Swap swap in Swaps)
{
Console.Write(swap + "\n");
}
}
public IgnoreRect(Vector2 min, Vector2 max)
{
min_ = min;
max_ = max;
Swaps.minMax(ref min_, ref max_);
}
// 一定距離以上離れてランダムに配置
// 領域の大きさに対してdistが長い場合、指定数が生成される保障は無くなります
// min : 領域の左下座標
// max : 領域の右上座標
// dist: 最短隔離距離
// num : 発生する点の数
static public List <Vector2> distanceBase(Vector2 min, Vector2 max, float dist, int num, List <IgnoreShape> ignoreShapes = null)
{
Swaps.minMax(ref min, ref max);
dist = Mathf.Clamp(dist, 0.00001f, Mathf.Abs(dist));
var list = new List <Vector2>();
if (ignoreShapes == null)
{
ignoreShapes = new List <IgnoreShape>();
}
var wh = max - min; // 幅高
// 領域の幅高どちらも最短距離に満たない場合は1点しか作成できない（とする）
if (max.x - min.x < dist && max.y - min.y < dist)
{
num = 1;
}
// 領域をdist幅以上の区画に分割
int sepNumX = ( int )(wh.x / dist);
int sepNumY = ( int )(wh.y / dist);
sepNumX = sepNumX == 0 ? 1 : sepNumX;
sepNumY = sepNumY == 0 ? 1 : sepNumY;
sepNumX = (sepNumX > 128 ? 128 : sepNumX);
sepNumY = (sepNumY > 128 ? 128 : sepNumY);
int totalCellNum = sepNumX * sepNumY;
float unitX = wh.x / sepNumX;
float unitY = wh.y / sepNumY;
var unit = new Vector2(unitX, unitY);
// 座標格納配列確保
var cells = new Vector2[sepNumX, sepNumY];
var pointExists = new bool[sepNumX, sepNumY];
// セル座標番号をシャッフル
var indices = new List <int>();
ListUtil.numbering(ref indices, sepNumX * sepNumY);
ListUtil.shuffle(ref indices);
var colPoints = new Vector2[9];
int checkNum = 0; // チェックしたセルの数
for (int i = 0; i < num && checkNum < indices.Count; ++i)
{
int e = indices[checkNum];
int x = e % sepNumX;
int y = e / sepNumX;
var cellPos = new Vector2(x * unitX, y * unitY);
// セルが除外範囲内の場合はスキップ
foreach (var ig in ignoreShapes)
{
if (ig.isIgnore(cellPos, cellPos + unit) == true)
{
// 除外
checkNum++;
i--;
}
}
// 周囲8セルの点の座標を取得
int sx = x > 0 ? x - 1 : x;
int ex = x < sepNumX - 1 ? x + 1 : x;
int sy = y > 0 ? y - 1 : y;
int ey = y < sepNumY - 1 ? y + 1 : y;
int cpn = 0; // 相手点の数
for (int oy = sy; oy <= ey; ++oy)
{
for (int ox = sx; ox <= ex; ++ox)
{
if (pointExists[ox, oy] == true)
{
colPoints[cpn] = cells[ox, oy];
cpn++;
}
}
}
const int tryNum = 16; // 試行回数
bool isValid = true;
Vector2 p = Vector2.zero;
for (int t = 0; t < tryNum; ++t)
{
// 設定点に対して周囲8セル内の点との距離をチェック
p = Vector2.one * cellPos + Randoms.Vec2.value(unitX, unitY);
for (int n = 0; n < cpn; ++n)
{
float d = (colPoints[n] - p).magnitude;
if (d < dist)
{
// 不採用
isValid = false;
break;
}
else
{
// 点が除外範囲内の場合は不採用
foreach (var ig in ignoreShapes)
{
if (ig.isIgnore(p) == true)
{
isValid = false;
break;
}
}
if (isValid == false)
{
break;
}
}
}
if (isValid == true)
{
break;
}
}
// 採用点が見つかったら格納
if (isValid == true)
{
list.Add(p);
cells[x, y] = p;
pointExists[x, y] = true;
}
else
{
i--; // numは進めない
}
checkNum++; // チェック数は採用不採用関係なくカウントアップ
}
return(list);
}
public static Vector2Int clamp( Vector2Int r, Vector2Int min, Vector2Int max ) {
Swaps.minMax( ref min, ref max );
return new Vector2Int( Int.clamp( r.x, min.x, max.x ), Int.clamp( r.y, min.y, max.y ) );
}
public static void Main(string[] args)
{
var a = new int[] { 2, 3, 4, 1, 5 };
Console.WriteLine(Swaps.minimumSwaps(a));
}