BZOJ3498 : PA2009 Cakes
生活随笔
收集整理的這篇文章主要介紹了
BZOJ3498 : PA2009 Cakes
小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
令三元環(huán)(i,j,k)中i>j>k,則每條邊只需要從大點(diǎn)連向小點(diǎn)
設(shè)d[x]表示從x連出的邊的條數(shù)
從1到n枚舉點(diǎn)i,然后枚舉所有與i相連的邊(i,x)(x<i)
如果$d[x]\leq\sqrt{m}$,則依次判斷與x相連的邊(x,y)(y<x)中的y是否與i相連
否則,依次判斷與i相連的邊(i,y)(y<x)中的y是否與x相連
用Hash表支持$O(1)$詢問
時(shí)間復(fù)雜度$O(m\sqrt{m})$
?
#include<cstdio> #include<algorithm> using namespace std; const int N=100010,M=250010,B=(1<<23)-1,BUF=5000000; struct edge{int v;edge*nxt;}epool[M],*ecur=epool,*g[N],*j,*k; struct Edge{int x,y;Edge*nxt;}Epool[M],*Ecur=Epool,*G[B+1],*l; int n,m,i,a[N],d[N],x,y,lim,hash;long long ans;char Buf[BUF],*buf=Buf;pair<int,int>e[M]; inline void read(int&a){for(a=0;*buf<48;buf++);while(*buf>47)a=a*10+*buf++-48;} inline int vis(int x,int y){for(l=G[(x<<8|y)&B];l;l=l->nxt)if(l->x==x&&l->y==y)return 1;return 0;} int main(){fread(Buf,1,BUF,stdin),read(n),read(m);while(lim*lim<m)lim++;for(i=1;i<=n;i++)read(a[i]);for(i=0;i<m;i++){read(x),read(y);if(x<y)swap(x,y);e[i].first=x,e[i].second=y;}for(sort(e,e+m),i=0;i<m;i++){d[x=e[i].first]++;ecur->v=y=e[i].second;ecur->nxt=g[x];g[x]=ecur++;Ecur->x=x;Ecur->y=y;Ecur->nxt=G[hash=(x<<8|y)&B];G[hash]=Ecur++;}for(i=3;i<=n;i++)for(j=g[i];j;j=j->nxt){y=a[i]>a[x=j->v]?a[i]:a[x];if(d[x]<=lim){for(k=g[x];k;k=k->nxt)if(vis(i,k->v))ans+=y>a[k->v]?y:a[k->v];}else for(k=j->nxt;k;k=k->nxt)if(vis(x,k->v))ans+=y>a[k->v]?y:a[k->v];}return printf("%lld",ans),0; }
?
總結(jié)
以上是生活随笔為你收集整理的BZOJ3498 : PA2009 Cakes的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UncaughtExceptionHan
- 下一篇: 《Objective-c》-(OC中含有