結果

問題 No.2017 Mod7 Parade
ユーザー HIcoder
提出日時 2023-07-16 19:12:39
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 330 ms / 2,000 ms
コード長 1,609 bytes
コンパイル時間 1,390 ms
コンパイル使用メモリ 111,680 KB
最終ジャッジ日時 2025-02-15 15:17:20
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<set> 
#include<algorithm>
#include<vector>
#include<string>
#include<set>
#include<map>
#include<numeric>
#include<queue>
#include<cmath>
#include <unordered_set>
using namespace std;
typedef long long ll;
const ll INF=1LL<<60;
typedef pair<ll,ll> P;
typedef pair<int,P> PP;
const ll MOD=1e9+7;

ll mod_pow(ll x,ll y,ll mod){
    ll res=1;
    x%=mod;
    while(y>0){
        if(y&1){
            res*=x;
            res%=mod;
        }
        x*=x;
        x%=mod;
        y/=2;
    }
    return res;
}

int main(){

    int K;
    cin>>K;
    vector<int> D(K),L(K);
    for(int i=0;i<K;i++){
        cin>>D[i]>>L[i];
    }    

    /*
    dp[i][x]=i個めまで考えた際に mod7=xとなる通り数
    */

    vector dp(K+2,vector<ll>(7,0));

    const ll mod7=7;
    const ll inv9mod7=mod_pow(9,7-2,7);//9のmod7での逆元 4になる. 9*4=1 (mod7)

    dp[0][0]=1;
    for(int i=0;i<K;i++){
        for(int k=0;k<mod7;k++){
            //i個目を使わない
            dp[i+1][k]+=dp[i][k];
            dp[i+1][k]%=MOD;


            ll nx=mod_pow(10,L[i],mod7)*k%mod7;//10^L[i]

            ll d=(mod_pow(10,L[i],mod7)-1+mod7)%mod7;// 9999..9(9がL[i]-1個) のmod7 
            d*=inv9mod7%mod7;//これでtは 1111..1(1がL[i]-1個) のmod7 
            d%=mod7;
            ll t = d*D[i]%mod7;
            t%=mod7;

            nx+=t;
            nx%=mod7;

            dp[i+1][nx]+=dp[i][k];
            dp[i+1][nx]%=MOD;

        }
    }

    ll ans=0;

    for(int k=0;k<mod7;k++){
        ans+=k*dp[K][k]%MOD;
        ans%=MOD;
    }
    cout<<ans<<endl;



}
0