結果

問題 No.122 傾向と対策：門松列（その３）
コンテスト
ユーザー chocorusk
提出日時 2019-12-29 12:10:16
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 138 ms / 5,000 ms
コード長 2,505 bytes
コンパイル時間 1,113 ms
コンパイル使用メモリ 122,712 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2025-01-01 17:51:43
合計ジャッジ時間 2,029 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge7
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 8
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <time.h>
#include <stack>
#include <array>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
const ll MOD=1e9+7;
ll solve(vector<int> l, vector<int> r){
  int m=l.size();
  vector<int> v(2*m);
  for(int i=0; i<m; i++){
    if(l[i]>=r[i]) return 0;
    v[2*i]=l[i], v[2*i+1]=r[i];
  }
  sort(v.begin(), v.end());
  v.erase(unique(v.begin(), v.end()), v.end());
  int n=v.size();
  for(int i=0; i<m; i++){
    l[i]=lower_bound(v.begin(), v.end(), l[i])-v.begin();
    r[i]=lower_bound(v.begin(), v.end(), r[i])-v.begin();
  }
  int a[5];
  ll ret=0;
  auto dfs=[&](auto dfs, int k)->void{
    if(k==m){
      int cnt[10]={};
      for(int i=0; i<m; i++){
        cnt[a[i]]++;
      }
      ll p=1;
      for(int i=0; i<n-1; i++){
        for(int j=0; j<cnt[i]; j++){
          p*=(v[i+1]-v[i]-j);
          p%=MOD;
        }
      }
      ret+=p;
      if(ret>=MOD) ret-=MOD;
      return;
    }
    for(int i=l[k]; i<r[k]; i++){
      a[k]=i;
      dfs(dfs, k+1);
    }
  };
  dfs(dfs, 0);
  return ret;
}
int main()
{
	int a[7], b[7];
  for(int i=0; i<7; i++) cin>>a[i]>>b[i];
  ll ans=0;
  for(int i=2; i<=20000; i++){
    vector<int> l(3), r(3);
    for(int j=0; j<3; j++){
      l[j]=max(i, a[2*j+1]);
      r[j]=b[2*j+1]+1;
    }
    ll s=solve(l, r);
    for(int j=0; j<3; j++){
      l[j]=max(i+1, a[2*j+1]);
      r[j]=b[2*j+1]+1;
    }
    s=(s-solve(l, r)+MOD)%MOD;
    if(s==0) continue;
    l.resize(4); r.resize(4);
    for(int j=0; j<4; j++){
      l[j]=a[2*j];
      r[j]=min(i, b[2*j]+1);
    }
    (s*=solve(l, r))%=MOD;
    (ans+=s)%=MOD;
  }
  for(int i=2; i<=20000; i++){
    vector<int> l(4), r(4);
    for(int j=0; j<4; j++){
      l[j]=max(i, a[2*j]);
      r[j]=b[2*j]+1;
    }
    ll s=solve(l, r);
    for(int j=0; j<4; j++){
      l[j]=max(i+1, a[2*j]);
      r[j]=b[2*j]+1;
    }
    s=(s-solve(l, r)+MOD)%MOD;
    if(s==0) continue;
    l.resize(3); r.resize(3);
    for(int j=0; j<3; j++){
      l[j]=a[2*j+1];
      r[j]=min(i, b[2*j+1]+1);
    }
    (s*=solve(l, r))%=MOD;
    (ans+=s)%=MOD;
  }
  cout<<ans<<endl;
	return 0;
}
0