P7735

考虑将一条边的贡献记在深度较大的电上,一次修改要改的点太多。

想到如果维护每个点的时间,那么一条边是实边当且仅当两个点的时间相同。

轻重链剖分,对于一次操作,维护每个节点的时间,维护重链上的所有点的点权。

对于一次询问,重边直接把贡献加起来,轻边考虑询问两边的时间,因为轻边个数只有 $\log n$ 个,所以复杂度是对的。

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
#define SegmentTree ST
namespace SegmentTree{
#define ls(k) k<<1
#define rs(k) k<<1|1
struct Node{
int col,val,sum,tcv,tcc,len;
}p[N<<2];
inline void Cc(int k,int co){
p[k].col=co;p[k].tcc=co;
}
inline void Cv(int k,int va){
p[k].val=va;p[k].tcv=va;
p[k].sum=va*p[k].len;
}
inline void PushDown(int k){
if(p[k].tcc){Cc(ls(k),p[k].tcc);Cc(rs(k),p[k].tcc);p[k].tcc=0;}
if(p[k].tcv!=-1){Cv(ls(k),p[k].tcv);Cv(rs(k),p[k].tcv);p[k].tcv=-1;}
}
inline void PushUp(int k){p[k].sum=p[ls(k)].sum+p[rs(k)].sum;}
inline void ChangeV(int k,int l,int r,int z,int y,int x){
if(z<=l&&r<=y){Cv(k,x);return;}PushDown(k);int mid=(l+r)>>1;
if(z<=mid) ChangeV(ls(k),l,mid,z,y,x);
if(mid<y) ChangeV(rs(k),mid+1,r,z,y,x);PushUp(k);
}
inline void ChangeC(int k,int l,int r,int z,int y,int x){
if(z<=l&&r<=y){Cc(k,x);return;}PushDown(k);int mid=(l+r)>>1;
if(z<=mid) ChangeC(ls(k),l,mid,z,y,x);
if(mid<y) ChangeC(rs(k),mid+1,r,z,y,x);PushUp(k);
}
inline void Build(int k,int l,int r){
p[k].tcv=-1;
int mid=(l+r)>>1;p[k].len=r-l+1;if(l==r) return;
Build(ls(k),l,mid);Build(rs(k),mid+1,r);
}
inline int Ask(int k,int l,int r,int z,int y){
if(z<=l&&r<=y) return p[k].sum;PushDown(k);int mid=(l+r)>>1,ans=0;
if(z<=mid) ans+=Ask(ls(k),l,mid,z,y);
if(mid<y) ans+=Ask(rs(k),mid+1,r,z,y);PushUp(k);return ans;
}
inline int Gc(int k,int l,int r,int w){
if(l==r) return p[k].col;PushDown(k);int mid=(l+r)>>1;
if(w<=mid) return Gc(ls(k),l,mid,w);else return Gc(rs(k),mid+1,r,w);PushUp(k);
}
inline void Clear(){mset(p,0);}
}

int t,n,m,id[N],rk[N],tot,fa[N],son[N],siz[N],dep[N],top[N];
vi e[N];

inline void dfs(int k,int fat){
dep[k]=dep[fat]+1;fa[k]=fat;siz[k]=1;
for(int to:e[k]) if(to!=fat){dfs(to,k);siz[k]+=siz[to];if(siz[son[k]]<siz[to]) son[k]=to;}
}
inline void Dfs(int k,int t){
id[k]=++tot;rk[tot]=k;top[k]=t;if(son[k]) Dfs(son[k],t);
for(int to:e[k]) if(to!=son[k]&&to!=fa[k]) Dfs(to,to);
}
inline void Clear(){
ST::Clear();rep(i,1,n) e[i].clear(),id[i]=rk[i]=fa[i]=son[i]=siz[i]=dep[i]=top[i]=0;tot=0;
}
inline void Change(int a,int b,int v){
int A=a,B=b,l=-1;
while(top[a]!=top[b]){
if(dep[top[a]]<dep[top[b]]) swap(a,b);
ST::ChangeC(1,1,n,id[top[a]],id[a],v);
if(top[a]!=a) ST::ChangeV(1,1,n,id[son[top[a]]],id[a],1);
a=fa[top[a]];
if(son[a]) ST::ChangeV(1,1,n,id[son[a]],id[son[a]],0);
}
if(dep[a]<dep[b]) swap(a,b);l=b;
if(l==B){if(son[A]) ST::ChangeV(1,1,n,id[son[A]],id[son[A]],0);}
else if(l==A){if(son[B]) ST::ChangeV(1,1,n,id[son[B]],id[son[B]],0);}
else {if(son[A]) ST::ChangeV(1,1,n,id[son[A]],id[son[A]],0);if(son[B]) ST::ChangeV(1,1,n,id[son[B]],id[son[B]],0);}
// printf("son[A]=%d\n",son[A]);
// printf("son[B]=%d\n",son[B]);
// printf("Ask(5)=%d\n",ST::Ask(1,1,n,id[5],id[5]));
// printf("id[%d]=%d id[%d]=%d\n",b,id[b],a,id[a]);
ST::ChangeC(1,1,n,id[b],id[a],v);
if(a!=b) ST::ChangeV(1,1,n,id[son[b]],id[a],1);
ST::ChangeV(1,1,n,id[b],id[b],0);
// printf("Ask(5)=%d\n",ST::Ask(1,1,n,id[5],id[5]));
}
inline int Ask(int a,int b){
int Ans=0;
// printf("Ask(2)=%d Ask(7)=%d\n",ST::Ask(1,1,n,id[2],id[2]),ST::Ask(1,1,n,id[7],id[7]));
// printf("Ask(5)=%d\n",ST::Ask(1,1,n,id[5],id[5]));
while(top[a]!=top[b]){
if(dep[top[a]]<dep[top[b]]) swap(a,b);
if(a!=top[a]) Ans+=ST::Ask(1,1,n,id[son[top[a]]],id[a]);
int t=top[a];a=fa[top[a]];
int c1=ST::Gc(1,1,n,id[t]),c2=ST::Gc(1,1,n,id[a]);
if(c1==c2&&c1!=0){
// puts("here");
Ans++;
}
}
// printf("Ans=%d\n",Ans);
if(dep[a]<dep[b]) swap(a,b);
if(a!=b) Ans+=ST::Ask(1,1,n,id[son[b]],id[a]);
return Ans;
}

int main(){
// freopen("my.in","r",stdin);
// freopen("my.out","w",stdout);
read(t);
while(t--){
// printf("t=%d\n",t);
read(n);read(m);
rep(i,1,n-1){
int u,v;read(u);read(v);e[u].pb(v);e[v].pb(u);
}
ST::Build(1,1,n);
dfs(1,0);Dfs(1,0);
// rep(i,1,n) printf("son[%d]=%d\n",i,son[i]);
rep(i,1,m){
int op,a,b;read(op);read(a);read(b);
if(op==1){
Change(a,b,i);
}
else{
printf("%d\n",Ask(a,b));
}
}
Clear();
}
return 0;
}

犯的错误:

  • 多测不清空。
  • 跳重链的时候忘记修改儿子。
  • 下传标记覆盖没有注意到 $0$。
  • 讨论初始点儿子的修改出错。

改进:

  • 调试之前先沉住气把代码整体阅读一遍,可以省下很多时间。

P7736

考虑 $k=2$ 不难发现答案就是邻接矩阵的行列式。

这是因为对于路径 $(i,p_i)$,$p$ 是一个排列,而交点个数的奇偶性等同于排列的奇偶性。这提示我们往行列式的方向去想。

考虑 $k>2$,但是 $n_i=n_1$。设 $n_i,n_{i+1}$ 组成的邻接矩阵为 $M_i$,考虑 $\prod |M_i|$ 的结果,不难发现就是问题的答案。这是因为考虑合并两个邻接矩阵 $|M_i||M_{i+1}|$,相当于 $(\sum(-1)^{\mathbb{sgn(p)}})(\sum(-1)^{\mathbb{sgn(q)}})$,正好相当于奇偶性的一个合并。

考虑一般情况,注意到上面的那种情况同样也可以写成 $|\prod M_i|=|M|$,而邻接矩阵的乘法是有意义的,$M_{i,j}$ 就表示第一层图的 $i$ 到第 $k$ 层图的 $j$ 的方案数是多少。考虑一般情况下的 $|M|$,合并两个邻接矩阵,方案的奇偶性不变,所以 $|M|$ 奇偶性的部分是正确的,考虑路径是否可能有交,假设两条路径有交,起点终点分别为 $x_i<x_j,y_i<y_j$,那么合并之后 $(x_i,y_i),(x_j,y_j)$ 是不合法的,但是与之对应的,$(x_i,y_j),(x_j,y_i)$ 是一条交点个数与上面奇偶性不同的一条路径,所以恰好容斥掉了。

所以答案就是相邻邻接矩阵乘积的行列式。

复杂度为 $O(kn^2+n^3\log n)$。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
inline int ksm(int a,int b,int mod){
int res=1;while(b){if(b&1){res=1ll*res*a%mod;}a=1ll*a*a%mod;b>>=1;}return res;
}
inline int inv(int x){return ksm(x,mod-2,mod);}

struct Matrix{
int n,m,a[N][N];
inline void Clear(){mset(a,0);n=m=0;}
inline Matrix operator * (const Matrix &b)const{
assert(m==b.n);
Matrix c;c.Clear();c.n=n;c.m=b.m;
rep(i,1,c.n)rep(j,1,c.m)rep(k,1,m){
c.a[i][j]=(c.a[i][j]+a[i][k]*b.a[k][j]%mod)%mod;
}
return c;
}
inline void Print(){
rep(i,1,n){
rep(j,1,m) printf("%d ",a[i][j]);
puts("");
}puts("");
}
inline int Det(){
assert(n==m);
int x=1;
rep(i,1,n){
if(!a[i][i]){
rep(j,i+1,n) if(a[j][i]){
rep(k,1,n) swap(a[j][k],a[i][k]);
x*=-1;break;
}
}
if(!a[i][i]) return 0;
int nowinv=inv(a[i][i]);
rep(j,i+1,n){
int K=a[j][i]*nowinv%mod;
rep(k,1,n) a[j][k]=((a[j][k]-K*a[i][k]%mod)%mod+mod)%mod;
}
}
// (*this).Print();
int ans=1;
rep(i,1,n) ans=ans*a[i][i]%mod;return (ans*x%mod+mod)%mod;
}
};

int t,n[N],m[N],K;

signed main(){
// freopen("my.in","r",stdin);
// freopen("my.out","w",stdout);
read(t);
while(t--){
read(K);
rep(i,1,K) read(n[i]);rep(i,1,K-1) read(m[i]);
Matrix A;A.Clear();A.n=n[1];A.m=n[1];rep(i,1,A.n) A.a[i][i]=1;
rep(i,1,K-1){
Matrix B;B.Clear();
rep(j,1,m[i]){
int u,v;read(u);read(v);
B.a[u][v]=1;
}
B.n=n[i];B.m=n[i+1];
// B.Print();
A=A*B;
// A.Print();
}
printf("%lld\n",A.Det());
}
return 0;
}

犯的错误:

  • 高斯消元忘记换一个不为 $0$ 的行上来。
  • 空间开小。

P7913

首先想到三分,但是显然是假的。

考虑能否处理出 $f_i$ 表示国内放 $i$ 个桥的最优值是多少,容易想到一个桥对应的那些区间一定是不交的区间,而且如果新加入一个桥,对以前的桥没有影响,我们可以看做是在还没选的区间里新选一些无交的区间,其中没选的区间的第一个区间是一定要选的。

P7915

选一个数后,在另一个也是这个数的位置标号,表示在最后的区间中这个数应该在哪里。比如 1234155324,选了第一个 $1$ 之后再第二个 $1$ 上标 $10$。

首先考虑我们需要保证前 $n$ 个数是一个排列,且对于标为 $k$ 的数往左往右至少有一边小于 $k$。也就是说,对于一种合法方案,设已经标号的区间为 $[z,y]$,假设我们接下来要选的数为 $k$,那么一定要保证 $z-1,y+1$ 这两个位置其中一个为 $k$。

同时,注意到在左右都可以选的情况下先选左边和先选右边不影响情况的合法性。所以贪心即可。

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
int t,a[N],n,st,b[N];
P p[N];
char s[N];

inline bool Check(int l,int r,int z,int y){
// puts("Here");
int cnt=1;
rep(i,1,2*n) b[i]=0;
b[z]=2*n;
while(y-z+1<n){
// printf("l=%d r=%d\n",l,r);
// printf("p[a[l]].se=%d\n",p[a[l]].se);
while(!b[l]&&(p[a[l]].se==z-1||p[a[l]].se==y+1)&&y-z+1<n){
cnt++;
if(p[a[l]].se==z-1) z--,b[z]=2*n-cnt+1;else y++,b[y]=2*n-cnt+1;
s[++st]='L';l++;
// puts("L");
}
if(y-z+1>=n) break;
if(!b[r]&&(p[a[r]].fi==z-1||p[a[r]].fi==y+1)){
cnt++;
if(p[a[r]].fi==z-1) z--,b[z]=2*n-cnt+1;else y++,b[y]=2*n-cnt+1;
s[++st]='R';r--;
// puts("R");
}
else return 0;
}
// rep(i,1,2*n) printf("%d ",b[i]);
int now=n+1;
while(z<=y){
if(b[z]==now) s[++st]='L',z++;
else if(b[y]==now) s[++st]='R',y--;
else assert(0);
now++;
}
return 1;
}
inline void Print(){
rep(i,1,st) putchar(s[i]);puts("");
}

int main(){
freopen("my.in","r",stdin);
freopen("my.out","w",stdout);
read(t);
while(t--){
read(n);
rep(i,1,2*n) read(a[i]);
rep(i,1,2*n) p[a[i]]=mp(-1,-1);
rep(i,1,2*n){
if(p[a[i]].fi==-1){p[a[i]].fi=i;}else p[a[i]].se=i;
}
bool op=0;
st=0;s[++st]='L';
if(Check(2,2*n,p[a[1]].se,p[a[1]].se)){op=1;Print();continue;}
st=0;s[++st]='R';
if(Check(1,2*n-1,p[a[2*n]].fi,p[a[2*n]].fi)){op=1;Print();continue;}
if(!op){puts("-1");}
}
}

P7914

显然是一个区间 DP,设 $f_{l,r}$ 表示一个合法超级括号序列,设 $A$ 是具有 $(C)$ 形式的合法超级括号序列,其中 $C$ 为任意字符串,$S$ 如题意所示,$B$ 是合法括号序列,那么 $f_{l,r}$ 有以下部分组成:

  • $\text{ASB}$
  • $\text{AB}$
  • $\text{A}$

设 $g_{l,r}$ 表示 $AS$ 的个数,$h{l,r}$ 表示 $A$ 的个数。

考虑计算 $h_{l,r}$,不难发现其由以下部分组成:

  • $\text{(B)}$
  • $\text{(S)}$
  • $\text{(SB)}$
  • $\text{(BS)}$

枚举 $S$ 的长度转移即可。

考虑计算 $g_{l,r}$,直接枚举 $S$ 长度即可。

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
int f[N][N],g[N][N][2],h[N][N],n,K;
char s[N];

int main(){
// freopen("my.in","r",stdin);
// freopen("my.out","w",stdout);
read(n);read(K);scanf("%s",s+1);
rep(i,1,n-1){
if((s[i]=='('||s[i]=='?')&&(s[i+1]==')'||s[i+1]=='?')){
f[i][i+1]=h[i][i+1]=1;
}
}
rep(i,3,n){
rep(j,1,n-i+1){
int l=j,r=j+i-1;
// rep(k,1,K){
// if(s[l+k-1]!='*'&&s[l+k-1]!='?') break;
// if(l+k-1>=r-1) break;
// g[l][r][0]=(g[l][r][0]+h[l+k][r])%mod;
// }
rep(k,1,K){
if(s[r-k+1]!='*'&&s[r-k+1]!='?') break;
if(l+1>=r-k+1) break;
g[l][r][1]=(g[l][r][1]+h[l][r-k])%mod;
}
if(s[l]!='*'&&s[l]!=')'&&s[r]!='*'&&s[r]!='('){
int L=l+1,R=r-1;
bool op=1;
rep(k,L,R) if(s[k]!='*'&&s[k]!='?'){op=0;break;}
if(R-L+1<=K&&op) h[l][r]++;
h[l][r]=(h[l][r]+f[L][R])%mod;
rep(k,L,R-1){
if(s[k]!='*'&&s[k]!='?') break;
if(k-L+1>K) break;
h[l][r]=(h[l][r]+f[k+1][R])%mod;
}
dec(k,L+1,R){
if(s[k]!='*'&&s[k]!='?') break;
if(R-k+1>K) break;
h[l][r]=(h[l][r]+f[L][k-1])%mod;
}
}
else h[l][r]=0;
rep(k,l,r-1){
f[l][r]=(f[l][r]+1ll*g[l][k][1]*f[k+1][r]%mod)%mod;
f[l][r]=(f[l][r]+1ll*h[l][k]*f[k+1][r]%mod)%mod;
}
f[l][r]=(f[l][r]+h[l][r])%mod;
// printf("g[l][r][1]=%d h[l][r]=%d f[l][r]=%d l=%d r=%d\n",g[l][r][1],h[l][r],f[l][r],l,r);
}
}
printf("%d\n",f[1][n]);
return 0;
}

P7916

参考 Piwry 的博客

首先显然有一个最小割的模型,但是网络流时间复杂度不够优。考虑网格图为平面图,而 $k=2$ 的时候可以利用平面图转对偶图,考虑如何扩展到一般情况。

利用同样的思路,先建出原图的对偶图,考虑我们只需要考虑那些相邻不用颜色附加点之间的点,把这些点两两配对即可。容易想到路径是不会相交的,如果相交,改变配对的方式,一定不会比原先的方案劣,所以我们只需要建出对偶图,然后跑最短路,最后做一个 DP 去求解最优配对方案即可,做 DP 的时候记得有路径不交的限制。

注意,相邻相同颜色附加点之间的点也是有意义的,可能会出现以下情况:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174

vc<P> e[N*N];
int c[N][N],r[N][N],n,m,T,d[N*N],f[N][N],xu[N*N],cnt[N*N],w[N][N];
bool vis[N*N];
vi v;
P id[N<<2];
struct Point{
int w,id,co;
inline bool operator < (const Point &b)const{return id<b.id;}
}p[N*N];

inline int ID(int a,int b){
return (a-1)*(m-1)+b;
}
struct Node{
int id,v;
inline bool operator < (const Node &b)const{return v>b.v;}
};
priority_queue<Node> q;
inline void dij(int id,int s){
fill(d+1,d+1000*1000+1,INF);q.push({s,0});d[s]=0;mset(vis,0);
while(q.size()){
Node top=q.top();q.pop();
if(vis[top.id]) continue;
vis[top.id]=1;
for(P to:e[top.id]){
if(vis[to.fi]||d[to.fi]<=d[top.id]+to.se) continue;
d[to.fi]=d[top.id]+to.se;
q.push({to.fi,d[to.fi]});
}
}
for(int i=0;i<(int)v.size();i++){
// printf("w[%d][%d=%d\n",i,id,d[v[i]]);
w[id][i]=w[i][id]=d[v[i]];
}
}
inline void Print(int a,int b,int c){
// printf("from=%lld to=%lld w=%lld\n",a,b,c);
cnt[a]++;cnt[b]++;
assert(cnt[a]<=e[a].size());
assert(cnt[b]<=e[b].size());
}

signed main(){
// freopen("my.in","r",stdin);
// freopen("my.out","w",stdout);
read(n);read(m);read(T);
rep(i,1,n-1)rep(j,1,m){
read(c[i][j]);
}
rep(i,1,n)rep(j,1,m-1){
read(r[i][j]);
}
rep(i,1,n-2)rep(j,1,m-1){
// Print(ID(i,j),ID(i+1,j),r[i+1][j]);
e[ID(i,j)].pb(mp(ID(i+1,j),r[i+1][j]));
e[ID(i+1,j)].pb(mp(ID(i,j),r[i+1][j]));
}
// puts("Finish colu");
rep(i,1,n-1)rep(j,1,m-2){
// Print(ID(i,j),ID(i,j+1),c[i][j+1]);
e[ID(i,j)].pb(mp(ID(i,j+1),c[i][j+1]));
e[ID(i,j+1)].pb(mp(ID(i,j),c[i][j+1]));
}
rep(i,1,m-1) id[i]=mp(1,i);id[m]=mp(-1,-1);rep(i,m+1,m+n-1) id[i]=mp(i-m,m-1);id[m+n]=mp(-1,-1);
rep(i,m+n+1,2*m+n-1) id[i]=mp(n-1,2*m+n-i);id[2*m+n]=mp(-1,-1);
rep(i,2*m+n+1,2*m+2*n-1) id[i]=mp(2*m+2*n-i,1);id[2*m+2*n]=mp(-1,-1);
// puts("Finish");
//矩形内部建边。
while(T--){
// printf("T=%d\n",T);
int pc;read(pc);
if(pc==1){puts("0");continue;}
rep(i,1,pc){
int w,id,co;read(w);read(id);read(co);
p[i].w=w;p[i].id=id;p[i].co=co;
}
sort(p+1,p+pc+1);
int tot=ID(n-1,m-1);
v.clear();
rep(i,2,pc){
// printf("i=%lld\n",i);
++tot;
rep(j,p[i-1].id,p[i].id-1){
if(id[j]==mp(-1ll,-1ll)) continue;
int w=-1;
if(j<=m){w=r[id[j].fi][id[j].se];}
else if(j<=m+n){w=c[id[j].fi][id[j].se+1];}
else if(j<=2*m+n){w=r[id[j].fi+1][id[j].se];}
else if(j<=2*m+2*n){w=c[id[j].fi][id[j].se];}
e[tot].pb(mp(ID(id[j].fi,id[j].se),w));
e[ID(id[j].fi,id[j].se)].pb(mp(tot,w));
Print(tot,ID(id[j].fi,id[j].se),w);
}
if(i!=2){
e[tot].pb(mp(tot-1,p[i-1].w));
e[tot-1].pb(mp(tot,p[i-1].w));
Print(tot,tot-1,p[i-1].w);
}
if(p[i].co!=p[i-1].co) v.pb(tot);
}
++tot;
e[tot-1].pb(mp(tot,p[pc].w));
e[tot].pb(mp(tot-1,p[pc].w));
Print(tot-1,tot,p[pc].w);
rep(j,p[pc].id,2*m+2*n){
if(id[j]==mp(-1ll,-1ll)) continue;
int w=-1;
if(j<=m){w=r[id[j].fi][id[j].se];}
else if(j<=m+n){w=c[id[j].fi][id[j].se+1];}
else if(j<=2*m+n){w=r[id[j].fi+1][id[j].se];}
else if(j<=2*m+2*n){w=c[id[j].fi][id[j].se];}
// printf("id[%lld].fi=%lld if[%lld].se=%lld w=%lld\n",j,id[j].fi,j,id[j].se,w);
e[tot].pb(mp(ID(id[j].fi,id[j].se),w));
e[ID(id[j].fi,id[j].se)].pb(mp(tot,w));
Print(tot,ID(id[j].fi,id[j].se),w);
}
// exit(0);
rep(j,1,p[1].id-1){
if(id[j]==mp(-1ll,-1ll)) continue;
int w=-1;
if(j<=m){w=r[id[j].fi][id[j].se];}
else if(j<=m+n){w=c[id[j].fi][id[j].se+1];}
else if(j<=2*m+n){w=r[id[j].fi+1][id[j].se];}
else if(j<=2*m+2*n){w=c[id[j].fi][id[j].se];}
e[tot].pb(mp(ID(id[j].fi,id[j].se),w));
e[ID(id[j].fi,id[j].se)].pb(mp(tot,w));
Print(tot,ID(id[j].fi,id[j].se),w);
}
// exit(0);
if(p[1].co!=p[pc].co) v.pb(tot);
e[tot].pb(mp(ID(n-1,m-1)+1,p[1].w));
e[ID(n-1,m-1)+1].pb(mp(tot,p[1].w));
Print(tot,ID(n-1,m-1)+1,p[1].w);
if(v.empty()){
puts("0");rep(i,1,tot){
while(cnt[i]) e[i].pop_back(),cnt[i]--;
}continue;
}
//初始化完成
for(int i=0;i<(int)v.size()-1;i++) dij(i,v[i]);
//最短路
rep(i,0,(int)v.size()-1) xu[i+1]=i;rep(i,v.size()+1,2*v.size()) xu[i]=xu[i-v.size()];
int len=2*v.size();
// rep(i,1,len/2)rep(j,1,len/2) printf("w[%lld][%lld]=%lld\n",xu[i],xu[j],w[xu[i]][xu[j]]);
rep(i,1,len-1) f[i][i+1]=w[xu[i]][xu[i+1]];
for(int i=4;i<=len;i+=2){
rep(j,1,len-i+1){
int l=j,r=j+i-1;
f[l][r]=f[l+1][r-1]+w[xu[l]][xu[r]];
rep(k,l,r-1){
if((k-l+1)&1) continue;
f[l][r]=min(f[l][r],f[l][k]+f[k+1][r]);
}
// printf("f[%d][%d]=%d\n",l,r,f[l][r]);
}
}
// printf("e[9307].size()=%d\n",e[9307].size());
int ans=INF;
rep(i,1,len/2+1) ans=min(ans,f[i][i+len/2-1]);
printf("%lld\n",ans);
rep(i,1,len) xu[i]=0;
rep(i,1,tot){
// assert(cnt[i]>=0);
while(cnt[i]>0){
if(e[i].empty()) assert(0);
e[i].pop_back();
cnt[i]--;
}
}
// exit(0);
}wx0
}

犯的错误:

  • 建边时下标弄错。
  • 建图出错。
  • 小于等于写成小于。

CF1194G Another Meme Problem

首先考虑可以枚举 $a,b$,然后对 $ax,bx$ 中的 $x$ 来 DP,只需要记录进位,是否存在 $a,b$,以及是否小于 $n$。但是有可能一对 $A,B$ 可能对应多个 $a,b$,所以我们规定 $a,b$ 互质,但是有可能一个 $A,B$,对应的 $a,b$ 并不互质,而 $a,b$ 只有可能是 $a’,b’$ $1$ 到 $4$ 倍,所以我们改一下 DP 状态,记录是否存在 $a,b$ 和它们的 $1$ 至 $4$ 倍即可。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
#include<bits/stdc++.h>
#define mset(a,b) memset((a),(b),sizeof((a)))
#define rep(i,l,r) for(int i=(l);i<=(r);i++)
#define dec(i,l,r) for(int i=(r);i>=(l);i--)
#define cmax(a,b) (((a)<(b))?(a=b):(a))
#define cmin(a,b) (((a)>(b))?(a=b):(a))
#define Next(k) for(int x=head[k];x;x=li[x].next)
#define vc vector
#define ar array
#define pi pair
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define N 110
#define M number
using namespace std;

typedef double dd;
typedef long double ld;
typedef long long ll;
typedef unsigned int uint;
typedef unsigned long long ull;
#define int long long
typedef pair<int,int> P;
typedef vector<int> vi;

const int INF=0x3f3f3f3f;
const dd eps=1e-9;
const int mod=998244353;

template<typename T> inline void read(T &x){
x=0; int f=1;
char c=getchar();
for(;!isdigit(c);c=getchar()) if(c == '-') f=-f;
for(;isdigit(c);c=getchar()) x=x*10+c-'0';
x*=f;
}

int c[N],n,f[N][21][21][4][4][4][4][2][2],a,b;
string s;

inline int cm(int i,int j,int co){
if(i<j) return 0;
else if(i>j) return 1;
else return co;
}
inline int gcd(int a,int b){return b==0?a:gcd(b,a%b);}

inline int dfs(int posi,int axj,int bxj,int s1,int s2,int s3,int s4,int lim1,int lim2){
if(posi==n+1){
if(!axj&&!bxj&&!lim1&&!lim2&&(s1==3||s2==3||s3==3||s4==3)) return 1;return 0;
}
int &now=f[posi][axj][bxj][s1][s2][s3][s4][lim1][lim2];
if(now!=-1) return now;
int ans=0;
rep(i,0,9){
int nowa=(i*a+axj)%10;
int nowaxj=(i*a+axj)/10;
int nowb=(i*b+bxj)%10;
int nowbxj=(i*b+bxj)/10;
int nows1=s1|((nowa==a)<<1)|(nowb==b);
int nows2=s2|((nowa==a*2)<<1)|(nowb==b*2);
int nows3=s3|((nowa==a*3)<<1)|(nowb==b*3);
int nows4=s4|((nowa==a*4)<<1)|(nowb==b*4);
ans=(ans+dfs(posi+1,nowaxj,nowbxj,nows1,nows2,nows3,nows4,cm(nowa,c[posi],lim1),cm(nowb,c[posi],lim2)))%mod;
}
now=ans;
return ans;
}

signed main(){
// freopen("my.in","r",stdin);
// freopen("my.out","w",stdout);
cin>>s;int len=s.length()-1;
dec(i,0,len) c[len-i]=s[i]-'0';n=len;
int ans=0;
int now=1;
rep(i,0,n){
ans=(ans+c[i]*now%mod)%mod;now=now*10%mod;
}
// printf("ans=%d\n",ans);
rep(i,1,9)rep(j,i+1,9){
if(gcd(i,j)!=1) continue;
mset(f,-1);
a=i;
b=j;
// ans=(ans+dfs(0,0,0,0,0,0,0,0,0))%mod;
int nowans=dfs(0,0,0,0,0,0,0,0,0);
// printf("a=%d b=%d ans=%d\n",a,b,nowans);
ans=(ans+nowans*2%mod)%mod;
}
printf("%lld\n",ans);
return 0;
}