結果
| 問題 | No.1222 -101 |
| コンテスト | |
| ユーザー |
sugarrr
|
| 提出日時 | 2020-09-04 22:54:12 |
| 言語 | C++17(clang) (17.0.6 + boost 1.89.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 5,611 bytes |
| 記録 | |
| コンパイル時間 | 2,565 ms |
| コンパイル使用メモリ | 166,340 KB |
| 実行使用メモリ | 19,900 KB |
| 最終ジャッジ日時 | 2024-11-26 20:49:42 |
| 合計ジャッジ時間 | 6,925 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 17 WA * 18 |
ソースコード
//#include <bits/stdc++.h>
#include "bits/stdc++.h"
using namespace std;
typedef long long ll;
/*
#include "boost/multiprecision/cpp_int.hpp"
#include "boost/multiprecision/cpp_dec_float.hpp"
namespace mp = boost::multiprecision;
typedef mp::cpp_int LL;
typedef mp::number<mp::cpp_dec_float<1024>> DD;// 仮数部が1024ビットの浮動小数点数型(TLEしたら小さくする)
*/
typedef long double dd;
#define i_7 (ll)(1E9+7)
//#define i_7 998244353
#define i_5 i_7-2
ll mod(ll a){
ll c=a%i_7;
if(c>=0)return c;
return c+i_7;
}
typedef pair<ll,ll> l_l;
typedef pair<dd,dd> d_d;
ll inf=(ll)1E16;
#define rep(i,l,r) for(ll i=l;i<=r;i++)
#define rrep(i,r,l) for(ll i=r;i>=l;i--)
#define pb push_back
ll max(ll a,ll b){if(a<b)return b;else return a;}
ll min(ll a,ll b){if(a>b)return b;else return a;}
void Max(ll &pos,ll val){pos=max(pos,val);}//Max(dp[n],dp[n-1]);
void Min(ll &pos,ll val){pos=min(pos,val);}
void Add(ll &pos,ll val){pos=mod(pos+val);}
dd EPS=1E-12;
#define fastio ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
#define fi first
#define se second
#define endl "\n" //インタラクティブで消す!!!!!!!!!!!!!!!!!!!!!
#define SORT(v) sort(v.begin(),v.end())
#define ERASE(v) v.erase(unique(v.begin(),v.end()),v.end())
#define POSL(v,x) (lower_bound(v.begin(),v.end(),x)-v.begin())
#define POSU(v,x) (upper_bound(v.begin(),v.end(),x)-v.begin())
template<class T,class S>
inline bool chmax(T &a, S b) {
if(a < b) {
a = b;
return true;
}
return false;
}
template<class T,class S>
inline bool chmin(T &a, S b) {
if(a > b) {
a = b;
return true;
}
return false;
}
void mod_print(ll k){
ll P=50000;
rep(y,1,P){
ll x=mod(y*k);
if(abs(x)<=P||x+P>=i_7){
if(x+P>=i_7){
x-=i_7;
}
cout<<x<<"/"<<y<<endl;
//cout<<setprecision(5)<<(dd)x/(dd)y;
return;
}
}
cout<<"nun"<<endl;
}
#define all(c) c.begin(),c.end()
typedef vector<ll> vl;
//////////////////////////
template<typename T, typename Func_mv>
struct Segtree {
ll n, n_org;
vector<T> dat;
Func_mv merge_values; //dataとdataの演算
T te; //dataの単位元かつ初期値
Segtree(){}
Segtree(ll n_org,
Func_mv merge_values,
T te):
n_org(n_org),
merge_values(merge_values),
te(te){
n = 1;
while(n < n_org) n <<= 1;
dat.resize(2*n-1, te);
}
void build(vector<T>& A){
for(ll k=0; k<ll(A.size()); k++) dat[k+n-1] = A[k];
for(ll k=n-2; k>=0; k--) dat[k] = merge_values(dat[2*k+1], dat[2*k+2]);
}
void update(ll k, T x){ //位置kの要素を更新
k+=n-1;
dat[k]=x;
while(k>0){
k=(k-1)/2;
dat[k]=merge_values(dat[2*k+1],dat[2*k+2]);
}
}
void add(ll k,T x){//位置kの要素に加算
k+=n-1;
dat[k]+=x;
while(k>0){
k=(k-1)/2;
dat[k]=merge_values(dat[2*k+1],dat[2*k+2]);
}
}
T get(ll a, ll b){ //[a,b]から(最小値とか和とかを)get
return query(a, b+1, 0, 0, n);
}
private:
T query(ll a, ll b, ll k, ll lb, ll rb){
if(rb<=a || b<=lb) return te;
if(a<=lb && rb<=b) return dat[k];
ll mb = (lb+rb)>>1;
T vl = query(a, b, 2*k+1, lb, mb);
T vr = query(a, b, 2*k+2, mb, rb);
return merge_values(vl, vr);
}
};
auto make_segtree = [](ll N){
//点更新と区間和
using T=ll;
auto merge_values=[](T&a,T&b){return mod(a+b);};
T te=0;
return Segtree<T, decltype(merge_values)>
(N, merge_values, te);
};
/*
auto st = make_segtree(n);
vector<ll> init(n);
st.build(init);
*/
struct task{
ll l,r,p;
};
ll po(ll i,ll p){
if(p==0)return 1;
else{
i=mod(i);
if(p==1)return i;
if(p%2==0)return po(mod(i*i),p/2);
return mod(i*po(i,p-1));
}
}
ll bunbo(ll n){
return po(n,i_5);
}
int main(){fastio
ll n,m;cin>>n>>m;
vector<task>va,vb;
rep(zz,1,m){
task t;cin>>t.l>>t.r>>t.p;
if(t.p==0)va.pb(t);
else vb.pb(t);
}
ll im[n+5];memset(im,0,sizeof(im));
for(auto t:vb){
im[t.l]++;
im[t.r+1]--;
}
rep(i,1,n+4)im[i]+=im[i-1];
ll ans1,ans2=0;
{
ll cnt=0;
rep(i,1,n)if(im[i]>0)cnt++;
ans1=po(2,cnt-vb.size());
}
{
vector<ll>comp;comp.pb(0);
rep(i,1,n)if(im[i]==0)comp.pb(i);
SORT(comp);ERASE(comp);
for(auto &x:va){
x.l=POSL(comp,x.l);
x.r=POSL(comp,x.r+1)-1;
if(x.l>x.r){
cout<<0<<endl;return 0;
}
}
/*
for(auto x:va){
cout<<x.l<<" "<<x.r<<endl;
}*/
ll cs=comp.size();
ll bnd[cs];memset(bnd,0,sizeof(bnd));
for(auto x:va){
chmax(bnd[x.r],x.l);
}
auto ep = make_segtree(cs);
vector<ll> init(cs);
ep.build(init);
ep.update(0,1);
rep(i,1,cs-1){
ep.update(i,ep.get(bnd[i-1],i-1)*bunbo(2));
}
rep(i,0,cs-1){
ll ad=mod(ep.get(i,i)*po(2,cs-1));
// cout<<i<<":"<<ad<<endl;
if(i==0&&va.size()!=0)continue;
Add(ans2,ad);
}
}
cout<<mod(ans1*ans2)<<endl;
return 0;
}
/*
10 6
2 3 1
6 8 1
7 9 1
1 4 0
3 5 0
1 10 0
*/
sugarrr