結果
| 問題 |
No.112 ややこしい鶴亀算
|
| コンテスト | |
| ユーザー |
 tu-sa
|
| 提出日時 | 2018-06-29 05:58:02 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
CE
(最新)
AC
(最初)
|
| 実行時間 | - |
| コード長 | 6,312 bytes |
| コンパイル時間 | 1,273 ms |
| コンパイル使用メモリ | 124,280 KB |
| 最終ジャッジ日時 | 2024-11-14 20:31:02 |
| 合計ジャッジ時間 | 1,842 ms |
|
ジャッジサーバーID （参考情報） |
judge5 / judge2 |
（要ログイン）
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。（リジャッジ後のコンパイルエラーは公開されます）
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
コンパイルメッセージ
main.cpp: In function 'std::ostream& operator<<(std::ostream&, const iter_pair_t<T>&)':
main.cpp:82:33: error: 'ostream_iterator' was not declared in this scope
82 | std::copy(v.beg, v.end, ostream_iterator<typename decltype(v.beg)::reference>(out, " "));
| ^~~~~~~~~~~~~~~~
main.cpp:29:1: note: 'std::ostream_iterator' is defined in header '<iterator>'; did you forget to '#include <iterator>'?
28 | #include <array>
+++ |+#include <iterator>
29 | using namespace std;
main.cpp:82:85: error: expected '(' before '>' token
82 | std::copy(v.beg, v.end, ostream_iterator<typename decltype(v.beg)::reference>(out, " "));
| ^
| (
ソースコード
////////////////////////////////////////
/// tu3 pro-con template ///
////////////////////////////////////////
#include <cassert>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <iostream>
#include <iomanip>
#include <sstream>
#include <algorithm>
#include <numeric>
#include <functional>
#include <vector>
#include <queue>
#include <string>
#include <complex>
#include <stack>
#include <set>
#include <map>
#include <list>
#include <unordered_map>
#include <unordered_set>
#include <bitset>
#include <regex>
#include <type_traits>
#include <mutex>
#include <array>
using namespace std;
//// MACRO ////
#define countof(a) (sizeof(a)/sizeof(a[0]))
#define REP(i,n) for (int i = 0; i < (n); i++)
#define RREP(i,n) for (int i = (n)-1; i >= 0; i--)
#define FOR(i,s,n) for (int i = (s); i < (n); i++)
#define RFOR(i,s,n) for (int i = (n)-1; i >= (s); i--)
#define pos(c,i) c.being() + (i)
#define allof(c) c.begin(), c.end()
#define aallof(a) a, countof(a)
#define partof(c,i,n) c.begin() + (i), c.begin() + (i) + (n)
#define apartof(a,i,n) a + (i), a + (i) + (n)
typedef unsigned int uint;
typedef long long llong;
typedef unsigned long long ullong;
#define long long long
#define EPS 1e-9
#define INF (1L << 28)
#define LINF (1LL << 60)
#define PREDICATE(t,a,exp) [&](const t & a) -> bool { return exp; }
#define COMPARISON_T(t) bool(*)(const t &, const t &)
#define COMPARISON(t,a,b,exp) [&](const t & a, const t & b) -> bool { return exp; }
#define CONVERTER(TSrc,t,TDest,exp) [&](const TSrc &t)->TDest { return exp; }
inline int sign_of(double x) { return abs(x) < EPS ? 0 : x > 0 ? 1 : -1; }
inline bool inRange(int val, int min, int max) { return val >= min && val < max; }
inline bool inRange(double val, double min, double max) { return val - min > -EPS && val - max < EPS; }
inline bool inRange(int x, int y, int W, int H) { return x >= 0 && x < W && y >= 0 && y < H; } // W,H含まない
template <class T> struct vevector : public vector<vector<T>> { vevector(int n = 0, int m = 0, const T &initial = T()) : vector<vector<T>>(n, vector<T>(m, initial)) { } };
template <class T> struct vevevector : public vector<vevector<T>> { vevevector(int n = 0, int m = 0, int l = 0, const T &initial = T()) : vector<vevector<T>>(n, vevector<T>(m, l, initial)) { } };
template <class T> struct vevevevector : public vector<vevevector<T>> { vevevevector(int n = 0, int m = 0, int l = 0, int k = 0, const T &initial = T()) : vector<vevevector<T>>(n, vevevector<T>(m, l, k, initial)) { } };
//// i/o helper ////
namespace std {
template <class T1, class T2> inline istream & operator >> (istream & in, pair<T1, T2> &p) { in >> p.first >> p.second; return in; }
template <class T1, class T2> inline ostream & operator << (ostream &out, const pair<T1, T2> &p) { out << p.first << " " << p.second; return out; }
}
template <class T> T read() { T t; cin >> t; return t; }
template <class T> vector<T> read(int n) { vector<T> v; v.reserve(n); REP(i, n) { v.push_back(read<T>()); } return v; }
template <class T> vevector<T> read(int n, int m) { vevector<T> v; REP(i, n) v.push_back(read<T>(m)); return v; }
template <class T> vector<T> readjag() { return read<T>(read<int>()); }
template <class T> vevector<T> readjag(int n) { vevector<T> v; v.reserve(n); REP(i, n) v.push_back(readjag<T>()); return v; }
template <class T> struct iter_pair_t { T beg, end; };
template <class T> iter_pair_t<T> iter_pair(T beg, T end) { return iter_pair_t<T>{beg, end}; }
template <class T> ostream & operator << (ostream &out, const iter_pair_t<T> &v)
{
std::copy(v.beg, v.end, ostream_iterator<typename decltype(v.beg)::reference>(out, " "));
return out;
}
template <class T1> ostream & operator << (ostream &out, const vector<T1> &v) { return out << iter_pair(begin(v), end(v)); }
template <class T1> ostream & operator << (ostream &out, const set<T1> &v) { return out << iter_pair(begin(v), end(v)); }
template <class T1, class T2> ostream & operator << (ostream &out, const map<T1, T2> &v) { return out << iter_pair(begin(v), end(v)); }
struct _Reader { istream &cin; template <class T> _Reader operator ,(T &rhs) { cin >> rhs; return *this; } };
struct _Writer { ostream &cout; bool f{ false }; template <class T> _Writer operator ,(const T &rhs) { cout << (f ? " " : "") << rhs; f = true; return *this; } };
#define READ(t,...) t __VA_ARGS__; (_Reader{cin}), __VA_ARGS__
#define WRITE(...) (_Writer{cout}), __VA_ARGS__; cout << '\n'
#define DEBUG(...) (_Writer{cerr}), __VA_ARGS__; cerr << '\n'
void solve();
int main()
{
cin.tie(0);
ios_base::sync_with_stdio(false);
cout << setprecision(std::numeric_limits<double>::max_digits10);
solve();
return 0;
}
enum struct GaussJordanResult
{
One,
Inconsistent,
Ambiguous,
};
GaussJordanResult gaussJordan(vevector<double> &m)
{
int R = m.size();
int C = m[0].size() - 1;
int pi = 0;
REP(pj, C)
{
// ピボット選択
FOR(i, pi + 1, R)
{
if (abs(m[i][pj]) > abs(m[pi][pj]))
{
m[pi].swap(m[i]);
}
}
if (m[pi][pj])
{
{
double r = m[pi][pj];
REP(j, C + 1) { m[pi][j] /= r; }
}
REP(i, R)
{
if (i == pi) { continue; }
double r = m[i][pj];
REP(j, C + 1) { m[i][j] = m[i][j] - r * m[pi][j]; }
}
pi++;
}
}
{
// 0 = !0 のような式が残ってると解なし。
FOR(i, pi, R)
{
if (abs(m[i][C]) > EPS)
{
return GaussJordanResult::Inconsistent;
}
}
return pi == C ? GaussJordanResult::One : GaussJordanResult::Ambiguous;
}
}
template<class T>
vector<T> gaussJordanAns(const vevector<T> &m)
{
// 行列のランクが列数より低い → 解が一意に定まらない
// 1 d 0 | a
// 0 0 1 | b のようなときは、x[0] = a, x[1] = 0, x[2] = b が解の1つ
// 0 0 0 | 0
int C = m[0].size() - 1;
vector<T> ans(C);
int i = 0;
REP(j, C)
{
while (abs(m[i][j]) < EPS) { j++; }
ans[j] = m[i][C];
i++;
}
return ans;
}
////////////////////
/// template end ///
////////////////////
void solve()
{
READ(int, N);
auto A = read<int>(N);
vevector<double> m(N, N + 1, 1);
REP(i, N)
{
m[i][i] = 0;
m[i][N] = A[i];
}
auto r = gaussJordan(m);
if (r == GaussJordanResult::One)
{
map<int, int> ma;
for(auto d : gaussJordanAns(m))
{
ma[int(d + EPS)]++;
}
WRITE(ma[2], ma[4]);
}
else
{
throw r;
}
}