結果

問題 No.492 IOI数列
ユーザー torus711torus711
提出日時 2017-03-10 23:27:41
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 1,000 ms
コード長 5,534 bytes
コンパイル時間 1,296 ms
コンパイル使用メモリ 109,900 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-09-06 14:21:05
合計ジャッジ時間 2,147 ms
ジャッジサーバーID
(参考情報)
judge14 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 2 ms
4,384 KB
testcase_04 AC 2 ms
4,384 KB
testcase_05 AC 1 ms
4,384 KB
testcase_06 AC 2 ms
4,380 KB
testcase_07 AC 1 ms
4,380 KB
testcase_08 AC 1 ms
4,380 KB
testcase_09 AC 1 ms
4,384 KB
testcase_10 AC 2 ms
4,384 KB
testcase_11 AC 1 ms
4,380 KB
testcase_12 AC 1 ms
4,380 KB
testcase_13 AC 1 ms
4,380 KB
testcase_14 AC 1 ms
4,384 KB
testcase_15 AC 2 ms
4,380 KB
testcase_16 AC 1 ms
4,384 KB
testcase_17 AC 1 ms
4,384 KB
testcase_18 AC 2 ms
4,380 KB
testcase_19 AC 1 ms
4,380 KB
testcase_20 AC 1 ms
4,384 KB
testcase_21 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <iomanip>
#include <sstream>
#include <vector>
#include <string>
#include <set>
#include <unordered_set>
#include <map>
#include <unordered_map>
#include <stack>
#include <queue>
#include <deque>
#include <algorithm>
#include <functional>
#include <iterator>
#include <limits>
#include <numeric>
#include <utility>
#include <type_traits>
#include <cmath>
#include <cassert>
#include <cstdio>

using namespace std;
using namespace placeholders;

using LL = long long;
using ULL = unsigned long long;
using VI = vector< int >;
using VVI = vector< vector< int > >;
using VS = vector< string >;
using ISS = istringstream;
using OSS = ostringstream;
using PII = pair< int, int >;
using VPII = vector< pair< int, int > >;
template < typename T = int > using VT = vector< T >;
template < typename T = int > using VVT = vector< vector< T > >;
template < typename T = int > using LIM = numeric_limits< T >;

template < typename T > inline istream& operator>>( istream &s, vector< T > &v ){ for ( T &t : v ) { s >> t; } return s; }
template < typename T > inline ostream& operator<<( ostream &s, const vector< T > &v ){ for ( int i = 0; i < int( v.size() ); ++i ){ s << ( " " + !i ) << v[i]; } return s; }

template < typename T > struct getv_fmt;
template <> struct getv_fmt<       int >{ static constexpr const char *fmt = "%d"; };
template <> struct getv_fmt< long long >{ static constexpr const char *fmt = "%lld"; };
template < typename T > void getv( std::vector< T > &v ){ for_each( begin( v ), end( v ), []( T &a ){ scanf( getv_fmt< T >::fmt, &a ); } ); };

template < typename T > inline T fromString( const string &s ) { T res; istringstream iss( s ); iss >> res; return res; }
template < typename T > inline string toString( const T &a ) { ostringstream oss; oss << a; return oss.str(); }

#define NUMBERED( name, number ) NUMBERED2( name, number )
#define NUMBERED2( name, number ) name ## _ ## number
#define REP1( n ) REP2( NUMBERED( REP_COUNTER, __LINE__ ), n )
#define REP2( i, n ) REP3( i, 0, n )
#define REP3( i, m, n ) for ( int i = ( int )( m ); i < ( int )( n ); ++i )
#define GET_REP( a, b, c, F, ... ) F
#define REP( ... ) GET_REP( __VA_ARGS__, REP3, REP2, REP1 )( __VA_ARGS__ )
#define FOR( e, c ) for ( auto &&e : c )
#define ALL( c ) begin( c ), end( c )
#define AALL( a ) ( remove_all_extents< decltype( a ) >::type * )a, ( remove_all_extents< decltype( a ) >::type * )a + sizeof( a ) / sizeof( remove_all_extents< decltype( a ) >::type )
#define DRANGE( c, p ) begin( c ), begin( c ) + ( p ), end( c )

#define SZ( v ) ( (int)( v ).size() )
#define EXIST( c, e ) ( ( c ).find( e ) != ( c ).end() )

template < typename T > inline bool chmin( T &a, const T &b ){ if ( b < a ) { a = b; return true; } return false; }
template < typename T > inline bool chmax( T &a, const T &b ){ if ( a < b ) { a = b; return true; } return false; }

#define PB push_back
#define EM emplace
#define EB emplace_back
#define BI back_inserter

#define MP make_pair
#define fst first
#define snd second

#define DUMP( x ) cerr << #x << " = " << ( x ) << endl

template < typename Number > using Vector = vector< Number >;
template < typename Number > using Matrix = vector< vector< Number > >;

// 単位行列
// O( N )
template < typename Number >
Matrix< Number > identify_matrix( const int N )
{
	Matrix< Number > res( N, Vector< Number >( N ) );
	for ( int i = 0; i < N; ++i )
	{
		res[i][i] = 1;
	}
	return res;
}

// 行列×ベクトル
template < typename Number >
Vector< Number > operator*( const Matrix< Number > &m, const Vector< Number > &v )
{
	if ( m[0].size() != v.size() )
	{
		throw "Matrix * Vector : size error";
	}

	Vector< Number > res( m.size() );
	for ( int i = 0; i < static_cast<int>( m.size() ); ++i )
	{
		for ( int k = 0; k < static_cast<int>( v.size() ); ++k )
		{
			res[i] += m[i][k] * v[k];
		}
	}
	return res;
}

// 行列同士の積
// O( N^3 )
template < typename Number >
Matrix< Number > matrix_multiply( const Matrix< Number > &a, const Matrix< Number > &b, unsigned int mod = 0 )
{
	assert( a[0].size() == b.size() );

	Matrix< Number > res( a.size(), Vector< Number >( b[0].size() ) );
	for ( int i = 0; i < static_cast<int>( a.size() ); ++i )
	{
		for ( int j = 0; j < static_cast<int>( b[0].size() ); ++j )
		{
			for ( int k = 0; k < static_cast<int>( b.size() ); ++k )
			{
				( res[i][j] += a[i][k] * b[k][j] ) %= mod;
			}
		}
	}
	return res;
}
template < typename Number >
inline const Matrix< Number > operator*( const Matrix< Number > &a, const Matrix< Number > &b )
{
	return matrix_multiply( a, b );
}

// 行列累乗
// OO( N^3 log n )
template < typename Number >
Matrix< Number > power( Matrix< Number > a, long long n, long long mod = numeric_limits< long long >::max() )
{
	Matrix< Number > res = identify_matrix< Number >( a.size() );

	while ( n )
	{
		if ( n & 1 )
		{
			res = matrix_multiply( res, a, mod 	);
		}
		a = matrix_multiply( a, a, mod );
		n >>= 1;
	}

	return res;
}

template < typename Number >
Number solve( const Matrix< Number > &A, const LL N, LL MOD )
{
	const Matrix< Number > AA = power( A, N, MOD );
	return ( AA * VT< LL >( { 1, 1 } ) )[1] % MOD;
}

int main()
{
	cin.tie( 0 );
	ios::sync_with_stdio( false );
	cout << setprecision( 12 ) << fixed;

	LL N;
	cin >> N;

	const Matrix< LL > A( {
			{ 1LL, 0LL },
			{ 1LL, 100LL } } );
	
	cout << solve( A, N - 1, 1000000007 ) << endl;

	N %= 11;
	if ( N == 0 )
	{
		cout << 0 << endl;
		return 0;
	}

	REP( i, N - 1 )
	{
		cout << 10;
	}
	cout << 1 << endl;

	return 0;
}
0