結果

問題 No.1261 数字集め
ユーザー chocoruskchocorusk
提出日時 2020-10-21 00:23:56
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 2,502 bytes
コンパイル時間 5,714 ms
コンパイル使用メモリ 191,680 KB
実行使用メモリ 63,544 KB
最終ジャッジ日時 2023-09-28 13:55:13
合計ジャッジ時間 32,841 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 AC 11 ms
4,376 KB
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 AC 33 ms
4,380 KB
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
testcase_38 RE -
testcase_39 RE -
testcase_40 RE -
testcase_41 RE -
testcase_42 RE -
testcase_43 RE -
testcase_44 RE -
testcase_45 RE -
testcase_46 RE -
testcase_47 RE -
testcase_48 RE -
testcase_49 RE -
testcase_50 RE -
testcase_51 RE -
testcase_52 WA -
testcase_53 RE -
testcase_54 RE -
testcase_55 RE -
testcase_56 RE -
testcase_57 RE -
testcase_58 RE -
testcase_59 RE -
testcase_60 RE -
testcase_61 RE -
testcase_62 RE -
testcase_63 RE -
testcase_64 RE -
testcase_65 RE -
testcase_66 RE -
testcase_67 RE -
testcase_68 RE -
testcase_69 RE -
testcase_70 RE -
testcase_71 RE -
testcase_72 RE -
testcase_73 RE -
testcase_74 RE -
testcase_75 RE -
testcase_76 RE -
testcase_77 RE -
testcase_78 AC 55 ms
14,500 KB
testcase_79 RE -
testcase_80 WA -
testcase_81 WA -
testcase_82 RE -
testcase_83 RE -
testcase_84 RE -
testcase_85 RE -
testcase_86 RE -
testcase_87 RE -
testcase_88 RE -
testcase_89 RE -
testcase_90 RE -
testcase_91 RE -
testcase_92 RE -
testcase_93 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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>
#include <list>
#include <atcoder/all>
#define popcount __builtin_popcount
using namespace std;
using namespace atcoder;
typedef long long ll;
typedef pair<int, int> P;
using mint=modint1000000007;
bool e[11000][11000];
int main()
{
    int n, m;
    cin>>n>>m;
    mint a[100010];
    for(int i=2; i<=n; i++){
        int ai; cin>>ai;
        ai--;
        a[i]=mint(ai);
    }
    int q; cin>>q;
    vector<int> v;
    for(int i=2; i<n; i++){
        if(n%i==0) v.push_back(i);
    }
    mint ap[100010];
    for(int i=2; i<=n; i++) ap[i]=a[i].pow(m-1);
    auto calc=[&](int i, int j){
        return (ap[i]*a[n]-ap[i]*a[j]+ap[j]*a[i]-ap[j]*a[n]+ap[n]*a[j]-ap[n]*a[i])/((a[i]-a[n])*(a[i]-a[j])*(a[n]-a[j]));
    };
    mint ans=0;
    for(auto x:v){
        for(auto y:v){
            if(x<y && y%x==0){
                ans+=calc(x, y);
            }
        }
    }
    cout<<ans.val()<<endl;
    v.push_back(1);
    v.push_back(n);
    int x[10010], y[10010];
    for(int i=0; i<q; i++){
        cin>>x[i]>>y[i];
        v.push_back(x[i]);
        v.push_back(y[i]);
    }
    sort(v.begin(), v.end());
    v.erase(unique(v.begin(), v.end()), v.end());
    int l=v.size();
    for(int i=0; i<l; i++){
        for(int j=0; j<i; j++){
            if(v[i]%v[j]==0){
                e[j][i]=1;
            }
        }
    }
    for(int i=0; i<q; i++){
        int xn=lower_bound(v.begin(), v.end(), x[i])-v.begin();
        int yn=lower_bound(v.begin(), v.end(), y[i])-v.begin();
        if(x[i]==1 && y[i]<n){
            for(int j=yn+1; j<l; j++){
                if(e[yn][j] && e[j][l-1]){
                    ans+=calc(y[i], v[j]);
                }
            }
        }else if(x[i]<n && y[i]==n){
            for(int j=1; j<xn; j++){
                if(e[0][j] && e[j][xn]){
                    ans+=calc(v[j], x[i]);
                }
            }
        }else if(1<x[i] && y[i]<n){
            if(e[0][xn] && e[yn][l-1]){
                ans+=calc(x[i], y[i]);
            }
        }
        e[xn][yn]=1;
        cout<<ans.val()<<endl;
    }
    return 0;
}
0