#include <bits/stdc++.h>
using namespace std;
#define pb push_back
#define eb emplace_back
#define mp make_pair
#define f first
#define s second
#define all(a) (a).begin(),(a).end()
#define For(i,a,b) for(auto i=(a);i<(b);i++)
#define FOR(i,b) For(i,0,b)
#define Rev(i,a,b) for(auto i=(a);i>(b);i--)
#define REV(i,a) Rev(i,a,-1)
#define sz(a) ((int)(a).size())
using ll=long long;using ld=long double;using uint=unsigned int;using ull=unsigned long long;
using pii=pair<int,int>;using pll=pair<ll,ll>;using pill=pair<int,ll>;using plli=pair<ll,int>;using pdd=pair<double,double>;using pld=pair<ld,ld>;
constexpr const char nl='\n',sp=' ';constexpr const int INT_INF=0x3f3f3f3f;constexpr const ll LL_INF=0x3f3f3f3f3f3f3f3f;
constexpr const double D_INF=numeric_limits<double>::infinity();constexpr const ld LD_INF=numeric_limits<ld>::infinity();constexpr const double EPS=1e-9;
template<typename T1,typename T2>bool feq(const T1&x,const T2&y){return is_floating_point<common_type_t<T1,T2>>::value?abs(x-y)<=EPS:x==y;}
template<typename T1,typename T2>bool flt(const T1&x,const T2&y){return is_floating_point<common_type_t<T1,T2>>::value?x<y-EPS:x<y;}
template<typename T1,typename T2>bool fgt(const T1&x,const T2&y){return is_floating_point<common_type_t<T1,T2>>::value?x>y+EPS:x>y;}
template<typename T>constexpr const T&_min(const T&x,const T&y){return x<y?x:y;}template<typename T>constexpr const T&_max(const T&x,const T&y){return x<y?y:x;}
template<typename T,typename...Ts>constexpr const T&_min(const T&x,const Ts&...xs){return _min(x,_min(xs...));}
template<typename T,typename...Ts>constexpr const T&_max(const T&x,const Ts&...xs){return _max(x,_max(xs...));}
template<typename T,typename...Ts>void MIN(T&x,const Ts&...xs){x=_min(x,xs...);}template<typename T,typename...Ts>void MAX(T&x,const Ts&...xs){x=_max(x,xs...);}
template<typename T>constexpr const T&_clamp(const T&v,const T&lo,const T&hi){return v<lo?lo:hi<v?hi:v;}template<typename T>void CLAMP(T&v,const T&lo,const T&hi){v=_clamp(v,lo,hi);}
template<typename T,typename...Args>unique_ptr<T>_make_unique(Args&&...args){return unique_ptr<T>(new T(forward<Args>(args)...));}
template<typename T,typename...Args>shared_ptr<T>_make_shared(Args&&...args){return shared_ptr<T>(new T(forward<Args>(args)...));}
#define min _min
#define max _max
#define clamp _clamp
#define make_unique _make_unique
#define make_shared _make_shared
template<typename...Ts>using uset=unordered_set<Ts...>;template<typename...Ts>using umap=unordered_map<Ts...>;template<typename...Ts>using pq=priority_queue<Ts...>;
template<typename T>using minpq=pq<T,vector<T>,greater<T>>;template<typename T>using maxpq=pq<T,vector<T>,less<T>>;
template<typename T1,typename T2,typename H1=hash<T1>,typename H2=hash<T2>>struct pair_hash{size_t operator()(const pair<T1,T2>&p)const{return 31*H1{}(p.first)+H2{}(p.second);}};
seed_seq seq {
(uint64_t)chrono::duration_cast<chrono::nanoseconds>(chrono::high_resolution_clock::now().time_since_epoch()).count(),
(uint64_t)__builtin_ia32_rdtsc(),(uint64_t)(uintptr_t)make_unique<char>().get()
};
mt19937 rng(seq);
constexpr int _bufferSize=4096,_maxNumLength=128;
char _inputBuffer[_bufferSize+1],*_inputPtr=_inputBuffer,_outputBuffer[_bufferSize],_c,_sign,*_tempInputBuf=nullptr,_numBuf[_maxNumLength],_tempOutputBuf[_maxNumLength],_fill=' ';
FILE*_input=stdin,*_output=stdout,*_error=stderr;const char*_delimiter=" ";int _cur,_outputPtr=0,_numPtr=0,_precision=6,_width=0,_tempOutputPtr=0,_cnt;ull _precisionBase=1000000;
#define _peekchar() (*_inputPtr?*_inputPtr:(_inputBuffer[fread(_inputPtr=_inputBuffer,1,_bufferSize,_input)]='\0',*_inputPtr))
#define _getchar() (*_inputPtr?*_inputPtr++:(_inputBuffer[fread(_inputPtr=_inputBuffer,1,_bufferSize,_input)]='\0',*_inputPtr++))
#define _hasNext() (*_inputPtr||!feof(_input))
#define _readSignAndNum(x) do{(x)=_getchar();}while((x)<=' ');_sign=(x)=='-';if(_sign)(x)=_getchar();for((x)-='0';(_c=_getchar())>='0';(x)=(x)*10+_c-'0')
#define _readFloatingPoint(x,T) for(T _div=1.0;(_c=_getchar())>='0';(x)+=(_c-'0')/(_div*=10))
#define rc(x) do{do{(x)=_getchar();}while((x)<=' ');}while(0)
#define ri(x) do{_readSignAndNum(x);if(_sign)(x)=-(x);}while(0)
#define rd(x) do{_readSignAndNum(x);if(_c=='.')_readFloatingPoint(x,double);if(_sign)(x)=-(x);}while(0)
#define rld(x) do{_readSignAndNum(x);if(_c=='.')_readFloatingPoint(x,ld);if(_sign)(x)=-(x);}while(0)
#define rcs(x) do{_cur=0;do{_c=_getchar();}while(_c<=' ');do{(x)[_cur++]=_c;}while((_c=_getchar())>' ');(x)[_cur]='\0';}while(0)
#define rs(x) do{if(!_tempInputBuf)assert(0);rcs(_tempInputBuf);(x)=string(_tempInputBuf,_cur);}while(0)
#define rcln(x) do{_cur=0;do{_c=_getchar();}while(_c<=' ');do{(x)[_cur++]=_c;}while((_c=_getchar())>=' ');(x)[_cur]='\0';}while(0)
#define rln(x) do{if(!_tempInputBuf)assert(0);rcln(_tempInputBuf);(x)=string(_tempInputBuf,_cur);}while(0)
#define setLength(x) do{if(_tempInputBuf)delete[](_tempInputBuf);_tempInputBuf=new char[(x)+1];}while(0)
void read(int&x){ri(x);}void read(uint&x){ri(x);}void read(ll&x){ri(x);}void read(ull&x){ri(x);}void read(double&x){rd(x);}void read(ld&x){rld(x);}
void read(char&x){rc(x);}void read(char*x){rcs(x);}void read(string&x){rs(x);}
template<typename T1,typename T2>void read(pair<T1,T2>&x){read(x.first);read(x.second);}template<typename T>void read(complex<T>&x){T _re,_im;read(_re);read(_im);x.real(_re);x.imag(_im);}
bool hasNext(){while(_hasNext()&&_peekchar()<=' ')_getchar();return _hasNext();}
template<typename T,typename...Ts>void read(T&x,Ts&...xs){read(x);read(xs...);}
void setInput(FILE*file){*_inputPtr='\0';_input=file;}void setInput(const char*s){*_inputPtr='\0';_input=fopen(s,"r");}void setInput(const string&s){*_inputPtr='\0';_input=fopen(s.c_str(),"r");}
#define _flush() do{_flushBuf();fflush(_output);}while(0)
#define _flushBuf() (fwrite(_outputBuffer,1,_outputPtr,_output),_outputPtr=0)
#define _putchar(x) (_outputBuffer[_outputPtr==_bufferSize?_flushBuf():_outputPtr]=(x),_outputPtr++)
#define _writeTempBuf(x) (_tempOutputBuf[_tempOutputPtr++]=(x))
#define _writeOutput() for(int _i=0;_i<_tempOutputPtr;_putchar(_tempOutputBuf[_i++]));_tempOutputPtr=0
#define _writeNum(x,T,digits) _cnt=0;for(T _y=(x);_y;_y/=10,_cnt++)_numBuf[_numPtr++]='0'+_y%10;for(;_cnt<digits;_cnt++)_numBuf[_numPtr++]='0';_flushNumBuf();
#define _writeFloatingPoint(x,T) ull _I=(ull)(x);ull _F=((x)-_I)*_precisionBase+(T)(0.5);if(_F>=_precisionBase){_I++;_F=0;}_writeNum(_I,ull,1);_writeTempBuf('.');_writeNum(_F,ull,_precision)
#define _checkFinite(x) if(std::isnan(x)){wcs("NaN");}else if(std::isinf(x)){wcs("Inf");}
#define _flushNumBuf() for(;_numPtr;_writeTempBuf(_numBuf[--_numPtr]))
#define _fillBuf(x) for(int _i=0;_i<(x);_i++)_putchar(_fill)
#define _flushTempBuf() int _tempLen=_tempOutputPtr;_fillBuf(_width-_tempLen);_writeOutput();_fillBuf(-_width-_tempLen)
#define wb(x) do{if(x)_writeTempBuf('1');else _writeTempBuf('0');_flushTempBuf();}while(0)
#define wc(x) do{_writeTempBuf(x); _flushTempBuf();}while(0)
#define wi(x) do{if((x)<0){_writeTempBuf('-');_writeNum(-(x),uint,1);}else{_writeNum(x,uint,1);}_flushTempBuf();}while(0)
#define wll(x) do{if((x)<0){_writeTempBuf('-');_writeNum(-(x),ull,1);}else{_writeNum(x,ull,1);}_flushTempBuf();}while(0)
#define wd(x) do{_checkFinite(x)else if((x)<0){_writeTempBuf('-');_writeFloatingPoint(-(x),double);}else{_writeFloatingPoint(x,double);}_flushTempBuf();}while(0)
#define wld(x) do{_checkFinite(x)else if((x)<0){_writeTempBuf('-');_writeFloatingPoint(-(x),ld);}else{_writeFloatingPoint(x,ld);}_flushTempBuf();}while(0)
#define wcs(x) do{int _slen=strlen(x);_fillBuf(_width-_slen);for(const char*_p=(x);*_p;_putchar(*_p++));_fillBuf(-_width-_slen);}while(0)
#define ws(x) do{_fillBuf(_width-(int)(x).length());for(int _i=0;_i<(int)(x).length();_putchar(x[_i++]));_fillBuf(-_width-(int)(x).length());}while(0)
#define setPrecision(x) do{_precision=(x);_precisionBase=1;for(int _i=0;_i<(x);_i++,_precisionBase*=10);}while(0)
#define setDelimiter(x) do{_delimiter=(x);}while(0)
#define setWidth(x) do{_width=(x);}while(0)
#define setFill(x) do{_fill=(x);}while(0)
void writeDelimiter(){for(const char*_p=_delimiter;*_p;_putchar(*_p++));}
void write(const bool&x){wb(x);}void write(const int&x){wi(x);}void write(const uint&x){wi(x);}void write(const ll&x){wll(x);}void write(const ull&x){wll(x);}
void write(const double&x){wd(x);}void write(const ld&x){wld(x);}void write(const char&x){wc(x);}void write(const char*x){wcs(x);}void write(const string&x){ws(x);}
template<typename T1,typename T2>void write(const pair<T1,T2>&x){write(x.first);writeDelimiter();write(x.second);}
template<typename T>void write(const complex<T>&x){write(x.real());writeDelimiter();write(x.imag());}
template<typename T>void write(const T&x){bool _first=1;for(auto&&_i:x){if(_first)_first=0;else writeDelimiter();write(_i);}}
template<typename T,typename...Ts>void write(const T&x,const Ts&...xs){write(x);writeDelimiter();write(xs...);}
void writeln(){_putchar('\n');}template<typename...Ts>void writeln(const Ts&...xs){write(xs...);_putchar('\n');}
void flush(){_flush();}class IOManager{public:~IOManager(){flush();}};unique_ptr<IOManager>iomanager=make_unique<IOManager>();
void setOutput(FILE*file){flush();_output=file;}void setOutput(const char*s){flush();_output=fopen(s,"w");}void setOutput(const string&s){flush();_output=fopen(s.c_str(),"w");}
template<typename...Ts>void debug(const Ts&...xs){FILE*_temp=_output;setOutput(_error);write(xs...);setOutput(_temp);}
template<typename...Ts>void debugln(const Ts&...xs){FILE*_temp=_output;setOutput(_error);writeln(xs...);setOutput(_temp);}
void setError(FILE*file){flush();_error=file;}void setError(const char*s){flush();_error=fopen(s,"w");}void setError(const string&s){flush();_error=fopen(s.c_str(),"w");}
#define MAXN 100005
int N, M, low[MAXN], pre[MAXN], cnt = 0, SZ[MAXN * 2];
ll ans = 0;
bool vis[MAXN];
vector<int> adj[MAXN], adj2[MAXN * 2];
vector<vector<int>> comp;
stack<pii> s;
void bcc(int v, int prev) {
int children = 0;
pre[v] = low[v] = cnt++;
for (int w : adj[v]) {
if (pre[w] == -1) {
children++;
s.emplace(v, w);
bcc(w, v);
MIN(low[v], low[w]);
if ((prev == v && children > 1) || (prev != v && low[w] >= pre[v])) {
comp.eb();
while (s.top().f != v && s.top().s != w) {
comp.back().pb(s.top().f);
comp.back().pb(s.top().s);
s.pop();
}
comp.back().pb(s.top().f);
comp.back().pb(s.top().s);
s.pop();
}
} else if (pre[w] < pre[v]) {
MIN(low[v], pre[w]);
s.emplace(v, w);
}
}
}
void dfs1(int v, int prev) {
if (v < N) vis[v] = true;
SZ[v] = v < N;
for (int w : adj2[v]) {
if (w == prev) continue;
dfs1(w, v);
SZ[v] += SZ[w];
}
}
void dfs2(int v, int r, int prev) {
if (v < N) for (int w : adj2[v]) ans -= (w == prev ? SZ[v] * ll(SZ[v] - 1) : (SZ[r] - SZ[w]) * ll(SZ[r] - SZ[w] - 1)) * ll(sz(comp[w - N]) - 1);
for (int w : adj2[v]) if (w != prev) dfs2(w, r, v);
}
int main() {
// setInput("in.txt");
// setOutput("out.txt");
// setError("err.txt");
read(N, M);
int a, b;
FOR(i, M) {
read(a, b);
adj[--a].pb(--b);
adj[b].pb(a);
}
FOR(i, N) low[i] = pre[i] = -1;
FOR(i, N) {
if (pre[i] == -1) bcc(i, i);
if (!s.empty()) {
comp.eb();
while (!s.empty()) {
comp.back().pb(s.top().f);
comp.back().pb(s.top().s);
s.pop();
}
}
}
FOR(i, sz(comp)) {
sort(all(comp[i]));
comp[i].erase(unique(all(comp[i])), comp[i].end());
for (int v : comp[i]) {
adj2[v].pb(N + i);
adj2[N + i].pb(v);
}
}
FOR(i, N) vis[i] = false;
FOR(i, N) {
if (vis[i]) continue;
dfs1(i, -1);
ans += SZ[i] * ll(SZ[i] - 1) * ll(SZ[i] - 2);
dfs2(i, i, -1);
}
writeln(ans);
return 0;
}
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
10 ms |
7416 KB |
Output is correct |
| 2 |
Correct |
11 ms |
7464 KB |
Output is correct |
| 3 |
Correct |
11 ms |
7520 KB |
Output is correct |
| 4 |
Correct |
9 ms |
7572 KB |
Output is correct |
| 5 |
Correct |
9 ms |
7576 KB |
Output is correct |
| 6 |
Correct |
8 ms |
7672 KB |
Output is correct |
| 7 |
Correct |
10 ms |
7704 KB |
Output is correct |
| 8 |
Execution timed out |
1121 ms |
815592 KB |
Time limit exceeded |
| 9 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
10 ms |
7416 KB |
Output is correct |
| 2 |
Correct |
11 ms |
7464 KB |
Output is correct |
| 3 |
Correct |
11 ms |
7520 KB |
Output is correct |
| 4 |
Correct |
9 ms |
7572 KB |
Output is correct |
| 5 |
Correct |
9 ms |
7576 KB |
Output is correct |
| 6 |
Correct |
8 ms |
7672 KB |
Output is correct |
| 7 |
Correct |
10 ms |
7704 KB |
Output is correct |
| 8 |
Execution timed out |
1121 ms |
815592 KB |
Time limit exceeded |
| 9 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
150 ms |
815592 KB |
Output is correct |
| 2 |
Correct |
172 ms |
815592 KB |
Output is correct |
| 3 |
Correct |
201 ms |
815592 KB |
Output is correct |
| 4 |
Correct |
172 ms |
815592 KB |
Output is correct |
| 5 |
Correct |
161 ms |
815592 KB |
Output is correct |
| 6 |
Correct |
153 ms |
815592 KB |
Output is correct |
| 7 |
Correct |
153 ms |
815592 KB |
Output is correct |
| 8 |
Correct |
212 ms |
815592 KB |
Output is correct |
| 9 |
Correct |
245 ms |
815592 KB |
Output is correct |
| 10 |
Correct |
189 ms |
815592 KB |
Output is correct |
| 11 |
Correct |
154 ms |
815592 KB |
Output is correct |
| 12 |
Correct |
179 ms |
815592 KB |
Output is correct |
| 13 |
Correct |
149 ms |
815592 KB |
Output is correct |
| 14 |
Correct |
136 ms |
815592 KB |
Output is correct |
| 15 |
Correct |
113 ms |
815592 KB |
Output is correct |
| 16 |
Correct |
100 ms |
815592 KB |
Output is correct |
| 17 |
Correct |
11 ms |
815592 KB |
Output is correct |
| 18 |
Correct |
12 ms |
815592 KB |
Output is correct |
| 19 |
Correct |
12 ms |
815592 KB |
Output is correct |
| 20 |
Correct |
11 ms |
815592 KB |
Output is correct |
| 21 |
Correct |
12 ms |
815592 KB |
Output is correct |
| 22 |
Correct |
13 ms |
815592 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
10 ms |
815592 KB |
Output is correct |
| 2 |
Correct |
9 ms |
815592 KB |
Output is correct |
| 3 |
Correct |
11 ms |
815592 KB |
Output is correct |
| 4 |
Correct |
9 ms |
815592 KB |
Output is correct |
| 5 |
Correct |
10 ms |
815592 KB |
Output is correct |
| 6 |
Correct |
10 ms |
815592 KB |
Output is correct |
| 7 |
Correct |
11 ms |
815592 KB |
Output is correct |
| 8 |
Correct |
11 ms |
815592 KB |
Output is correct |
| 9 |
Correct |
9 ms |
815592 KB |
Output is correct |
| 10 |
Correct |
10 ms |
815592 KB |
Output is correct |
| 11 |
Correct |
9 ms |
815592 KB |
Output is correct |
| 12 |
Correct |
9 ms |
815592 KB |
Output is correct |
| 13 |
Correct |
10 ms |
815592 KB |
Output is correct |
| 14 |
Correct |
11 ms |
815592 KB |
Output is correct |
| 15 |
Correct |
9 ms |
815592 KB |
Output is correct |
| 16 |
Correct |
10 ms |
815592 KB |
Output is correct |
| 17 |
Correct |
9 ms |
815592 KB |
Output is correct |
| 18 |
Correct |
12 ms |
815592 KB |
Output is correct |
| 19 |
Correct |
11 ms |
815592 KB |
Output is correct |
| 20 |
Correct |
9 ms |
815592 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
260 ms |
815592 KB |
Output is correct |
| 2 |
Correct |
191 ms |
815592 KB |
Output is correct |
| 3 |
Correct |
216 ms |
815592 KB |
Output is correct |
| 4 |
Correct |
196 ms |
815592 KB |
Output is correct |
| 5 |
Correct |
194 ms |
815592 KB |
Output is correct |
| 6 |
Correct |
200 ms |
815592 KB |
Output is correct |
| 7 |
Correct |
199 ms |
815592 KB |
Output is correct |
| 8 |
Correct |
215 ms |
815592 KB |
Output is correct |
| 9 |
Correct |
206 ms |
815592 KB |
Output is correct |
| 10 |
Correct |
205 ms |
815592 KB |
Output is correct |
| 11 |
Correct |
231 ms |
815592 KB |
Output is correct |
| 12 |
Correct |
230 ms |
815592 KB |
Output is correct |
| 13 |
Correct |
215 ms |
815592 KB |
Output is correct |
| 14 |
Correct |
278 ms |
815592 KB |
Output is correct |
| 15 |
Correct |
166 ms |
815592 KB |
Output is correct |
| 16 |
Correct |
94 ms |
815592 KB |
Output is correct |
| 17 |
Correct |
137 ms |
815592 KB |
Output is correct |
| 18 |
Correct |
177 ms |
815592 KB |
Output is correct |
| 19 |
Correct |
171 ms |
815592 KB |
Output is correct |
| 20 |
Correct |
175 ms |
815592 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
12 ms |
815592 KB |
Output is correct |
| 2 |
Correct |
11 ms |
815592 KB |
Output is correct |
| 3 |
Execution timed out |
1148 ms |
1049600 KB |
Time limit exceeded |
| 4 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Runtime error |
239 ms |
1049600 KB |
Memory limit exceeded: We have a known bug that the memory usage is measured incorrectly (possibly because of Meltdown/Spectre patch), so your solution may be correct. Please submit again. Sorry for the inconvenience. |
| 2 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
10 ms |
7416 KB |
Output is correct |
| 2 |
Correct |
11 ms |
7464 KB |
Output is correct |
| 3 |
Correct |
11 ms |
7520 KB |
Output is correct |
| 4 |
Correct |
9 ms |
7572 KB |
Output is correct |
| 5 |
Correct |
9 ms |
7576 KB |
Output is correct |
| 6 |
Correct |
8 ms |
7672 KB |
Output is correct |
| 7 |
Correct |
10 ms |
7704 KB |
Output is correct |
| 8 |
Execution timed out |
1121 ms |
815592 KB |
Time limit exceeded |
| 9 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
10 ms |
7416 KB |
Output is correct |
| 2 |
Correct |
11 ms |
7464 KB |
Output is correct |
| 3 |
Correct |
11 ms |
7520 KB |
Output is correct |
| 4 |
Correct |
9 ms |
7572 KB |
Output is correct |
| 5 |
Correct |
9 ms |
7576 KB |
Output is correct |
| 6 |
Correct |
8 ms |
7672 KB |
Output is correct |
| 7 |
Correct |
10 ms |
7704 KB |
Output is correct |
| 8 |
Execution timed out |
1121 ms |
815592 KB |
Time limit exceeded |
| 9 |
Halted |
0 ms |
0 KB |
- |
