結果

問題 No.329 全射
ユーザー te-sh
提出日時 2017-07-04 14:36:47
言語 D
(dmd 2.109.1)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 3,053 bytes
コンパイル時間 554 ms
コンパイル使用メモリ 134,404 KB
最終ジャッジ日時 2024-11-14 20:05:47
合計ジャッジ時間 1,801 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
/home/linuxbrew/.linuxbrew/opt/dmd/include/dlang/dmd/std/format/internal/write.d(143): Error: cannot implicitly convert expression `obj` of type `const(FactorRing!(1000000007, true))` to `int`
/home/linuxbrew/.linuxbrew/opt/dmd/include/dlang/dmd/std/format/write.d(1239): Error: template instance `std.format.internal.write.formatValueImpl!(LockingTextWriter, FactorRing!(1000000007, true), char)` error instantiating
/home/linuxbrew/.linuxbrew/opt/dmd/include/dlang/dmd/std/format/write.d(632):        instantiated from here: `formatValue!(LockingTextWriter, FactorRing!(1000000007, true), char)`
/home/linuxbrew/.linuxbrew/opt/dmd/include/dlang/dmd/std/stdio.d(1759):        instantiated from here: `formattedWrite!(LockingTextWriter, char, FactorRing!(1000000007, true))`
/home/linuxbrew/.linuxbrew/opt/dmd/include/dlang/dmd/std/stdio.d(4277):        instantiated from here: `write!(FactorRing!(1000000007, true), char)`
Main.d(30):        instantiated from here: `writeln!(FactorRing!(1000000007, true))`

ソースコード

diff #
プレゼンテーションモードにする

import std.algorithm, std.conv, std.range, std.stdio, std.string;
import std.container; // SList, DList, BinaryHeap
const p = 10 ^^ 9 + 7;
alias Graph!(int, size_t) graph;
alias FactorRing!(p, true) mint;
void main()
{
auto rd1 = readln.split.to!(size_t[]), n = rd1[0], m = rd1[1];
auto w = readln.split.to!(int[]);
auto g = new int[][](n, n);
foreach (_; 0..m) {
auto rd2 = readln.split.to!(size_t[]), i = rd2[0]-1, j = rd2[1]-1;
g[i][j] = w[i];
}
auto d = graph.floydWarshal(g);
auto ft = factTable!mint(1000);
auto st = starling2!mint(1000);
auto r = mint(0);
foreach (i; 0..n)
foreach (j; 0..n)
if (i == j || d[i][j] < graph.inf && d[i][j] >= w[j])
r += ft[w[j]] * st[w[i]][w[j]];
writeln(r);
}
pure T[] factTable(T)(size_t n)
{
auto t = new T[](n + 1);
t[0] = T(1);
foreach (i; 1..n+1)
t[i] = t[i-1] * i;
return t;
}
pure T[][] starling2(T)(size_t n)
{
auto t = new T[][](n + 1);
t[0] = new T[](1);
t[0][0] = 1;
foreach (i; 1..n+1) {
t[i] = new T[](i + 1);
t[i][0] = 0;
t[i][$-1] = 1;
foreach (j; 1..i)
t[i][j] = t[i - 1][j - 1] + j.to!T * t[i - 1][j];
}
return t;
}
template Graph(Wt, Node, Wt _inf = 10 ^^ 9, Node _sent = Node.max)
{
import std.algorithm, std.array, std.conv;
const inf = _inf, sent = _sent;
Wt[][] floydWarshal(Wt[][] g)
{
Wt[][] dist;
Node[][] inter;
floydWarshal(g, dist, inter);
return dist;
}
void floydWarshal(Wt[][] g, out Wt[][] dist, out Node[][] inter)
{
auto n = g.length;
dist = g.map!(i => i.dup).array;
inter = new Node[][](n, n);
foreach (i; 0..n) inter[i][] = sent;
foreach (k; 0..n)
foreach (i; 0..n)
foreach (j; 0..n)
if (dist[i][k] && dist[k][j] && (!dist[i][j] || dist[i][j] < min(dist[i][k], dist[k][j]))) {
dist[i][j] = min(dist[i][k], dist[k][j]);
inter[i][j] = k.to!Node;
}
}
}
struct FactorRing(int m, bool pos = false)
{
long v;
@property int toInt() { return v.to!int; }
alias toInt this;
this(T)(T _v) { v = mod(_v); }
ref FactorRing!(m, pos) opAssign(int _v)
{
v = mod(_v);
return this;
}
pure auto mod(long _v) const
{
static if (pos) return _v % m;
else return (_v % m + m) % m;
}
pure auto opBinary(string op: "+")(long rhs) const { return FactorRing!(m, pos)(v + rhs); }
pure auto opBinary(string op: "-")(long rhs) const { return FactorRing!(m, pos)(v - rhs); }
pure auto opBinary(string op: "*")(long rhs) const { return FactorRing!(m, pos)(v * rhs); }
pure auto opBinary(string op)(FactorRing!(m, pos) rhs) const
if (op == "+" || op == "-" || op == "*") { return opBinary!op(rhs.v); }
auto opOpAssign(string op: "+")(long rhs) { v = mod(v + rhs); }
auto opOpAssign(string op: "-")(long rhs) { v = mod(v - rhs); }
auto opOpAssign(string op: "*")(long rhs) { v = mod(v * rhs); }
auto opOpAssign(string op)(FactorRing!(m, pos) rhs)
if (op == "+" || op == "-" || op == "*") { return opOpAssign!op(rhs.v); }
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0