結果
| 問題 |
No.895 MESE
|
| コンテスト | |
| ユーザー |
mban
|
| 提出日時 | 2019-09-27 22:48:30 |
| 言語 | C#(csc) (csc 3.9.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 6,577 bytes |
| コンパイル時間 | 1,188 ms |
| コンパイル使用メモリ | 112,908 KB |
| 実行使用メモリ | 29,568 KB |
| 最終ジャッジ日時 | 2024-09-25 00:14:05 |
| 合計ジャッジ時間 | 12,281 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 WA * 1 |
| other | AC * 2 WA * 24 |
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc) Copyright (C) Microsoft Corporation. All rights reserved.
ソースコード
using System;
using System.Collections.Generic;
using System.Numerics;
public class Program
{
private int A, B, C;
private ModInt[] Fact;
public void Solve()
{
var sc = new Scanner();
A = sc.NextInt();
B = sc.NextInt();
C = sc.NextInt();
int sum = A + B + C;
// aの上bit > b > c
Fact = new ModInt[300001];
Fact[0] = 1;
for (int i = 1; i <= 300000; i++)
{
Fact[i] = Fact[i - 1] * i;
}
ModInt ans = 0;
// Zの最上位bit
for (int i = C - 1; i < sum - 2; i++)
{
// Zの最上位bitの和
ModInt mostZ = ModInt.Pow(2 ,i);
ModInt mostZSum = mostZ * Com(i, C - 1);
// Zの最上位以外の和
ModInt t;
if (C == 1)
t = 0;
else
{
t = Com(i - 1, C - 2) * ModInt.Pow(2, i - 1);
}
ModInt ZSum = mostZSum + t;
// abのパターン数
ModInt patXY = Com(sum - C - 1, A - 1) - Com(i,B);
ans += ZSum * patXY;
//Console.WriteLine($"{i} {ZSum} {patXY}");
}
Console.WriteLine(ans);
}
private ModInt Com(int n, int m)
{
if (n < m) return 0;
return Fact[n] * ModInt.Inverse(Fact[n - m] * Fact[m]);
}
public static void Main(string[] args)
{
new Program().Solve();
}
}
#region ModInt
/// <summary>
/// [0,<see cref="Mod"/>) までの値を取るような数
/// </summary>
public struct ModInt
{
/// <summary>
/// 剰余を取る値.
/// </summary>
public const long Mod = (int) 1e9 + 7;
/// <summary>
/// 実際の数値.
/// </summary>
public long num;
/// <summary>
/// 値が <paramref name="n"/> であるようなインスタンスを構築します.
/// </summary>
/// <param name="n">インスタンスが持つ値</param>
/// <remarks>パフォーマンスの問題上,コンストラクタ内では剰余を取りません.そのため,<paramref name="n"/> ∈ [0,<see cref="Mod"/>) を満たすような <paramref name="n"/> を渡してください.このコンストラクタは O(1) で実行されます.</remarks>
public ModInt(long n)
{
num = n;
}
/// <summary>
/// このインスタンスの数値を文字列に変換します.
/// </summary>
/// <returns>[0,<see cref="Mod"/>) の範囲内の整数を 10 進表記したもの.</returns>
public override string ToString()
{
return num.ToString();
}
public static ModInt operator +(ModInt l, ModInt r)
{
l.num += r.num;
if (l.num >= Mod) l.num -= Mod;
return l;
}
public static ModInt operator -(ModInt l, ModInt r)
{
l.num -= r.num;
if (l.num < 0) l.num += Mod;
return l;
}
public static ModInt operator *(ModInt l, ModInt r)
{
return new ModInt(l.num * r.num % Mod);
}
public static implicit operator ModInt(long n)
{
n %= Mod;
if (n < 0) n += Mod;
return new ModInt(n);
}
/// <summary>
/// 与えられた 2 つの数値からべき剰余を計算します.
/// </summary>
/// <param name="v">べき乗の底</param>
/// <param name="k">べき指数</param>
/// <returns>繰り返し二乗法により O(N log N) で実行されます.</returns>
public static ModInt Pow(ModInt v, long k)
{
return Pow(v.num, k);
}
/// <summary>
/// 与えられた 2 つの数値からべき剰余を計算します.
/// </summary>
/// <param name="v">べき乗の底</param>
/// <param name="k">べき指数</param>
/// <returns>繰り返し二乗法により O(N log N) で実行されます.</returns>
public static ModInt Pow(long v, long k)
{
long ret = 1;
for (k %= Mod - 1; k > 0; k >>= 1, v = v * v % Mod)
if ((k & 1) == 1)
ret = ret * v % Mod;
return new ModInt(ret);
}
/// <summary>
/// 与えられた数の逆元を計算します.
/// </summary>
/// <param name="v">逆元を取る対象となる数</param>
/// <returns>逆元となるような値</returns>
/// <remarks>法が素数であることを仮定して,フェルマーの小定理に従って逆元を O(log N) で計算します.</remarks>
public static ModInt Inverse(ModInt v)
{
return Pow(v, Mod - 2);
}
}
#endregion
class Scanner
{
public Scanner()
{
_pos = 0;
_line = new string[0];
}
const char Separator = ' ';
private int _pos;
private string[] _line;
#region スペース区切りで取得
public string Next()
{
if (_pos >= _line.Length)
{
_line = Console.ReadLine().Split(Separator);
_pos = 0;
}
return _line[_pos++];
}
public int NextInt()
{
return int.Parse(Next());
}
public long NextLong()
{
return long.Parse(Next());
}
public double NextDouble()
{
return double.Parse(Next());
}
#endregion
#region 型変換
private int[] ToIntArray(string[] array)
{
var result = new int[array.Length];
for (int i = 0; i < array.Length; i++)
{
result[i] = int.Parse(array[i]);
}
return result;
}
private long[] ToLongArray(string[] array)
{
var result = new long[array.Length];
for (int i = 0; i < array.Length; i++)
{
result[i] = long.Parse(array[i]);
}
return result;
}
private double[] ToDoubleArray(string[] array)
{
var result = new double[array.Length];
for (int i = 0; i < array.Length; i++)
{
result[i] = double.Parse(array[i]);
}
return result;
}
#endregion
#region 配列取得
public string[] Array()
{
if (_pos >= _line.Length)
_line = Console.ReadLine().Split(Separator);
_pos = _line.Length;
return _line;
}
public int[] IntArray()
{
return ToIntArray(Array());
}
public long[] LongArray()
{
return ToLongArray(Array());
}
public double[] DoubleArray()
{
return ToDoubleArray(Array());
}
#endregion
}
mban