結果

問題 No.2272 多項式乗算 mod 258280327
ユーザー 👑 p-adic
提出日時 2023-04-09 16:16:40
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 2,197 bytes
コンパイル時間 663 ms
コンパイル使用メモリ 69,360 KB
最終ジャッジ日時 2025-02-12 04:15:10
ジャッジサーバーID
(参考情報)
judge2 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 27 TLE * 6
権限があれば一括ダウンロードができます

ソースコード

diff #

// 誤解法(O(NM)愚直解)チェック
#include <iostream>
#include <stdio.h>
#include <stdint.h>
#include <cassert>
using namespace std;

using ll = long long;

#define MAIN main
#define TYPE_OF( VAR ) remove_const<remove_reference<decltype( VAR )>::type >::type
#define UNTIE ios_base::sync_with_stdio( false ); cin.tie( nullptr )
#define CEXPR( LL , BOUND , VALUE ) constexpr 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 FOREQ( VAR , INITIAL , FINAL ) for( TYPE_OF( FINAL ) VAR = INITIAL ; VAR <= FINAL ; VAR ++ )
#define QUIT return 0
#define COUT( ANSWER ) cout << ( ANSWER ) << "\n";
#define RETURN( ANSWER ) COUT( ANSWER ); QUIT

int MAIN()
{
  UNTIE;
  CEXPR( int , power3_max , 531441 ); // 400000を超える最小の3羃3^12
  CEXPR( int , bound_deg , 200000 );
  CEXPR( ll , bound_coef , 1000000000000000000 );
  CEXPR( ll , P , 258280327 );
  CIN_ASSERT( deg_F , 0 , bound_deg );
  ll F[power3_max];
  FOREQ( i , 0 , deg_F ){
    CIN_ASSERT( F_i , 0 , bound_coef );
    F[i] = move( F_i );
  }
  CIN_ASSERT( deg_G , 0 , bound_deg );
  ll G[power3_max];
  FOREQ( i , 0 , deg_G ){
    CIN_ASSERT( G_i , 0 , bound_coef );
    G[i] = move( G_i );
  }
  assert( deg_F > 0 ? F[deg_F] > 0 : true );
  assert( deg_G > 0 ? G[deg_G] > 0 : true );
  if( deg_F == 0 ){
    if( F[0] == 0 ){
      cout << 0 << "\n";
      cout << 0 << "\n";
      QUIT;
    }
  }
  if( deg_G == 0 ){
    if( G[0] == 0 ){
      cout << 0 << "\n";
      cout << 0 << "\n";
      QUIT;
    }
  }
  FOREQ( i , 0 , deg_F ){
    F[i] %= P;
  }
  FOREQ( i , 0 , deg_G ){
    G[i] %= P;
  }
  int deg_FG = deg_F + deg_G;
  cout << deg_FG << "\n";
  FOR( d , 0 , deg_FG ){
    ll answer = 0;
    int i_min = max( 0 , d - deg_G );
    int i_max = min( deg_F , d );
    FOREQ( i , i_min , i_max ){
      ( answer += F[i] * G[d - i] ) %= P;
    }
    cout << answer << " ";
  }
  RETURN( ( F[deg_F] * G[deg_G] ) % P);
}
0