結果
| 問題 |
No.453 製薬会社
|
| ユーザー |
 yupiteru_kun
|
| 提出日時 | 2023-07-17 02:39:43 |
| 言語 | C# (.NET 8.0.404) |
| 結果 |
AC
|
| 実行時間 | 55 ms / 2,000 ms |
| コード長 | 28,080 bytes |
| コンパイル時間 | 11,357 ms |
| コンパイル使用メモリ | 169,660 KB |
| 実行使用メモリ | 189,684 KB |
| 最終ジャッジ日時 | 2024-09-17 21:19:39 |
| 合計ジャッジ時間 | 10,145 ms |
|
ジャッジサーバーID （参考情報） |
judge2 / judge4 |
（要ログイン）
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 9 |
コンパイルメッセージ
復元対象のプロジェクトを決定しています... /home/judge/data/code/main.csproj を復元しました (96 ms)。 MSBuild のバージョン 17.9.6+a4ecab324 (.NET) main -> /home/judge/data/code/bin/Release/net8.0/main.dll main -> /home/judge/data/code/bin/Release/net8.0/publish/
ソースコード
using System;
using System.Collections.Generic;
using System.IO;
using System.Linq;
using static System.Math;
using System.Text;
using System.Threading;
using System.Globalization;
using System.Runtime.CompilerServices;
using System.Runtime.InteropServices;
using Library;
namespace Program
{
public static class ProblemA
{
static bool SAIKI = false;
static public int numberOfRandomCases = 0;
static public void MakeTestCase(List<string> _input, List<string> _output, ref Func<string[], bool> _outputChecker)
{
}
static public void Solve()
{
var C = NN;
var D = NN;
var simplex = new LIB_LPSimplex();
simplex.AddSubject(new LIB_Fraction(3, 4), new LIB_Fraction(2, 7)).LessThan(C);
simplex.AddSubject(new LIB_Fraction(1, 4), new LIB_Fraction(5, 7)).LessThan(D);
Console.WriteLine((double)simplex.Maximize(1000, 2000));
}
class Printer : StreamWriter
{
public override IFormatProvider FormatProvider { get { return CultureInfo.InvariantCulture; } }
public Printer(Stream stream) : base(stream, new UTF8Encoding(false, true)) { base.AutoFlush = false; }
public Printer(Stream stream, Encoding encoding) : base(stream, encoding) { base.AutoFlush = false; }
}
static LIB_FastIO fastio = new LIB_FastIODebug();
static string[] args;
static public void Main(string[] args_t) { args = args_t; if (args_t.Length == 0) { fastio = new LIB_FastIO(); Console.SetOut(new Printer(Console.OpenStandardOutput())); } if (SAIKI) { var t = new Thread(Solve, 134217728); t.Start(); t.Join(); } else Solve(); Console.Out.Flush(); }
static long NN => fastio.Long();
static double ND => fastio.Double();
static string NS => fastio.Scan();
static long[] NNList(long N) => Repeat(0, N).Select(_ => NN).ToArray();
static double[] NDList(long N) => Repeat(0, N).Select(_ => ND).ToArray();
static string[] NSList(long N) => Repeat(0, N).Select(_ => NS).ToArray();
static long Count<T>(this IEnumerable<T> x, Func<T, bool> pred) => Enumerable.Count(x, pred);
static IEnumerable<T> Repeat<T>(T v, long n) => Enumerable.Repeat<T>(v, (int)n);
static IEnumerable<int> Range(long s, long c) => Enumerable.Range((int)s, (int)c);
static IOrderedEnumerable<T> OrderByRand<T>(this IEnumerable<T> x) => Enumerable.OrderBy(x, _ => xorshift);
static IOrderedEnumerable<T> OrderBy<T>(this IEnumerable<T> x) => Enumerable.OrderBy(x.OrderByRand(), e => e);
static IOrderedEnumerable<T1> OrderBy<T1, T2>(this IEnumerable<T1> x, Func<T1, T2> selector) => Enumerable.OrderBy(x.OrderByRand(), selector);
static IOrderedEnumerable<T> OrderByDescending<T>(this IEnumerable<T> x) => Enumerable.OrderByDescending(x.OrderByRand(), e => e);
static IOrderedEnumerable<T1> OrderByDescending<T1, T2>(this IEnumerable<T1> x, Func<T1, T2> selector) => Enumerable.OrderByDescending(x.OrderByRand(), selector);
static IOrderedEnumerable<string> OrderBy(this IEnumerable<string> x) => x.OrderByRand().OrderBy(e => e, StringComparer.OrdinalIgnoreCase);
static IOrderedEnumerable<T> OrderBy<T>(this IEnumerable<T> x, Func<T, string> selector) => x.OrderByRand().OrderBy(selector, StringComparer.OrdinalIgnoreCase);
static IOrderedEnumerable<string> OrderByDescending(this IEnumerable<string> x) => x.OrderByRand().OrderByDescending(e => e, StringComparer.OrdinalIgnoreCase);
static IOrderedEnumerable<T> OrderByDescending<T>(this IEnumerable<T> x, Func<T, string> selector) => x.OrderByRand().OrderByDescending(selector, StringComparer.OrdinalIgnoreCase);
static string Join<T>(this IEnumerable<T> x, string separator = "") => string.Join(separator, x);
static uint xorshift { get { _xsi.MoveNext(); return _xsi.Current; } }
static IEnumerator<uint> _xsi = _xsc();
static IEnumerator<uint> _xsc() { uint x = 123456789, y = 362436069, z = 521288629, w = (uint)(DateTime.Now.Ticks & 0xffffffff); while (true) { var t = x ^ (x << 11); x = y; y = z; z = w; w = (w ^ (w >> 19)) ^ (t ^ (t >> 8)); yield return w; } }
static bool Chmax<T>(this ref T lhs, T rhs) where T : struct, IComparable<T> { if (lhs.CompareTo(rhs) < 0) { lhs = rhs; return true; } return false; }
static bool Chmin<T>(this ref T lhs, T rhs) where T : struct, IComparable<T> { if (lhs.CompareTo(rhs) > 0) { lhs = rhs; return true; } return false; }
static void Fill<T>(this T[] array, T value) => array.AsSpan().Fill(value);
static void Fill<T>(this T[,] array, T value) => MemoryMarshal.CreateSpan(ref array[0, 0], array.Length).Fill(value);
static void Fill<T>(this T[,,] array, T value) => MemoryMarshal.CreateSpan(ref array[0, 0, 0], array.Length).Fill(value);
static void Fill<T>(this T[,,,] array, T value) => MemoryMarshal.CreateSpan(ref array[0, 0, 0, 0], array.Length).Fill(value);
}
}
namespace Library {
struct LIB_Fraction : IEquatable<LIB_Fraction>, IEquatable<long>, IComparable<LIB_Fraction>, IComparable<long>
{
public long Numerator { get; private set; }
public long Denominator { get; private set; }
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public LIB_Fraction(long nume = 0, long deno = 1)
{
Numerator = nume;
Denominator = deno;
Yakubun();
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
void Yakubun()
{
if (Denominator < 0) { Numerator *= -1; Denominator *= -1; }
var gcd = LIB_Math.GCD(Abs(Numerator), Denominator);
if (gcd == 0) { Numerator = 0; Denominator = 1; }
else { Numerator /= gcd; Denominator /= gcd; }
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public implicit operator LIB_Fraction(long x) => new LIB_Fraction(x);
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public explicit operator double(LIB_Fraction x) => (double)x.Numerator / x.Denominator;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public void Add(LIB_Fraction x)
{
Numerator = Numerator * x.Denominator + x.Numerator * Denominator;
Denominator *= x.Denominator;
Yakubun();
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public void Sub(LIB_Fraction x)
{
Numerator = Numerator * x.Denominator - x.Numerator * Denominator;
Denominator *= x.Denominator;
Yakubun();
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public void Mul(LIB_Fraction x)
{
Numerator *= x.Numerator;
Denominator *= x.Denominator;
Yakubun();
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public void Div(LIB_Fraction x)
{
Numerator *= x.Denominator;
Denominator *= x.Numerator;
Yakubun();
}
public LIB_Fraction Abs
{
[MethodImpl(MethodImplOptions.AggressiveInlining)]
get => new LIB_Fraction(Abs(Numerator), Denominator);
private set { }
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public LIB_Fraction operator +(LIB_Fraction x, LIB_Fraction y) { x.Add(y); return x; }
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public LIB_Fraction operator -(LIB_Fraction x, LIB_Fraction y) { x.Sub(y); return x; }
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public LIB_Fraction operator *(LIB_Fraction x, LIB_Fraction y) { x.Mul(y); return x; }
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public LIB_Fraction operator /(LIB_Fraction x, LIB_Fraction y) { x.Div(y); return x; }
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public bool operator <(LIB_Fraction x, LIB_Fraction y) => x.Denominator == y.Denominator ? x.Numerator < y.Numerator : x.Numerator * y.Denominator < y.Numerator * x.Denominator;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public bool operator >(LIB_Fraction x, LIB_Fraction y) => x.Denominator == y.Denominator ? x.Numerator > y.Numerator : x.Numerator * y.Denominator > y.Numerator * x.Denominator;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public bool operator <=(LIB_Fraction x, LIB_Fraction y) => x.Denominator == y.Denominator ? x.Numerator <= y.Numerator : x.Numerator * y.Denominator <= y.Numerator * x.Denominator;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public bool operator >=(LIB_Fraction x, LIB_Fraction y) => x.Denominator == y.Denominator ? x.Numerator >= y.Numerator : x.Numerator * y.Denominator >= y.Numerator * x.Denominator;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public bool operator ==(LIB_Fraction x, LIB_Fraction y) => x.Numerator == y.Numerator && x.Denominator == y.Denominator;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public bool operator !=(LIB_Fraction x, LIB_Fraction y) => x.Numerator != y.Numerator || x.Denominator != y.Denominator;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public int CompareTo(LIB_Fraction x)
{
if (Denominator == x.Denominator) return Numerator.CompareTo(x.Numerator);
var left = Numerator * x.Denominator;
var right = x.Numerator * Denominator;
if (left < right) return -1;
if (left > right) return 1;
return 0;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public int CompareTo(long x)
{
var right = x * Denominator;
if (Numerator < right) return -1;
if (Numerator > right) return 1;
return 0;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public bool Equals(LIB_Fraction x) => x.Numerator == Numerator && x.Denominator == Denominator;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public bool Equals(long x) => Numerator == x && Denominator == 1;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public override bool Equals(object x) => x == null ? false : Equals((LIB_Fraction)x);
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public override int GetHashCode()
{
unchecked
{
return (Numerator.GetHashCode() << 5 + Numerator.GetHashCode()) ^ Denominator.GetHashCode();
}
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public override string ToString() => Denominator == 1 ? Numerator.ToString() : $"{Numerator}/{Denominator}";
}
class LIB_Math
{
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public IEnumerable<long> Primes(long x)
{
if (x < 2) yield break;
yield return 2;
var halfx = x / 2;
var table = new bool[halfx + 1];
var max = (long)(Math.Sqrt(x) / 2);
for (long i = 1; i <= max; ++i)
{
if (table[i]) continue;
var add = 2 * i + 1;
yield return add;
for (long j = 2 * i * (i + 1); j <= halfx; j += add)
table[j] = true;
}
for (long i = max + 1; i <= halfx; ++i)
if (!table[i] && 2 * i + 1 <= x)
yield return 2 * i + 1;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public IEnumerable<long> Factors(long x)
{
if (x < 2) yield break;
while (x % 2 == 0)
{
x /= 2; yield return 2;
}
var max = (long)Math.Sqrt(x);
for (long i = 3; i <= max; i += 2)
{
while (x % i == 0)
{
x /= i; yield return i;
}
}
if (x != 1) yield return x;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public IEnumerable<long> Divisor(long x)
{
if (x < 1) yield break;
var max = (long)Math.Sqrt(x);
for (long i = 1; i <= max; ++i)
{
if (x % i != 0) continue;
yield return i;
if (i != x / i) yield return x / i;
}
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public long GCD(long a, long b)
{
while (b > 0)
{
var tmp = b;
b = a % b;
a = tmp;
}
return a;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public long LCM(long a, long b) => a / GCD(a, b) * b;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public long Pow(long x, long y)
{
long a = 1;
while (y != 0)
{
if ((y & 1) == 1) a *= x;
if (x < long.MaxValue / x) x *= x;
y >>= 1;
}
return a;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public decimal Sqrt(decimal x)
{
decimal prev, cur = (decimal)Math.Sqrt((double)x);
do
{
prev = cur;
if (prev == 0) return 0;
cur = (prev + x / prev) / 2;
} while (cur != prev);
return cur;
}
static List<long> _fact = new List<long>() { 1 };
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static void Build(long n)
{
if (n >= _fact.Count)
for (int i = _fact.Count; i <= n; ++i) _fact.Add(_fact[i - 1] * i);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public long Comb(long n, long k)
{
Build(n);
if (n == 0 && k == 0) return 1;
if (n < k || n < 0) return 0;
return _fact[(int)n] / _fact[(int)(n - k)] / _fact[(int)k];
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public long Perm(long n, long k)
{
Build(n);
if (n == 0 && k == 0) return 1;
if (n < k || n < 0) return 0;
return _fact[(int)n] / _fact[(int)(n - k)];
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public IEnumerable<List<int>> MakePermutation(long n, bool zeroIndexed = true)
{
if (n <= 0) throw new Exception();
var c = new int[n];
var a = new int[n];
if (!zeroIndexed) a[0] = 1;
for (var i = 1; i < n; i++) a[i] = a[i - 1] + 1;
yield return new List<int>(a);
for (var i = 0; i < n;)
{
if (c[i] < i)
{
if (i % 2 == 0)
{
var t = a[0]; a[0] = a[i]; a[i] = t;
}
else
{
var t = a[c[i]]; a[c[i]] = a[i]; a[i] = t;
}
yield return new List<int>(a);
++c[i];
i = 0;
}
else
{
c[i] = 0;
++i;
}
}
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public bool NextPermutation(long[] ary)
{
var n = ary.Length;
var i = n - 1;
while (i - 1 >= 0 && ary[i - 1] > ary[i]) --i;
if (i == 0) return false;
var j = i;
while (j + 1 < n && ary[i - 1] < ary[j + 1]) ++j;
var tmp = ary[i - 1]; ary[i - 1] = ary[j]; ary[j] = tmp;
var s = i;
var t = n - 1;
while (t - s > 0)
{
tmp = ary[t]; ary[t] = ary[s]; ary[s] = tmp;
++s; --t;
}
return true;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public (long, long) InvGCD(long a, long b)
{
a = (a + b) % b;
if (a == 0) return (b, 0);
var s = b;
var t = a;
var m0 = 0L;
var m1 = 1L;
while (t > 0)
{
var u = s / t; s -= t * u; m0 -= m1 * u;
var tmp = s; s = t; t = tmp;
tmp = m0; m0 = m1; m1 = tmp;
}
if (m0 < 0) m0 += b / s;
return (s, m0);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public (long, long) CRT(long[] r, long[] m)
{
var r0 = 0L;
var m0 = 1L;
for (var i = 0; i < m.Length; i++)
{
var m1 = m[i];
var r1 = (r[i] + m1) % m1;
if (m0 < m1)
{
r0 ^= r1; r1 ^= r0; r0 ^= r1;
m0 ^= m1; m1 ^= m0; m0 ^= m1;
}
if (m0 % m1 == 0)
{
if (r0 % m1 != r1) return (0, 0);
continue;
}
var gim = InvGCD(m0, m1);
var u1 = m1 / gim.Item1;
if ((r1 - r0) % gim.Item1 != 0) return (0, 0);
var x = (r1 - r0) / gim.Item1 % u1 * gim.Item2 % u1;
r0 += x * m0;
m0 *= u1;
if (r0 < 0) r0 += m0;
}
return (r0, m0);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static public long Phi(long n)
{
var ret = n;
for (var i = 2L; i * i <= n; ++i)
{
if (n % i == 0)
{
ret -= ret / i;
while (n % i == 0) n /= i;
}
}
if (n > 1) ret -= ret / n;
return ret;
}
[MethodImpl(MethodImplOptions.AggressiveInlining | MethodImplOptions.AggressiveOptimization)]
/// <summary>
/// a*i + b
/// i=0...n-1
/// 格子高さm
/// </summary>
static public long FloorSum(long n, long m, long a, long b)
{
unchecked
{
var ans = 0UL;
var un = (ulong)n;
var um = (ulong)m;
var ua = (ulong)a;
var ub = (ulong)b;
if (a < 0)
{
var am = a % m;
if (am < 0) am += m;
var a2 = (ulong)am;
ans -= un * (un - 1) / 2 * ((a2 + (ulong)(-a)) / um);
ua = a2;
}
if (b < 0)
{
var bm = b % m;
if (bm < 0) bm += m;
var b2 = (ulong)bm;
ans -= un * ((b2 + (ulong)(-b)) / um);
ub = b2;
}
while (true)
{
if (ua >= um)
{
ans += (un - 1) * un / 2 * (ua / um);
ua %= um;
}
if (ub >= um)
{
ans += un * (ub / um);
ub %= um;
}
var ymax = ua * un + ub;
if (ymax < um) return (long)ans;
un = ymax / um;
ub = ymax % um;
um ^= ua; ua ^= um; um ^= ua;
}
}
}
}
class LIB_FastIO
{
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public LIB_FastIO() { str = Console.OpenStandardInput(); }
readonly Stream str;
readonly byte[] buf = new byte[2048];
int len, ptr;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
byte read()
{
if (ptr >= len)
{
ptr = 0;
if ((len = str.Read(buf, 0, 2048)) <= 0)
{
return 0;
}
}
return buf[ptr++];
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
char Char()
{
byte b = 0;
do b = read();
while (b < 33 || 126 < b);
return (char)b;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
virtual public string Scan()
{
var sb = new StringBuilder();
for (var b = Char(); b >= 33 && b <= 126; b = (char)read())
sb.Append(b);
return sb.ToString();
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
virtual public long Long()
{
long ret = 0; byte b = 0; var ng = false;
do b = read();
while (b != '-' && (b < '0' || '9' < b));
if (b == '-') { ng = true; b = read(); }
for (; true; b = read())
{
if (b < '0' || '9' < b)
return ng ? -ret : ret;
else ret = (ret << 3) + (ret << 1) + b - '0';
}
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
virtual public double Double() { return double.Parse(Scan(), CultureInfo.InvariantCulture); }
}
class LIB_FastIODebug : LIB_FastIO
{
Queue<string> param = new Queue<string>();
[MethodImpl(MethodImplOptions.AggressiveInlining)]
string NextString() { if (param.Count == 0) foreach (var item in Console.ReadLine().Split(' ')) param.Enqueue(item); return param.Dequeue(); }
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public LIB_FastIODebug() { }
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public override string Scan() => NextString();
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public override long Long() => long.Parse(NextString());
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public override double Double() => double.Parse(NextString());
}
class LIB_LPSimplex
{
public class Subject
{
public LIB_Fraction[] ary;
public LIB_Fraction value;
public bool hasSlack;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public Subject(LIB_Fraction[] ary)
{
this.ary = ary;
value = 0;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
private void EqualTo(LIB_Fraction val)
{
// TODO public で呼ばれるとスラック変数が無いので基底決め打ちができない
value = val;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public void LessThan(LIB_Fraction val)
{
hasSlack = true;
EqualTo(val);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public void GreaterThan(LIB_Fraction val)
{
for (var i = 0; i < ary.Length; ++i) ary[i] *= -1;
LessThan(val * -1);
}
}
List<Subject> subs;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public LIB_LPSimplex()
{
subs = new List<Subject>();
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public Subject AddSubject(params LIB_Fraction[] ary)
{
var sub = new Subject(ary.ToArray());
subs.Add(sub);
return sub;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
private void NonNegative(params long[] idx)
{
// TODO 自由変数の取り扱い
if (idx.Length == 0) return;
var ary = Enumerable.Repeat(0, (int)idx.Max() + 1).Select(e => new LIB_Fraction(0, 1)).ToArray();
foreach (var item in idx)
{
ary[item] = 1;
AddSubject(ary).GreaterThan(0);
ary[item] = 0;
}
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public LIB_Fraction Minimize(params LIB_Fraction[] ary)
{
var maxlen = subs.Max(e => e.ary.Length);
var target = new LIB_Fraction[maxlen + subs.Count];
for (var i = 0; i < ary.Length; ++i)
{
target[i] = ary[i];
}
var targetConstant = new LIB_Fraction(0, 1);
var kitei = new long[subs.Count];
var subjects = new LIB_Fraction[subs.Count][];
var constants = new LIB_Fraction[subs.Count];
for (var i = 0; i < subs.Count; ++i)
{
var sub = subs[i];
kitei[i] = maxlen + i;
subjects[i] = new LIB_Fraction[maxlen + subs.Count];
constants[i] = sub.value;
for (var j = 0; j < sub.ary.Length; ++j)
{
subjects[i][j] = sub.ary[j];
}
subjects[i][kitei[i]] = 1;
}
while (true)
{
var minidx = -1;
var minvalue = new LIB_Fraction(0, 1);
for (var i = 0; i < target.Length; ++i)
{
if (minvalue > target[i])
{
minidx = i;
minvalue = target[i];
}
}
if (minidx == -1) break;
minvalue = new LIB_Fraction(1, 0);
var selectedSubIdx = -1;
for (var i = 0; i < subjects.Length; ++i)
{
var zouka = constants[i] / subjects[i][minidx];
if (minvalue > zouka)
{
minvalue = zouka;
selectedSubIdx = i;
}
}
kitei[selectedSubIdx] = minidx;
var selectedSub = subjects[selectedSubIdx];
var mul = 1 / selectedSub[minidx];
for (var i = 0; i < selectedSub.Length; ++i)
{
selectedSub[i] *= mul;
}
constants[selectedSubIdx] *= mul;
for (var i = 0; i < subjects.Length; ++i)
{
if (i == selectedSubIdx) continue;
mul = subjects[i][minidx];
for (var j = 0; j < subjects[i].Length; ++j)
{
subjects[i][j] -= selectedSub[j] * mul;
}
constants[i] -= constants[selectedSubIdx] * mul;
}
mul = target[minidx];
for (var i = 0; i < target.Length; ++i)
{
target[i] -= selectedSub[i] * mul;
}
targetConstant -= constants[selectedSubIdx] * mul;
}
return -1 * targetConstant;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public LIB_Fraction Maximize(params LIB_Fraction[] ary)
{
ary = ary.ToArray();
for (var i = 0; i < ary.Length; ++i) ary[i] *= -1;
return -1 * Minimize(ary);
}
}
}