目录
A.Aqours
记一个数组dp[u]表示点u到叶子节点距离的最小值。按照顺序对每个叶子节点暴力往上搜,因为题目给出的树保证了深的点一定比浅的大,也就保证了距离近的点一定会先更新。搜的和回溯的时候顺便更新dp数组就行了。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=5e6+7;
int p[N],dp[N],ans,n;
bool is[N];
int dfs(int u,int stp){
if(u==0){ans=-1;return N;}
if(~dp[u]){ans=dp[u]+stp;return dp[u];}
int tmp=dfs(p[u],stp+1);
dp[u]=min(stp,tmp+1);
return dp[u];
}
int main()
{
scanf("%d",&n);
is[1]=true;p[1]=0;dp[1]=-1;
for(int i=2;i<=n;i++){
is[i]=true;dp[i]=-1;
scanf("%d",&p[i]);
is[p[i]]=false;
}
for(int i=1;i<=n;i++){
if(is[i]){
printf("%d ",i);
dfs(i,0);
printf("%d\n",ans);
}
}
}
B.玖凛两开花
看错好久的题。
假设答案是x,也就是说0~x-1的点都一定都要与x~n-1的点相连,注意x~n-1的点可以不与任何点相连,但0~x-1的点必须都相连。这个可以做二分匹配来判断。所以只要二分答案x就行了。
也可以不二分,直接从小到大枚举答案,如果当前点被连接了,就让连接的点重连,当前点也要连。这样就少了一个log的复杂度。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e4+7;
struct Edge{
int v,nxt;
}e[N*2];
int p[N],edn;
void add(int u,int v){
e[++edn]=(Edge){v,p[u]};p[u]=edn;
}
int mt[N];
bool vis[N];
bool dfs(int u,int md){
for(int i=p[u];~i;i=e[i].nxt){
int v=e[i].v;
if(vis[v]||v<md) continue;
vis[v]=true;
if(mt[v]==-1||dfs(mt[v],md)){
mt[v]=u;
return true;
}
}
return false;
}
int main()
{
int n,m,u,v;
scanf("%d%d",&n,&m);
memset(p,-1,sizeof(p));edn=-1;
for(int i=1;i<=m;i++){
scanf("%d%d",&u,&v);
if(u>v) swap(u,v);
add(u,v);
}
memset(mt,-1,sizeof(mt));
int ans=0;
for(int i=0;i<n;i++){
if(mt[i]!=-1){
memset(vis,false,sizeof(vis));
if(!dfs(mt[i],i+1)) break;
}
memset(vis,false,sizeof(vis));
if(!dfs(i,i+1)) break;
ans=i+1;
}
printf("%d\n",ans);
}
C.御坂妹妹
三分可以卡过,正解求上下凸包,留坑。
D.吉良吉影的奇妙计划
暴力打表或者dp
dp[i][j][p][q]表示第i位,有了j个左,第i-1位是p,第i位是q的答案。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=998244353;
ll dp[50][50][5][5];
int main()
{
int n;
scanf("%d",&n);
dp[2][2][1][1]=1;
dp[2][0][2][2]=1;
dp[2][1][0][0]=1;
for(int i=3;i<=n*2;i++){
int hi=min(i,n);
int lo=max(0,i-n);
for(int j=lo;j<=hi;j++){
if(j>=1){
dp[i][j][0][0]=dp[i-1][j-1][2][2];
dp[i][j][0][1]=dp[i-1][j-1][0][0];
dp[i][j][1][1]=(dp[i-1][j-1][1][1]+dp[i-1][j-1][0][1])%mod;
}
if(i-1>=j){
dp[i][j][0][0]=(dp[i][j][0][0]+dp[i-1][j][1][1])%mod;
dp[i][j][0][2]=dp[i-1][j][0][0];
dp[i][j][2][2]=(dp[i-1][j][2][2]+dp[i-1][j][0][2])%mod;
}
}
}
int ans=0;
for(int i=0;i<3;i++){
for(int j=0;j<3;j++){
ans=(ans+dp[n*2][n][i][j])%mod;
}
}
printf("%d\n",ans);
}
E.Souls-like Game
矩阵求区间修改和区间积。
div2好像纯暴力也可以,div1线段树加预处理矩阵的幂次。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=1e4+7;
const ll mod=998244353;
struct rec{
ll p[3][3];
ll val;
void read(){
for(int i=0;i<3;i++) for(int j=0;j<3;j++) scanf("%lld",&p[i][j]);
}
ll getv(){
val=0; for(int i=0;i<3;i++) for(int j=0;j<3;j++) val=(val+p[i][j])%mod;
return val;
}
void operator = (const rec k){
for(int i=0;i<3;i++) for(int j=0;j<3;j++) p[i][j]=k.p[i][j];
}
rec operator * (const rec k)const{
rec res;
for(int i=0;i<3;i++){
for(int j=0;j<3;j++){
res.p[i][j]=0;
for(int t=0;t<3;t++){
res.p[i][j]=(res.p[i][j]+p[i][t]*k.p[t][j]%mod)%mod;
}
}
}
return res;
}
rec operator ^ (int k)const{
rec res,a;
for(int i=0;i<3;i++)
for(int j=0;j<3;j++)
res.p[i][j]=(i==j),a.p[i][j]=p[i][j];
while(k){
if(k&1) res=res*a;
a=a*a;k>>=1;
}
return res;
}
}t[N<<2],lz[N<<2];
void pu(int rt){
t[rt]=t[rt<<1]*t[rt<<1|1];
}
void pd(int rt,int l,int r){
if(lz[rt].val){
int m=l+r>>1;
t[rt<<1]=lz[rt]^(m-l+1);
t[rt<<1|1]=lz[rt]^(r-m);
lz[rt<<1].val=lz[rt<<1|1].val=1;
lz[rt<<1]=lz[rt];lz[rt<<1|1]=lz[rt];
lz[rt].val=0;
}
}
void bd(int rt,int l,int r){
lz[rt].val=0;
if(l==r){
t[rt].read();
return;
}
int m=l+r>>1;
bd(rt<<1,l,m);
bd(rt<<1|1,m+1,r);
pu(rt);
}
void upd(int rt,int l,int r,int L,int R,rec val){
if(L<=l&&r<=R){
t[rt]=val^(r-l+1);
lz[rt]=val;
lz[rt].val=1;
return;
}
pd(rt,l,r);
int m=l+r>>1;
if(L<=m) upd(rt<<1,l,m,L,R,val);
if(m<R) upd(rt<<1|1,m+1,r,L,R,val);
pu(rt);
}
rec que(int rt,int l,int r,int L,int R){
if(L<=l&&r<=R) return t[rt];
pd(rt,l,r);
int m=l+r>>1;
rec res;
for(int i=0;i<3;i++) for(int j=0;j<3;j++) res.p[i][j]=(i==j);
if(L<=m) res=res*que(rt<<1,l,m,L,R);
if(m<R) res=res*que(rt<<1|1,m+1,r,L,R);
return res;
}
int main()
{
int n,m;
scanf("%d%d",&n,&m);n--;
bd(1,1,n);
while(m--){
int t,l,r;
rec tmp;
scanf("%d%d%d",&t,&l,&r);
if(t==1){
tmp.read();
upd(1,1,n,l,r,tmp);
}
else{
printf("%lld\n",que(1,1,n,l,r-1).getv());
}
}
}
G.穗乃果的考试
统计每个1对答案的贡献,队友代码
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=998244353;
char s[2005];
int main()
{
ll n,m;
scanf("%lld%lld",&n,&m);
ll ans=0;
for(ll i=1;i<=n;i++)
{
scanf("%s",s+1);
for(ll j=1;j<=m;j++)
if(s[j]=='1')
ans=(ans+(i*j)%mod*(n-i+1)%mod*(m-j+1)%mod)%mod;
}
printf("%lld\n",ans%mod);
return 0;
}
J.去音乐会
两个范围其实都只有1e3,所以可以n^2暴力,然后就是一个拓展欧几里得解不定方程。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
void egcd(int a,int b,int &x,int &y){
if(!b){x=1;y=0;}
else{egcd(b,a%b,y,x);y-=x*(a/b);}
}
int main()
{
int t,a,b,c,d,x,y;
ll n,ans;
scanf("%d",&t);
while(t--){
scanf("%d%d%d%d%lld",&a,&b,&c,&d,&n);
ans=0;
for(int i=a+1;i<=b;i++){
for(int j=c+1;j<=d;j++){
int aa=-b,bb=d,cc=i-j;
int t=__gcd(aa,bb);
if(cc%t) continue;
aa/=t;bb/=t;cc/=t;
egcd(aa,bb,x,y);
if(bb<0) bb=-bb;
x=((1ll*cc*x)%bb+bb)%bb;
if(1ll*b*x+i<=n) ans+=((n-i)/b-x)/(bb)+1;
}
}
printf("%lld\n",ans);
}
}