#pragma GCC optimize ( "O3" ) #pragma GCC optimize( "unroll-loops" ) #pragma GCC target ( "sse4.2,fma,avx2,popcnt,lzcnt,bmi2" ) #include using namespace std; using ll = long long; #define TYPE_OF( VAR ) decay_t #define UNTIE ios_base::sync_with_stdio( false ); cin.tie( nullptr ) #define CEXPR( LL , BOUND , VALUE ) constexpr const LL BOUND = VALUE #define CIN( LL , A ) LL A; cin >> A #define ASSERT( A , MIN , MAX ) assert( ( MIN ) <= A && A <= ( MAX ) ) #define CIN_ASSERT( A , MIN , MAX ) CIN( TYPE_OF( MAX ) , A ); ASSERT( A , MIN , MAX ) #define FOR( VAR , INITIAL , FINAL_PLUS_ONE ) for( TYPE_OF( FINAL_PLUS_ONE ) VAR = INITIAL ; VAR < FINAL_PLUS_ONE ; VAR ++ ) #define QUIT return 0 #define COUT( ANSWER ) cout << ( ANSWER ) << "\n" #define RETURN( ANSWER ) COUT( ANSWER ); QUIT int main() { UNTIE; CEXPR( int , bound_N , 20 ); CIN_ASSERT( N , 1 , bound_N ); CEXPR( ll , bound_B , 1000000000 ); CIN_ASSERT( B , 1 , bound_B ); CEXPR( ll , bound_Aij , 1000000000000000000 ); ll A[bound_N][bound_N]; FOR( i , 0 , N ){ ll ( &Ai )[bound_N] = A[i]; FOR( j , 0 , N ){ CIN_ASSERT( Aij , 0 , bound_Aij ); Ai[j] = Aij %= B; } } CEXPR( uint , bound_power_N , 1 << bound_N ); uint power_N = 1 << N; ll det[2][bound_power_N] = {}; ll ( &e )[bound_power_N] = det[0]; ll ( &o )[bound_power_N] = det[1]; e[0] = 1; o[0] = 0; FOR( d , 1 , power_N ){ ll& ed = e[d] = 0; ll& od = o[d] = 0; int n = 0; uint d_copy = d; while( d_copy != 0 ){ d_copy ^= d_copy & -d_copy; n++; } d_copy = d; ll ( &Ai )[bound_N] = A[N - n]; bool even = true; FOR( j , 0 , N ){ uint digit = 1 << j; if( ( d & digit ) != 0 ){ uint d_sub = d ^ digit; ll& Aij = Ai[j]; if( even ){ ( ed += Aij * e[d_sub] ) %= B; ( od += Aij * o[d_sub] ) %= B; } else { ( ed += Aij * o[d_sub] ) %= B; ( od += Aij * e[d_sub] ) %= B; } even = ! even; } } } RETURN( o[power_N - 1] ); }