結果

問題 No.498 ワープクリスタル (給料日編)
ユーザー algon_320
提出日時 2017-03-25 02:34:00
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 8 ms / 2,000 ms
コード長 3,565 bytes
コンパイル時間 1,799 ms
コンパイル使用メモリ 171,920 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-06 03:26:14
合計ジャッジ時間 2,503 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define int long long
using ll=long long;
using vi=vector<int>;
using vl=vector<long long>;
using pii=pair<int,int>;
using pll=pair<long long,long long>;
#define ITR(i,c) for(auto i=begin(c);i!=end(c);++i)
#define FORE(x,c) for(auto &x:c)
#define FOR(i,a,n) for(int i=a,i##_len=(int)(n);i<i##_len;++i)
#define REP(i,n) FOR(i,0,n)
#define RREP(i,n) for(int i=(int)(n);i>=0;--i)
#define ALL(c) begin(c),end(c)
#define RALL(c) rbegin(c),rend(c)   // c++14
#define SZ(c) ((int)c.size())
#define EXIST(c,x) (c.find(x)!=end(c))
#define OUTOFRANGE(y,x,h,w) (y<0||x<0||y>=h||x>=w)
#define dump(...)
const int DX[9]={0,1,0,-1,1,1,-1,-1,0},DY[9]={-1,0,1,0,-1,1,1,-1,0};
#define INF (1001001001)
#define INFLL (1001001001001001001ll)
template<class T> ostream& operator << (ostream &os,const vector<T> &v) {
    ITR(i,v) os << *i << (i==end(v)-1 ? "" : "\n"); return os; }
template<class T> istream& operator >> (istream &is,vector<T> &v) {
    ITR(i,v) is >> * i; return is; }
template<class T> istream& operator >> (istream &is, pair<T,T> &p) {
        is >> p.first >> p.second; return is; }
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return 1;}return 0;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return 1;}return 0;}
//------------------------------------------------------------------------------
struct before_main_function {
    before_main_function() {
        #ifdef int
            #undef INF
            #define INF INFLL
            #define stoi stoll
        #endif
        cin.tie(0);ios::sync_with_stdio(false);
        cout<<setprecision(15)<<fixed;
    }
} before_main_function;
//------------------------------------------------------------------------------

const int MOD=1e9+7;
const int MAX_N=10000;
ll fact[MAX_N+1];
ll fact_inv[MAX_N+1];
ll mod_pow(ll a,ll b) {
    if(a==0) return 1;
    a%=MOD;
    ll product=1;
    while(b>0) {
        if(b&1) {
            product*=a;
            product%=MOD;
        }
        a*=a;
        a%=MOD;
        b>>=1;
    }
    return product%MOD;
}
void init_fact() {
    fact[0]=1, fact_inv[0]=1;
    FOR(i,1,MAX_N) {
        fact[i]=(fact[i-1]*i)%MOD;
        fact_inv[i]=mod_pow(fact[i],MOD-2);
    }
}
signed main() {
    init_fact();
    int gx,gy,K;
    cin>>gx>>gy>>K;
    vi x(K),y(K),N(5,0);
    REP(i,K) {
        cin>>x[i]>>y[i]>>N[i];
    }
    int ans=0;
    REP(i,N[0]+1) {
        REP(j,N[1]+1) {
            REP(k,N[2]+1) {
                REP(l,N[3]+1) {
                    REP(m,N[4]+1) {
                        int X=0,Y=0;
                        if(N[0]!=0) X+=x[0]*i, Y+=y[0]*i;
                        if(N[1]!=0) X+=x[1]*j, Y+=y[1]*j;
                        if(N[2]!=0) X+=x[2]*k, Y+=y[2]*k;
                        if(N[3]!=0) X+=x[3]*l, Y+=y[3]*l;
                        if(N[4]!=0) X+=x[4]*m, Y+=y[4]*m;
                        if(X==gx && Y==gy) {
                            int c=1;
                            c*=fact[i+j+k+l+m];
                            c*=fact_inv[i];
                            c%=MOD;
                            c*=fact_inv[j];
                            c%=MOD;
                            c*=fact_inv[k];
                            c%=MOD;
                            c*=fact_inv[l];
                            c%=MOD;
                            c*=fact_inv[m];
                            c%=MOD;
                            ans+=c;
                        }
                    }
                }
            }
        }
    }
    cout<<ans%MOD<<endl;
    return 0;
}
0