結果

問題 No.1334 Multiply or Add
ユーザー chocoruskchocorusk
提出日時 2021-01-09 20:38:45
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,669 bytes
コンパイル時間 1,403 ms
コンパイル使用メモリ 136,324 KB
実行使用メモリ 8,320 KB
最終ジャッジ日時 2024-04-29 09:17:35
合計ジャッジ時間 7,263 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 WA -
testcase_04 AC 46 ms
8,192 KB
testcase_05 AC 49 ms
8,192 KB
testcase_06 WA -
testcase_07 WA -
testcase_08 AC 49 ms
8,320 KB
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 AC 47 ms
8,192 KB
testcase_13 AC 46 ms
8,192 KB
testcase_14 AC 47 ms
8,064 KB
testcase_15 AC 46 ms
8,192 KB
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 AC 47 ms
8,192 KB
testcase_23 WA -
testcase_24 AC 48 ms
8,192 KB
testcase_25 AC 47 ms
8,320 KB
testcase_26 AC 47 ms
8,192 KB
testcase_27 AC 47 ms
8,192 KB
testcase_28 AC 46 ms
8,192 KB
testcase_29 WA -
testcase_30 AC 47 ms
8,192 KB
testcase_31 AC 47 ms
8,320 KB
testcase_32 WA -
testcase_33 AC 47 ms
8,064 KB
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 AC 46 ms
8,192 KB
testcase_39 AC 2 ms
5,376 KB
testcase_40 AC 20 ms
5,376 KB
testcase_41 AC 8 ms
5,376 KB
testcase_42 AC 9 ms
5,376 KB
testcase_43 AC 36 ms
6,912 KB
testcase_44 AC 21 ms
5,376 KB
testcase_45 AC 8 ms
5,376 KB
testcase_46 AC 22 ms
5,632 KB
testcase_47 AC 31 ms
6,528 KB
testcase_48 AC 22 ms
5,632 KB
testcase_49 AC 42 ms
7,680 KB
testcase_50 AC 82 ms
5,376 KB
testcase_51 AC 82 ms
5,376 KB
testcase_52 AC 82 ms
5,376 KB
testcase_53 AC 81 ms
5,376 KB
testcase_54 AC 80 ms
5,376 KB
testcase_55 AC 89 ms
5,376 KB
testcase_56 AC 2 ms
5,376 KB
testcase_57 AC 2 ms
5,376 KB
testcase_58 AC 2 ms
5,376 KB
testcase_59 WA -
testcase_60 AC 2 ms
5,376 KB
testcase_61 WA -
testcase_62 WA -
testcase_63 AC 2 ms
5,376 KB
testcase_64 AC 2 ms
5,376 KB
testcase_65 AC 2 ms
5,376 KB
testcase_66 AC 69 ms
5,376 KB
testcase_67 AC 44 ms
5,376 KB
testcase_68 AC 4 ms
5,376 KB
testcase_69 AC 75 ms
5,376 KB
testcase_70 AC 71 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

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 <stack>
#include <array>
#include <list>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
const ll MOD=1e9+7;
template<typename T, bool ismin=true>
struct ConvexHullTrickMonotone{
	using P=pair<T, T>;
	deque<P> deq;
	ConvexHullTrickMonotone(){}
	void add(T a, T b){
		if(!ismin) a*=(-1), b*=(-1);
		if(deq.empty()){
			deq.push_back({a, b});
			return;
		}
		if(deq.back().first<=a){
			if(deq.back().first==a && deq.back().second<=b) return;
			while(deq.size()>1){
				T a1=deq.back().first, b1=deq.back().second;
				T a2=deq[deq.size()-2].first, b2=deq[deq.size()-2].second;
				if((b1-b2)*a2*(a-a1)<(b-b1)*a*(a1-a2)) break;
				deq.pop_back();
			}
			deq.push_back({a, b});
		}else{
			while(deq.size()>1){
				T a1=deq.front().first, b1=deq.front().second;
				T a2=deq[1].first, b2=deq[1].second;
				if((b1-b2)*a2*(a-a1)>(b-b1)*a*(a1-a2)) break;
				deq.pop_front();
			}
			deq.push_front({a, b});
		}
	}
	T get(T x, bool isinc){
		if(isinc){
			while(deq.size()>1){
				T a1=deq.front().first, b1=deq.front().second;
				T a2=deq[1].first, b2=deq[1].second;
				if(x/a1+b1<x/a2+b2) break;
				deq.pop_front();
			}
			T ret=x/deq.front().first+deq.front().second;
			if(ismin) return ret;
			else return -ret;
		}else{
			while(deq.size()>1){
				T a1=deq.back().first, b1=deq.back().second;
				T a2=deq[deq.size()-2].first, b2=deq[deq.size()-2].second;
				if(x/a1+b1<x/a2+b2) break;
				deq.pop_back();
			}
			T ret=x/deq.back().first+deq.back().second;
			if(ismin) return ret;
			else return -ret;
		}
	}
	void clear(){
		deq.clear();
	}
};
int main()
{
    int n; cin>>n;
    ll a[200020];
    ll tot=1;
    for(int i=0; i<n; i++){
        cin>>a[i];
        tot=min(MOD, tot*a[i]);
    }
    if(tot==MOD){
        ll ans=1;
        for(int i=0; i<n; i++){
            ans=ans*a[i]%MOD;
        }
        cout<<ans<<endl;
        return 0;
    }
    ll p[200020];
    p[0]=1;
    for(int i=0; i<n; i++){
        p[i+1]=p[i]*a[i];
    }
    using lll=__int128_t;
    ConvexHullTrickMonotone<lll, false> cht;
    ll dp[200020];
    dp[0]=0;
    cht.add(1, 0);
    for(int i=0; i<n; i++){
        dp[i+1]=cht.get(p[i+1], true);
        cht.add(p[i+1], dp[i+1]);
    }
    cout<<dp[n]<<endl;
    return 0;
}
0