QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#260028#7794. 水果茶275307894a80 319ms479924kbC++144.3kb2023-11-21 18:29:532023-11-21 18:29:54

Judging History

你现在查看的是最新测评结果

  • [2023-11-21 18:29:54]
  • 评测
  • 测评结果:80
  • 用时:319ms
  • 内存:479924kb
  • [2023-11-21 18:29:53]
  • 提交

answer

#include<bits/stdc++.h>
#define Gc() getchar() 
#define Me(x,y) memset(x,y,sizeof(x))
#define Mc(x,y) memcpy(x,y,sizeof(x))
#define d(x,y) ((m)*(x-1)+(y))
#define R(n) (rnd()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound
#define fi first
#define se second
#define eb emplace_back
using namespace std;using ll=long long;using db=double;using lb=long db;using ui=unsigned;using ull=unsigned long long;using pii=pair<int,int>;
const int N=1e7+5,M=(1<<22)+5,K=(1<<25)+5,mod=1315105,Mod=mod-1;const db eps=1e-9;const int INF=1e9+7;mt19937 rnd(time(0));

namespace Rd{
	char ch,B0[1<<15],*S=B0,*T=B0;
	#define getc() (S==T&&(T=(S=B0)+fread(B0,1,1<<15,stdin),S==T)?0:*S++)
	inline int read(){
		int aa;
		while(ch=getc(),ch<'0'||ch>'9');aa=ch-'0';
		while(ch=getc(),ch>='0'&&ch<='9')aa=aa*10+ch-'0';
		return aa;
	}
}

int n,m,k,A[N];vector<int> G[N];
struct ljb{int head,h[N];pair<int,int> f[N];void add(int x,int y){f[++head]=make_pair(y,h[x]);h[x]=head;}}S;
int buf[N*8],*f[N],*F=buf;
void init(){
	int seed=Rd::read();
	mt19937_64 rnd(seed);
	int MAXM=1e9;
	for(int i=1;i<=n;i++) A[i]=rnd()%MAXM+1;
}
int Ct;
void con(int x,int y){G[x].emplace_back(y);}
namespace Tarjan{
	int dfn[N],low[N],dh,st[N],sh;
	void dfs(int x,int La){
		dfn[x]=low[x]=++dh;st[++sh]=x;
		for(int p=S.h[x],i=S.f[p].fi;p;p=S.f[p].se,i=S.f[p].fi) if(i^La){
			if(dfn[i]) {low[x]=min(low[x],dfn[i]);continue;}
			dfs(i,x);low[x]=min(low[x],low[i]);
			if(low[i]>=dfn[x]) {++Ct;con(x,Ct);while(st[sh+1]^i) con(Ct,st[sh--]);}
		}
	}
	void build(){
		Ct=n;dfs(1,0);
	}
}
int Le[N],son[N];
void dfs1(int x){
	for(int i:G[x]) dfs1(i);
	if(x<=n){
		for(int i:G[x]) Le[i]>Le[x]&&(son[x]=i,Le[x]=Le[i]);
	}else{
		int cir=G[x].size()+1;
		for(int i=0;i<G[x].size();i++) if(Le[G[x][i]]+min(i+1,cir-i-1)>Le[x]) Le[x]=Le[G[x][i]]+min(i+1,cir-i-1),son[x]=G[x][i];
	}
}
void dfs2(int x){
	f[x]=F;F+=Le[x]-Le[son[x]];if(!son[x]) F++;
	if(son[x]) dfs2(son[x]);
	for(int i:G[x]) if(i^son[x]) dfs2(i);
}
ll ans;
int g[N*6];
void add(int &x,int y){(x+=y)>=mod&&(x-=mod);}
vector<pii> clr;
void dfs3(int x){
	for(int i:G[x]) dfs3(i);
	if(x<=n){
		// cerr<<x<<' '<<Le[x]<<' '<<son[x]<<'\n';
		if(k<=Le[x]) ans+=1ll*A[x]*f[x][k]%mod;
		f[x][0]=A[x]%mod;
		for(int i:G[x]) if(i^son[x]){
			for(int j=0;j<=Le[i];j++) if(j+Le[x]>=k&&j<=k) ans+=1ll*f[i][j]*f[x][k-j]%mod;
			for(int j=0;j<=Le[i];j++) add(f[x][j],f[i][j]);
		}
	}else{
		int cir=G[x].size()+1;
		int id=0;for(int i=0;i<G[x].size();i++) if(G[x][i]==son[x]) id=i;
		for(int i=0;i<G[x].size();i++) if(G[x][i]^son[x]){
			int y=G[x][i],det=min(abs(i-id),cir-abs(i-id));
			for(int j=0;j<=Le[y];j++) if(j+det<=k&&j+det+Le[son[x]]>=k) ans+=1ll*f[y][j]*f[son[x]][k-det-j]%mod;
		}
		clr.clear();
		for(int i=0;i<cir/2;i++) {
			int y=G[x][i];
			if(y^son[x]){
				for(int j=0;j<=Le[y];j++) add(g[j+i],f[y][j]);
				clr.emplace_back(i,i+Le[y]);
			} 
		}
		for(int i=0;i<G[x].size();i++){
			int y=G[x][i];
			if(y^son[x]) {
				for(int j=0;j<=Le[y];j++) add(g[j+i],mod-f[y][j]);
			}
			int pt=(i+cir/2)%cir;
			if(pt!=cir-1&&G[x][pt]^son[x]){
				int z=G[x][pt];
				for(int j=0;j<=Le[z];j++) add(g[j+i+cir/2],f[z][j]);
				clr.emplace_back(i+cir/2,i+cir/2+Le[z]);
			}
			if(y^son[x]){
				for(int j=0;j<=Le[y];j++) if(k-j+i>=0) ans+=1ll*g[k-j+i]*f[y][j]%mod;
			}
		}
		for(auto i:clr) fill(g+i.fi,g+i.se+1,0);
		if(cir%2==0){
			for(int i=0;i<G[x].size();i++)if(G[x][i]^son[x]&&i+cir/2<cir-1){
				int y=G[x][i],z=G[x][i+cir/2];if(z==son[x]) continue;
				int det=cir/2;//cerr<<Le[y]<<' '<<Le[z]<<' '<<det<<'\n';
				if(Le[y]<=Le[z]){
					for(int h=0;h<=Le[y];h++) if(h+det<=k&&h+det+Le[z]>=k) ans-=1ll*f[y][h]*f[z][k-det-h]%mod;
				}else{
					for(int h=0;h<=Le[z];h++) if(h+det<=k&&h+det+Le[y]>=k) ans-=1ll*f[z][h]*f[y][k-det-h]%mod;
				}
			}
		}
		for(int i=0;i<G[x].size();i++) if(G[x][i]^son[x]){
			int y=G[x][i],det=min(i+1,cir-i-1);
			for(int j=0;j<=Le[y];j++) add(f[x][det+j],f[y][j]);
		}
	}
}
int main(){
	int i,j;n=Rd::read();m=Rd::read();k=Rd::read();
	while(m--){
		int x,y;x=Rd::read();y=Rd::read();
		S.add(x,y);S.add(y,x);
	}
	init();
	Tarjan::build();
	dfs1(1);
	dfs2(1);
	dfs3(1);
	printf("%lld\n",ans%mod);
	cerr<<clock()*1.0/CLOCKS_PER_SEC<<'\n';
}

Details

Tip: Click on the bar to expand more detailed information

Subtask #1:

score: 10
Accepted

Test #1:

score: 10
Accepted
time: 52ms
memory: 316844kb

input:

3000 3408 93
1 2
2 3
3 4
4 5
1 6
3 7
6 8
7 9
6 10
6 11
11 12
12 13
9 14
12 15
13 16
15 17
15 18
18 19
15 20
18 21
20 22
18 23
22 24
23 25
23 26
26 27
27 28
27 29
29 30
26 31
27 32
29 33
32 34
30 35
32 36
36 37
37 38
35 39
39 40
37 41
37 42
41 43
39 44
41 45
45 46
45 47
46 48
48 49
47 50
47 51
51 52
...

output:

929234

result:

ok 1 number(s): "929234"

Test #2:

score: 0
Accepted
time: 28ms
memory: 316848kb

input:

3000 3404 8
1 2
1 3
1 4
4 5
4 6
2 7
4 8
5 9
8 10
10 11
9 12
12 13
9 14
11 15
13 16
15 17
16 18
15 19
17 20
18 21
21 22
19 23
21 24
24 25
24 26
25 27
26 28
27 29
27 30
26 31
31 32
28 33
30 34
30 35
32 36
35 37
34 38
34 39
36 40
39 41
39 42
42 43
39 44
44 45
44 46
42 47
45 48
47 49
46 50
46 51
51 52
4...

output:

700402

result:

ok 1 number(s): "700402"

Test #3:

score: 0
Accepted
time: 44ms
memory: 316604kb

input:

2928 3260 15
1 2
1 3
1 4
4 5
2 6
1 7
1 8
2 9
9 10
6 11
8 12
3 13
6 14
11 15
3 16
1 17
16 18
7 19
6 20
1 21
3 22
18 23
4 24
11 25
19 26
21 27
15 28
1 29
18 30
25 31
27 32
18 33
29 34
33 35
16 36
13 37
19 38
13 39
18 40
18 41
34 42
28 43
42 44
4 45
22 46
28 47
36 48
46 49
12 50
39 51
15 52
5 53
45 54
...

output:

563777

result:

ok 1 number(s): "563777"

Test #4:

score: 0
Accepted
time: 27ms
memory: 316572kb

input:

2935 3257 17
1 2
1 3
2 4
4 5
5 6
4 7
7 8
8 9
7 10
10 11
4 12
1 13
4 14
1 15
3 16
7 17
14 18
10 19
15 20
20 21
7 22
4 23
14 24
7 25
22 26
22 27
4 28
11 29
13 30
11 31
23 32
32 33
22 34
31 35
10 36
11 37
25 38
8 39
30 40
6 41
1 42
4 43
18 44
28 45
11 46
11 47
14 48
42 49
36 50
28 51
15 52
16 53
18 54
...

output:

1076031

result:

ok 1 number(s): "1076031"

Test #5:

score: 0
Accepted
time: 24ms
memory: 316504kb

input:

2789 3104 17
1 2
2 3
3 4
4 5
4 6
3 7
2 8
6 9
2 10
6 11
8 12
3 13
5 14
11 15
4 16
12 17
13 18
15 19
6 20
9 21
17 22
20 23
23 24
20 25
17 26
13 27
11 28
2 29
16 30
12 31
21 32
16 33
33 34
14 35
22 36
12 37
37 38
31 39
32 40
12 41
4 42
42 43
18 44
29 45
23 46
19 47
35 48
47 49
34 50
17 51
29 52
13 53
4...

output:

1030647

result:

ok 1 number(s): "1030647"

Test #6:

score: 0
Accepted
time: 28ms
memory: 316624kb

input:

2773 3068 5
1 2
1 3
1 4
3 5
5 6
1 7
6 8
3 9
3 10
10 11
9 12
5 13
9 14
4 15
9 16
11 17
1 18
11 19
17 20
17 21
16 22
4 23
23 24
22 25
12 26
18 27
17 28
21 29
16 30
3 31
29 32
9 33
6 34
9 35
28 36
3 37
31 38
18 39
10 40
23 41
15 42
4 43
12 44
41 45
11 46
3 47
11 48
19 49
25 50
11 51
2 52
24 53
42 54
17...

output:

1019285

result:

ok 1 number(s): "1019285"

Test #7:

score: 0
Accepted
time: 39ms
memory: 316628kb

input:

2756 3058 17
1 2
2 3
3 4
3 5
1 6
6 7
1 8
3 9
2 10
2 11
11 12
12 13
13 14
14 15
5 16
16 17
5 18
10 19
14 20
18 21
17 22
22 23
9 24
7 25
11 26
19 27
12 28
1 29
6 30
4 31
25 32
5 33
23 34
1 35
11 36
1 37
35 38
34 39
32 40
12 41
20 42
12 43
31 44
1 45
34 46
20 47
3 48
44 49
27 50
2 51
35 52
34 53
53 54
...

output:

664550

result:

ok 1 number(s): "664550"

Subtask #2:

score: 10
Accepted

Test #8:

score: 10
Accepted
time: 63ms
memory: 326888kb

input:

99734 99733 32
1 2
1 3
2 4
3 5
4 6
5 7
5 8
5 9
4 10
8 11
10 12
2 13
5 14
14 15
11 16
11 17
1 18
5 19
9 20
16 21
20 22
1 23
23 24
5 25
24 26
25 27
1 28
3 29
15 30
14 31
24 32
3 33
29 34
31 35
11 36
9 37
30 38
4 39
26 40
3 41
25 42
36 43
34 44
43 45
13 46
42 47
5 48
48 49
14 50
30 51
2 52
47 53
53 54
...

output:

1163399

result:

ok 1 number(s): "1163399"

Test #9:

score: 0
Accepted
time: 64ms
memory: 326852kb

input:

99986 99985 77
1 2
2 3
3 4
1 5
3 6
6 7
3 8
6 9
7 10
7 11
9 12
4 13
2 14
6 15
3 16
16 17
15 18
5 19
19 20
7 21
13 22
5 23
22 24
24 25
22 26
9 27
26 28
10 29
18 30
29 31
19 32
17 33
23 34
30 35
11 36
32 37
34 38
27 39
32 40
38 41
18 42
28 43
43 44
2 45
7 46
28 47
10 48
34 49
9 50
18 51
33 52
32 53
3 5...

output:

954238

result:

ok 1 number(s): "954238"

Test #10:

score: 0
Accepted
time: 71ms
memory: 326900kb

input:

99712 99711 27
1 2
2 3
3 4
1 5
4 6
3 7
6 8
8 9
3 10
3 11
10 12
5 13
8 14
4 15
13 16
12 17
6 18
11 19
16 20
14 21
15 22
14 23
23 24
22 25
21 26
25 27
11 28
8 29
21 30
6 31
3 32
25 33
14 34
9 35
15 36
27 37
10 38
24 39
20 40
17 41
4 42
15 43
21 44
16 45
36 46
28 47
12 48
42 49
41 50
30 51
20 52
3 53
2...

output:

328290

result:

ok 1 number(s): "328290"

Test #11:

score: 0
Accepted
time: 56ms
memory: 326752kb

input:

99923 99922 587
1 2
2 3
2 4
1 5
3 6
6 7
3 8
6 9
6 10
2 11
7 12
11 13
7 14
14 15
7 16
8 17
17 18
13 19
17 20
9 21
11 22
11 23
21 24
7 25
25 26
21 27
17 28
5 29
4 30
14 31
1 32
10 33
33 34
7 35
19 36
30 37
15 38
21 39
11 40
39 41
1 42
10 43
33 44
33 45
44 46
44 47
15 48
16 49
40 50
18 51
2 52
17 53
6 ...

output:

271483

result:

ok 1 number(s): "271483"

Test #12:

score: 0
Accepted
time: 64ms
memory: 365448kb

input:

99877 99876 506
1 2
2 3
3 4
4 5
5 6
6 7
7 8
8 9
9 10
10 11
11 12
12 13
13 14
14 15
15 16
16 17
17 18
18 19
19 20
20 21
21 22
22 23
23 24
24 25
25 26
26 27
27 28
28 29
29 30
30 31
31 32
32 33
33 34
34 35
35 36
36 37
37 38
38 39
39 40
40 41
41 42
42 43
43 44
44 45
45 46
46 47
47 48
48 49
49 50
50 51
5...

output:

642658

result:

ok 1 number(s): "642658"

Subtask #3:

score: 10
Accepted

Test #13:

score: 10
Accepted
time: 67ms
memory: 327864kb

input:

100000 100000 229
1 2
2 3
1 4
2 5
1 6
4 7
3 8
3 9
7 10
6 11
2 12
10 13
12 14
5 15
7 16
9 17
16 18
13 19
13 20
15 21
14 22
14 23
16 24
22 25
16 26
17 27
22 28
27 29
28 30
30 31
24 32
29 33
30 34
26 35
32 36
32 37
30 38
37 39
35 40
31 41
32 42
38 43
36 44
35 45
40 46
40 47
42 48
47 49
46 50
48 51
46 5...

output:

910558

result:

ok 1 number(s): "910558"

Subtask #4:

score: 10
Accepted

Test #14:

score: 10
Accepted
time: 47ms
memory: 320204kb

input:

89874 89894 91
1 2
2 3
2 4
1 5
4 6
1 7
5 8
8 9
6 10
7 11
10 12
1 13
10 14
3 15
13 16
11 17
15 18
4 19
7 20
13 21
3 22
15 23
14 24
16 25
14 26
22 27
7 28
1 29
10 30
26 31
28 32
5 33
7 34
15 35
20 36
34 37
16 38
18 39
24 40
13 41
17 42
26 43
18 44
21 45
42 46
41 47
3 48
46 49
34 50
17 51
22 52
10 53
3...

output:

1049004

result:

ok 1 number(s): "1049004"

Subtask #5:

score: 20
Accepted

Test #15:

score: 20
Accepted
time: 55ms
memory: 330620kb

input:

90000 93021 19
1 2
1 3
3 4
1 5
3 6
5 7
5 8
5 9
4 10
5 11
3 12
7 13
10 14
5 15
6 16
9 17
8 18
16 19
17 20
12 21
20 22
22 23
15 24
22 25
22 26
19 27
26 28
27 29
22 30
29 31
23 32
28 33
32 34
26 35
29 36
36 37
28 38
36 39
36 40
34 41
33 42
35 43
41 44
36 45
40 46
43 47
46 48
39 49
44 50
47 51
46 52
52 ...

output:

618387

result:

ok 1 number(s): "618387"

Subtask #6:

score: 20
Accepted

Test #16:

score: 20
Accepted
time: 56ms
memory: 327560kb

input:

89849 95377 39
1 2
2 3
2 4
2 5
4 6
6 7
7 8
4 9
5 10
6 11
2 12
2 13
12 14
1 15
1 16
8 17
12 18
10 19
5 20
9 21
1 22
7 23
23 24
23 25
16 26
26 27
26 28
8 29
25 30
16 31
22 32
14 33
6 34
16 35
11 36
12 37
14 38
13 39
5 40
26 41
1 42
32 43
38 44
31 45
16 46
5 47
29 48
45 49
22 50
50 51
48 52
20 53
33 54...

output:

1122581

result:

ok 1 number(s): "1122581"

Subtask #7:

score: 0
Time Limit Exceeded

Test #17:

score: 20
Accepted
time: 319ms
memory: 479924kb

input:

3000002 4500001 2
2 3
2 1
3 1
4 5
4 1
5 1
6 7
6 1
7 1
8 9
8 1
9 1
10 11
10 1
11 1
12 13
12 1
13 1
14 15
14 1
15 1
16 17
16 1
17 1
18 19
18 1
19 1
20 21
20 1
21 1
22 23
22 1
23 1
24 25
24 1
25 1
26 27
26 1
27 1
28 29
28 1
29 1
30 31
30 1
31 1
32 33
32 1
33 1
34 35
34 1
35 1
36 37
36 1
37 1
38 39
38 1...

output:

527074

result:

ok 1 number(s): "527074"

Test #18:

score: -20
Time Limit Exceeded

input:

4999999 5000000 5847
1 2
2 3
3 4
4 5
5 6
6 7
7 8
8 9
9 10
10 11
11 12
12 13
13 14
14 15
15 16
16 17
17 18
18 19
19 20
20 21
21 22
22 23
23 24
24 25
25 26
26 27
27 28
28 29
29 30
30 31
31 32
32 33
33 34
34 35
35 36
36 37
37 38
38 39
39 40
40 41
41 42
42 43
43 44
44 45
45 46
46 47
47 48
48 49
49 50
50...

output:


result: