NEERC-2017

A. Archery Tournament

用线段树套set维护横坐标区间内的所有圆,查询时在$O(\log n)$个set中二分查找即可。

时间复杂度$O(n\log^2n)$。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<ctype.h>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); }
#define MS(x, y) memset(x, y, sizeof(x))
#define ls o<<1
#define rs o<<1|1
#define rt 1,1,g
#define mid (l+r>>1)
#define lson ls,l,mid
#define rson rs,mid+1,r
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1, class T2>inline void gmax(T1 &a, T2 b) { if (b > a)a = b; }
template <class T1, class T2>inline void gmin(T1 &a, T2 b) { if (b < a)a = b; }
const int N = 2e6 + 10, M = 0, Z = 1e9 + 7, inf = 0x3f3f3f3f;
template <class T1, class T2>inline void gadd(T1 &a, T2 b) { a = (a + b) % Z; }
int casenum, casei;
int n;
struct AAA
{
    int y, x;
}a[N];
LL K(LL x)
{
    return x * x;
}
bool in_circle(int o1, int o2)
{
    int y1 = a[o1].y;
    int x1 = a[o1].x;
    int y2 = a[o2].y;
    int x2 = a[o2].x;
    return K(y1 - y2) + K(x1 - x2) < K(y2);
}
struct A
{
    int ask;
    int id;
    bool operator < (const A & b)const
    {
        if(ask && b.ask)
        {
            if(a[id].y != a[b.id].y)return a[id].y < a[b.id].y;
            else return id < b.id;
        }
        else if(!ask && !b.ask)
        {
            if(a[id].y != a[b.id].y)return a[id].y < a[b.id].y;
            else return id < b.id;
        }
        else if(ask && !b.ask)
        {
            if(in_circle(id, b.id))return 1;
            else
            {
                if(a[id].y != a[b.id].y)return a[id].y < a[b.id].y;
                else return id < b.id;
            }
        }
        else //if(!ask && b.ask)
        {
            if(in_circle(b.id, id))return 0;
            else
            {
                if(a[id].y != a[b.id].y)return a[id].y < a[b.id].y;
                else return id < b.id;
            }
        }
    }
};
struct Q
{
    int tp, y, x;
}q[N];
int v[N * 3]; int g;
set<A>sot[N * 4];
void build(int o, int l, int r)
{
    sot[o].clear();
    if(l == r)
    {
        return;
    }
    build(lson);
    build(rson);
}
int O, L, R;
void ins(int o, int l, int r)
{
    if(L <= l && r <= R)
    {
        sot[o].insert({0, O});
        return;
    }
    if(L <= mid)ins(lson);
    if(R > mid)ins(rson);
}
void del(int o, int l, int r)
{
    if(L <= l && r <= R)
    {
        sot[o].erase({0, O});
        return;
    }
    if(L <= mid)del(lson);
    if(R > mid)del(rson);
}
int P;
int find(int o, int l, int r)
{
    auto it = sot[o].lower_bound({1, 0});
    if(it != sot[o].end())
    {
        int aim = it->id;
        if(in_circle(0, aim))return aim;
    }
    if(l == r)return 0;
    if(P <= mid)return find(lson);
    else return find(rson);
}
int main()
{
	while(~scanf("%d", &n))
	{
        g = 0;
        for(int i = 1; i <= n; ++i)
        {
            scanf("%d%d%d", &q[i].tp, &q[i].x, &q[i].y);
            if(q[i].tp == 1)
            {
                q[i].tp = i;
                a[i].y = q[i].y; a[i].x = q[i].x;
                v[++g] = q[i].x - q[i].y;
                v[++g] = q[i].x + q[i].y;
            }
            else
            {
                q[i].tp = 0;
                v[++g] = q[i].x;
            }
        }
        sort(v + 1, v + g + 1);
        g = unique(v + 1, v + g + 1) - v - 1;
        build(rt);
        for(int i = 1; i <= n; ++i)
        {
            if(q[i].tp)
            {
                O = q[i].tp;
                L = lower_bound(v + 1, v + g + 1, q[i].x - q[i].y) - v + 1;
                R = lower_bound(v + 1, v + g + 1, q[i].x + q[i].y) - v - 1;
                if(L <= R)
                {
                    ins(rt);
                }
            }
            else
            {
                P = lower_bound(v + 1, v + g + 1, q[i].x) - v;
                a[0].y = q[i].y; a[0].x = q[i].x;
                O = find(rt);
                if(O)
                {
                    printf("%d\n", O);
                    L = lower_bound(v + 1, v + g + 1, q[O].x - q[O].y) - v + 1;
                    R = lower_bound(v + 1, v + g + 1, q[O].x + q[O].y) - v - 1;
                    del(rt);
                }
                else puts("-1");
            }
        }
	}
	return 0;
}
/*
【trick&&吐槽】
4
1 0 12
1 24 10
1 12 3
2 16 14

8
1 0 12
2 -11 22
1 24 10
1 12 3
2 12 12
2 16 14
1 28 15
2 3 6

【题意】


【分析】


【时间复杂度&&优化】


*/

  

B. Box

分类讨论。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<ctype.h>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre() {  }
#define MS(x, y) memset(x, y, sizeof(x))
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1, class T2>inline void gmax(T1 &a, T2 b) { if (b > a)a = b; }
template <class T1, class T2>inline void gmin(T1 &a, T2 b) { if (b < a)a = b; }
const int N = 0, M = 0, Z = 1e9 + 7, inf = 0x3f3f3f3f;
template <class T1, class T2>inline void gadd(T1 &a, T2 b) { a = (a + b) % Z; }
int casenum, casei;
int a[3];
int w, h;
int ww, hh;
bool check4()
{
    return ww <= w && hh <= h;
}
bool check3()
{
    //1
    ww = a[0] * 0 + a[1] * 2 + a[2] * 2;
    hh = a[0] * 1 + a[1] * 0 + a[2] * 2;
    if(check4())return 1;

    //2
    ww = a[0] * 0 + a[1] * 2 + a[2] * 2;
    hh = a[0] * 1 + a[1] * 1 + a[2] * 1;
    if(check4())return 1;

    //3
    ww = a[0] * 1 + a[1] * 1 + a[2] * 3;
    hh = a[0] * 1 + a[1] * 1 + a[2] * 0;
    if(check4())return 1;

    //4
    ww = a[0] * 1 + a[1] * 1 + a[2] * 2;
    hh = a[0] * 1 + a[1] * 1 + a[2] * 1;
    if(check4())return 1;

    //5
    ww = a[0] * 1 + a[1] * 2 + a[2] * 1;
    hh = a[0] * 1 + a[1] * 0 + a[2] * 2;
    if(check4())return 1;

    //6
    ww = a[0] * 2 + a[1] * 1 + a[2] * 1;
    hh = a[0] * 1 + a[1] * 1 + a[2] * 1;
    if(check4())return 1;

    //7
    ww = a[0] * 2 + a[1] * 1 + a[2] * 0;
    hh = a[0] * 1 + a[1] * 1 + a[2] * 2;
    if(check4())return 1;

    return 0;
}
bool check2()
{
    sort(a, a + 3);
    do
    {
        if(check3())return 1;
    }while(next_permutation(a, a + 3));
    return 0;
}
bool check()
{
    if(check2())return 1;
    swap(w, h);
    if(check2())return 1;
    else return 0;
}
int main()
{
	while(~scanf("%d%d%d", &a[0], &a[1], &a[2]))
	{
	    scanf("%d%d", &w, &h);
	    puts(check() ? "Yes" : "No");
	}
	return 0;
}

/*
【trick&&吐槽】


【题意】


【分析】


【时间复杂度&&优化】


*/

  

C. Connections

考虑强连通分量的Kosaraju算法,会发现只有$2(n-1)$条边是有用的,符合题目中保留$2n$条边要求。

#include<cstdio>
const int N=300010;
int Case,n,m,i,x,y,cnt,use[N],g[2][N],nxt[2][N],v[2][N],ed,q[N],t,vis[N],e[N][2];
inline void add(int x,int y){
	v[0][++ed]=y;nxt[0][ed]=g[0][x];g[0][x]=ed;
	v[1][ed]=x;nxt[1][ed]=g[1][y];g[1][y]=ed;
}
void dfs1(int x){
	vis[x]=1;
	for(int i=g[0][x];i;i=nxt[0][i])if(!vis[v[0][i]])use[i]=1,dfs1(v[0][i]);
	q[++t]=x;
}
void dfs2(int x){
	vis[x]=0;
	for(int i=g[1][x];i;i=nxt[1][i])if(vis[v[1][i]])use[i]=1,dfs2(v[1][i]);
}
int main(){
	scanf("%d",&Case);
	while(Case--){
		scanf("%d%d",&n,&m);
		for(ed=0,i=1;i<=n;i++)vis[i]=g[0][i]=g[1][i]=0;
		for(i=1;i<=m;i++){
			scanf("%d%d",&x,&y);
			e[i][0]=x,e[i][1]=y;
			add(x,y);
			use[i]=0;
		}
		for(t=0,i=1;i<=n;i++)if(!vis[i])dfs1(i);
		for(i=n;i;i--)if(vis[q[i]])dfs2(q[i]);
		cnt=n*2;
		for(i=1;i<=m;i++)if(use[i])cnt--;
		for(i=1;i<=m;i++)if(!use[i]&&cnt>0)use[i]=1,cnt--;
		for(i=1;i<=m;i++)if(!use[i])printf("%d %d\n",e[i][0],e[i][1]);
	}
}
/*
111
4 9
1 2
1 3
2 3
2 4
3 2
3 4
4 1
4 2
4 3
*/

  

D. Designing the Toy

斜着构造然后再补充剩下部分即可。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<ctype.h>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre() {  }
#define MS(x, y) memset(x, y, sizeof(x))
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1, class T2>inline void gmax(T1 &a, T2 b) { if (b > a)a = b; }
template <class T1, class T2>inline void gmin(T1 &a, T2 b) { if (b < a)a = b; }
const int N = 105, M = 0, Z = 1e9 + 7, inf = 0x3f3f3f3f;
template <class T1, class T2>inline void gadd(T1 &a, T2 b) { a = (a + b) % Z; }
int casenum, casei;
int A, B, C;
int a[3];
bool v[N][N];
struct AAA
{
    AAA(int x, int y, int z)
    {
        v[0] = x;
        v[1] = y;
        v[2] = z;
    }
    int v[3];
};
vector<AAA>vt;
void SWAP(int x, int y)
{
    int g = vt.size();
    for(int i = 0; i < g; ++i)
    {
        swap(vt[i].v[x], vt[i].v[y]);
    }
}
int main()
{
	while(~scanf("%d%d%d", &A, &B, &C))
	{
        a[0] = A;
        a[1] = B;
        a[2] = C;
        sort(a, a + 3);
        if(a[0] * a[1] < a[2])
        {
            puts("-1");
            continue;
        }
        else
        {
            vt.clear();
            MS(v, 0);
            int cnt = a[2] - a[1];
            for(int i = 1; i <= a[0]; ++i)
            {
                //vt.push_back({i, i, 1});
                vt.push_back(AAA(i, i, 1));
                v[i][i] = 1;
            }
            for(int i = a[0] + 1; i <= a[1]; ++i)
            {
                //vt.push_back({1, i, 1});
                vt.push_back(AAA(1, i, 1));
                v[1][i] = 1;
            }
            for(int i = 1; i <= a[0]; ++i)
            {
                for(int j = 1; j <= a[1]; ++j)
                {
                    if(v[i][j] == 0 && cnt)
                    {
                        --cnt;
                        //vt.push_back({i, j, 1});
                        vt.push_back(AAA(i, j, 1));
                        v[i][j] = 1;
                    }
                }
            }
        }
        if(A >= max(B, C))
        {
            if(C >= B);
            else SWAP(0, 1);
        }
        else if(C >= max(A, B))
        {
            SWAP(0, 2);
            if(A >= B);
            else SWAP(1, 2);
        }
        else
        {
            SWAP(1, 2);
            if(C >= A);
            else SWAP(0, 2);
        }
        printf("%d\n", vt.size());
        for(auto it : vt)
        {
            printf("%d %d %d\n", it.v[0], it.v[1], it.v[2]);
        }
	}
	return 0;
}

/*
【trick&&吐槽】


【题意】


【分析】


【时间复杂度&&优化】


*/

  

E. Easy Quest

依次考虑每个事件,维护每种数字出现次数和待定的数字个数,当数字不够时确定一个数字。

#include<cstdio>
const int N=100000;
int n,i,x,c[N],cnt,m,q[N];
int main(){
  scanf("%d",&n);
  for(i=1;i<=n;i++){
    scanf("%d",&x);
    if(x>0)c[x]++;
    else if(x==0){
      cnt++;
    }else{
      if(c[-x])c[-x]--;
      else if(cnt){
        q[++m]=-x;
        cnt--;
      }else return puts("No"),0;
    }
  }
  puts("Yes");
  for(i=1;i<=m;i++)printf("%d ",q[i]);
  while(cnt--)printf("1 ");
}

  

F. The Final Level

贪心选离目的地较近的那一维坐标靠近即可。

 

G. The Great Wall

二分答案,转化为统计有多少方案的权值和不超过$mid$。

整理出式子后可以通过排序后双指针+树状数组维护。

时间复杂度$O(n\log^2n)$。

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=30010;
int n,len,i;
int m;
ll K,a[N],b[N],c[N],L,R,ANS,MID;
ll f[N],g[N];
int q[N],p[N],bit[N];
inline bool cmp1(int x,int y){return f[x]<f[y];}
inline bool cmp2(int x,int y){return g[x]<g[y];}
inline void add(int x){for(;x<=n;x+=x&-x)bit[x]++;}
inline int ask(int x){int t=0;for(;x;x-=x&-x)t+=bit[x];return t;}
inline int sum(int l,int r){
	l=max(l,1);
	r=min(r,n);
	if(l>r)return 0;
	return ask(r)-ask(l-1);
}
ll cal(){
	int i,j;
	m=n-len+1;
	for(i=1;i<=m;i++){
		f[i]=a[i-1]+b[i+len-1]-b[i-1]-a[i+len-1];
		q[i]=i;
	}
	sort(q+1,q+m+1,cmp1);
	ll ret=0;
	for(i=1;i<=n;i++)bit[i]=0;
	for(i=m,j=1;i;i--){
		while(j<=m&&f[q[i]]+f[q[j]]<=MID)add(q[j++]);
		ret+=sum(q[i]+len,m);
	}
	for(i=1;i<=m;i++){
		f[i]=a[i-1]-b[i-1]+c[i+len-1]-b[i+len-1];
		g[i]=b[i-1]-c[i-1]+b[i+len-1]-a[i+len-1];
		q[i]=p[i]=i;
	}
	sort(q+1,q+m+1,cmp1);
	sort(p+1,p+m+1,cmp2);
	for(i=1;i<=n;i++)bit[i]=0;
	for(i=m,j=1;i;i--){
		while(j<=m&&f[q[i]]+g[p[j]]<=MID)add(p[j++]);
		ret+=sum(q[i]+1,q[i]+len-1);
	}
	return ret;
}
int main(){
	scanf("%d%d%lld",&n,&len,&K);
	for(i=1;i<=n;i++)scanf("%lld",&a[i]),a[i]+=a[i-1];
	for(i=1;i<=n;i++)scanf("%lld",&b[i]),b[i]+=b[i-1];
	for(i=1;i<=n;i++)scanf("%lld",&c[i]),c[i]+=c[i-1];
	L=0,R=c[n];
	while(L<=R){
		MID=(L+R)>>1;
		ll tmp=cal();
		if(tmp>=K)R=(ANS=MID)-1;else L=MID+1;
	}
	printf("%lld",ANS+a[n]);
}
/*
4 2 3
1 2 3 4
3 3 5 5
7 7 7 7
*/

  

H. Hack

留坑。

 

I. Interactive Sort

考虑将偶数划分成若干不可区分大小的集合,将这些集合从小到大排序。一开始所有偶数形成一个不可区分大小的集合。

依次考虑每个奇数,将其在这些集合中二分查找,从对应集合中任取一个数进行比较,得到它大约落在哪两个集合中,再暴力遍历两个集合中所有数字,将对应集合划分成两半。同时,这个过程可以得出有多少偶数小于当前奇数,故可以确定出该奇数。

那么考虑完所有奇数后,所有偶数也就排好了序,因为数据随机,因此可以认为每次都是均分的,比较次数$O(n\log n)$。

#include<cstdio>
#include<vector>
using namespace std;
const int N=11111,K=2;
int n,m,i,j,k;
int cnt;
int cp;
vector<int>a[N];
int fin[N];
int q[N];
//+ -
vector<int>ans;
inline bool cmp(int x,int y){//x<y?
	if(x>0){
		printf("? %d %d\n",x,-y);
		fflush(stdout);
		char op[9];
		scanf("%s",op);
		return op[0]=='<';
	}else{
		printf("? %d %d\n",y,-x);
		fflush(stdout);
		char op[9];
		scanf("%s",op);
		return op[0]=='>';
	}
}
int flag,c0,c1;
int result[N];
int num[N];
inline void gao(int x,int o){
	c0=c1=0;
	int cur=a[q[x]].size();
	for(int i=0;i<cur;i++){
		int y=a[q[x]][i];
		num[i]=y;
		result[i]=cmp(y,o);
		if(result[i]){
			c0++;
		}else{
			c1++;
		}
	}
	if(!c1){
		flag=1;
		return;
	}
	if(!c0){
		flag=-1;
		return;
	}
	flag=0;
	cp++;
	for(int i=cp;i>x+1;i--)q[i]=q[i-1];
	q[x+1]=++cnt;
	a[cnt].clear();
	a[q[x]].clear();
	for(int i=0;i<cur;i++)if(result[i])a[q[x]].push_back(num[i]);
	else a[cnt].push_back(num[i]);
}
inline void gao2(int x,int o){
	int t=cmp(a[q[x]][0],o);
	if(t){
	  flag=1;
	  return;
	}else{
      flag=-1;
      return;
	}
}
int main(){
	scanf("%d",&n);
	if(n==1){
		puts("! 1");
		fflush(stdout);
		return 0;
	}
	if(n==2){
		puts("! 2 1");
		fflush(stdout);
		return 0;
	}
	//even
	//odd
	m=n/2;
	
	
	//use odd to sort even
	cnt=1;
	cp=1;
	a[1].clear();
	for(i=1;i<=m;i++)a[1].push_back(i);
	q[1]=1;
	for(i=1;i<=n-m;i++){
		int l=1,r=cp;
		while(l+K<=r){
			int mid=(l+r)>>1;
			gao2(mid,-i);
			if(flag<0)r=mid;else l=mid;
		}
		while(l<=r){
			int mid=(l+r)>>1;
			gao(mid,-i);
			if(flag==0)break;
			if(flag<0)r=mid-1;else l=mid+1;
		}
		l=max(l-2,1);
		while(l<=cp){
          gao2(l,-i);
          if(flag>0)l++;else break;
		}
		int res=0;
		for(j=1;j<l;j++)res+=a[q[j]].size();
		fin[m+i]=res*2+1;
	}
	int k=0;
	for(i=1;i<=cp;i++){
		for(j=0;j<a[q[i]].size();j++){
			k+=2;
			fin[a[q[i]][j]]=k;
		}
	}
	
	
	printf("!");
	for(i=1;i<=n;i++)printf(" %d",fin[i]);
	puts("");
	fflush(stdout);
}

  

J. Journey from Petersburg to Moscow

若经过至多$k$条边,则答案显然包含在$1$到$n$直接的最短路之中。

否则至少经过$k$条边,枚举最大$k$条边中最小的那一条边的权值$x$,将所有边权值重置为$\max(0,w-x)$,求出$1$到$n$的最短路$d$,则可以用$d+kx$更新答案。这是因为如果此时求出的最短路$d$中边数不够$k$,那么$d+kx$必然不优;而如果收费的边超过$k$条,则减去的$x$过小,也会造成答案不优,因此一定存在一个$x$使得刚好取到最优解。

时间复杂度$O(m^2\log m)$。

#include<cstdio>
#include<queue>
#include<algorithm>
#include<vector>
using namespace std;
typedef long long ll;
typedef pair<ll,int>P;
const int N=3010;
const ll inf=1LL<<60;
int n,m,K,i,g[N],v[N<<1],w[N<<1],nxt[N<<1],ed;
int lim,W;
ll d[N];
ll ans;
priority_queue<P,vector<P>,greater<P> >q;
inline void add(int x,int y,int z){
	v[++ed]=y;
	w[ed]=z;
	nxt[ed]=g[x];
	g[x]=ed;
}
struct E{int x,y,z;}e[N];
inline bool cmp(const E&a,const E&b){return a.z<b.z;}
inline void ext(int x,ll y){
	if(d[x]>y){
		q.push(P(d[x]=y,x));
	}
}
inline void dij(){
	int i,j;
	for(i=1;i<=n;i++)d[i]=inf;
	ext(1,0);
	while(!q.empty()){
		P t=q.top();q.pop();
		if(d[t.second]<t.first)continue;
		for(int i=g[t.second];i;i=nxt[i]){
			ext(v[i],t.first+max(w[i]-W,0));
		}
	}
}
int main(){
	scanf("%d%d%d",&n,&m,&K);
	for(i=1;i<=m;i++)scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].z);
	sort(e+1,e+m+1,cmp);
	ed=1;
	for(i=1;i<=m;i++)add(e[i].x,e[i].y,e[i].z),add(e[i].y,e[i].x,e[i].z);
	dij();
	ans=d[n];
	for(lim=1;lim<=m;lim++){
		W=e[lim].z;
		dij();
		ans=min(ans,d[n]+1LL*W*K);
	}
	printf("%lld",ans);
}

  

K. Knapsack Cryptosystem

留坑。

 

L. Laminar Family

首先将所有未被任何一条树链经过的树边删除,那么剩下的点度数都不能超过$2$,否则必然非法。

那么剩下的图每个连通块都只可能是链,是经典序列问题,把所有区间按左端点排序后单调栈判断即可。

时间复杂度$O(f\log n)$。

#include<cstdio>
#include<algorithm>
#include<cstdlib>
using namespace std;
const int N=100010;
int n,m,i,x,y,g[N],v[N<<1],nxt[N<<1],ed;
int f[N],son[N],size[N],top[N],d[N];
int e[N][2];
int tag[N];
int deg[N];
int G[N],V[N<<1],NXT[N<<1],ED;
bool vis[N];
int ge[N],vv[N],nxte[N],ee;
void NIE(){
	puts("No");
	exit(0);
}
inline void add(int x,int y){v[++ed]=y;nxt[ed]=g[x];g[x]=ed;}
inline void adde(int x,int y){vv[++ee]=y;nxte[ee]=ge[x];ge[x]=ee;}
void dfs(int x){
	d[x]=d[f[x]]+1;
	size[x]=1;
	for(int i=g[x];i;i=nxt[i])if(v[i]!=f[x]){
		f[v[i]]=x;
		dfs(v[i]);
		size[x]+=size[v[i]];
		if(size[v[i]]>size[son[x]])son[x]=v[i];
	}
}
void dfs2(int x,int y){
	top[x]=y;
	if(son[x])dfs2(son[x],y);
	for(int i=g[x];i;i=nxt[i])if(v[i]!=f[x]&&v[i]!=son[x])dfs2(v[i],v[i]);
}
inline void addedge(int x,int y){
	//printf("! %d %d\n",x,y);
	deg[x]++;deg[y]++;
	V[++ED]=y;NXT[ED]=G[x];G[x]=ED;
	V[++ED]=x;NXT[ED]=G[y];G[y]=ED;
}
void dfs3(int x){
	for(int i=g[x];i;i=nxt[i])if(v[i]!=f[x]){
		dfs3(v[i]);
		tag[x]+=tag[v[i]];
	}
	if(x>1&&tag[x])addedge(x,f[x]);
}
inline int lca(int x,int y){
	while(top[x]!=top[y]){
		if(d[top[x]]<d[top[y]])swap(x,y);
		x=f[top[x]];
	}
	return d[x]<d[y]?x:y;
}
int q[N],pos[N];
int all;
struct Seg{
	int l,r;
}seg[N],st[N];
inline bool cmp(const Seg&a,const Seg&b){
	if(a.l!=b.l)return a.l<b.l;
	return a.r>b.r;
}
inline void append(int l,int r){
	if(l>r)swap(l,r);
	seg[++all].l=l;
	seg[all].r=r;
}
inline void check(){
	if(all<=1)return;
	sort(seg+1,seg+all+1,cmp);
	int i;
	int t=0;
	for(i=1;i<=all;i++){
		while(t){
			if(seg[i].r<=st[t].r)break;
			if(seg[i].l<=st[t].r)NIE();
			t--;
		}
		st[++t]=seg[i];
	}
}
inline void solve(int S){
	int cnt=0;
	int i,j;
	while(1){
		vis[S]=1;
		q[++cnt]=S;
		int t=0;
		for(i=G[S];i;i=NXT[i])if(!vis[V[i]]){
			t=V[i];
		}
		if(!t)break;
		S=t;
	}
	//for(i=1;i<=cnt;i++)printf("%d ",q[i]);puts("");
	for(i=1;i<=cnt;i++)pos[q[i]]=i;
	all=0;
	for(i=1;i<=cnt;i++){
		for(j=ge[q[i]];j;j=nxte[j]){
			append(pos[q[i]],pos[vv[j]]);
		}
	}
	check();
}
int main(){
	scanf("%d%d",&n,&m);
	for(i=1;i<n;i++)scanf("%d%d",&x,&y),add(x,y),add(y,x);
	dfs(1);
	dfs2(1,1);
	for(i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		e[i][0]=x,e[i][1]=y;
		int z=lca(x,y);
		tag[x]++;
		tag[y]++;
		tag[z]-=2;
		adde(x,y);
	}
	dfs3(1);
	for(i=1;i<=n;i++)if(deg[i]>2)NIE();
	for(i=1;i<=n;i++)if(deg[i]==1&&!vis[i])solve(i);
	puts("Yes");
}
/*
4 2
1 2
2 3
2 4
1 2
4 2



6 5
1 2
2 3
3 4
5 6
5 2
2 1
6 6
1 4
3 4
4 1
*/

 

posted @ 2017-12-10 02:10  Claris  阅读(2572)  评论(1编辑  收藏  举报