玩具谜题 模拟
按照题意模拟就好
#include<cstdio>
struct node{
int fa;
char name[20];
}q[100100];
int main(){
freopen("testdata.in","r",stdin);
int n,m,i,f,num,h=0;
scanf("%d%d",&n,&m);
for(i=0;i<n;i++)scanf("%d%s",&q[i].fa,q[i].name);
//i从1开始 PAC
for(i=0;i<m;i++){
scanf("%d%d",&f,&num);
if(q[h].fa==f)h=(h-num+n)%n;
else h=(h+num)%n;
}
printf("%s",q[h].name);
return 0;
}
总结
简单模拟
天天爱跑步 LCA
学习了大佬题解,讲的非常清晰。
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;
#define PB(v) push_back(v)
template<class T>inline void read(T&x)
{
x=0;int f=0;char ch=getchar();
while(ch<'0'||ch>'9')f|=ch=='-',ch=getchar();
while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
if(f)x=-x;
}
const int N=3e5+10;
const int M=6e5+10;
int n,m,hd[N],tot,w[N],bucket[N+M],ans[N],deep[N],f[N][25],dis[N],spn[M];
bool vis[N];
vector<int>v1[M],v2[M],v3[M];
struct Edge{
int v,nx;
}e[M];
struct node{
int u,v,dis,lca;
}p[N];
inline void add(int u,int v)
{
e[++tot].v=v;
e[tot].nx=hd[u];
hd[u]=tot;
}
inline void dfs(int now,int dep)
{
vis[now]=true;
deep[now]=dep;
for(int i=1;i<=22;i++)
f[now][i]=f[f[now][i-1]][i-1];
for(int i=hd[now];i;i=e[i].nx)
if(!vis[e[i].v])
{
f[e[i].v][0]=now;
dis[e[i].v]=dis[now]+1;
dfs(e[i].v,dep+1);
}
vis[now]=false;
}
inline int LCA(int u,int v)
{
if(deep[u]<deep[v])swap(u,v);
for(int i=0,div=deep[u]-deep[v];i<=22;i++)
if(div&(1<<i))u=f[u][i];
for(int i=22;i>=0;i--)
if(f[u][i]!=f[v][i])u=f[u][i],v=f[v][i];
return u==v?u:f[u][0];
}
//u->lca(u,v)的路径,当且仅当deep[i]+w[i]==deep[u]时,u节点对i节点是有贡献的
//只要符合deep[i]+w[i]的全部是玩家起点的点,就能对i点产生贡献
inline void dfs1(int now)
{
vis[now]=true;
int prev=bucket[deep[now]+w[now]+N];//全部+N避免负数
for(int i=hd[now];i;i=e[i].nx)
if(!vis[e[i].v])dfs1(e[i].v);
bucket[deep[now]+N]+=spn[now];
ans[now]+=bucket[deep[now]+w[now]+N]-prev;
int len=v1[now].size();
for(int k=0;k<len;k++)
--bucket[deep[v1[now][k]]+N];
vis[now]=false;
}
//lca(u,v)->v的路径,当且仅当dis(u,v)-deep[v]=w[i]-deep[i]成立时有贡献
inline void dfs2(int now)
{
vis[now]=true;
int prev=bucket[w[now]-deep[now]+N];
for(int i=hd[now];i;i=e[i].nx)
if(!vis[e[i].v])dfs2(e[i].v);
int len=v2[now].size();
for(int k=0;k<len;k++)
++bucket[v2[now][k]+N];
ans[now]+=bucket[w[now]-deep[now]+N]-prev;
len=v3[now].size();
for(int k=0;k<len;k++)
--bucket[v3[now][k]+N];
vis[now]=false;
}
int main()
{
//freopen("in.txt","r",stdin);
read(n);read(m);
for(int i=1,u,v;i<n;i++)
read(u),read(v),add(u,v),add(v,u);
for(int i=1;i<=n;i++)
read(w[i]);
f[1][0]=1;dfs(1,0);
for(int i=1,u,v;i<=m;i++)
{
read(u);read(v);
p[i].u=u;p[i].v=v;
p[i].lca=LCA(u,v);
p[i].dis=dis[u]+dis[v]-dis[p[i].lca]*2;
spn[u]++;//以x为路径起点的路径条数
v1[p[i].lca].PB(u);//以x为lca的路径的起点的集合
v2[v].PB(p[i].dis-deep[p[i].v]);//以x为终点的路径的起点集合
v3[p[i].lca].PB(p[i].dis-deep[p[i].v]);//以x为LCA的路径的终点的集合
}
dfs1(1);//从下至上
dfs2(1);//从上至下
for(int i=1;i<=m;i++)//如果lca能观察到这条路径上的人,需要去重
if(deep[p[i].u]==deep[p[i].lca]+w[p[i].lca])ans[p[i].lca]--;
for(int i=1;i<=n;i++)
printf("%d ",ans[i]);
printf("\n");
return 0;
}
总结
难度挺大的题
换教室 期望DP+最短路
好久以前做的了
#include<cstdio>
#include<cstring>
#include<algorithm>
#define INF 0x3f3f3f3f
using namespace std;
const int N1=2010;
const int N2=310;
double dp[N1][N1][2];
//dp[i][j][0]表示前i个时间段用j次申请(第i个不申请)所需要的最小体力值
//dp[i][j][1]表示申请
double k[N1];//通过概率
int c[N1];//安排教室
int d[N1];//另外教室
int dis[N1][N1];//边权
int main()
{
//freopen("in.txt","r",stdin);
int n,m,v,e,a,b;
scanf("%d%d%d%d",&n,&m,&v,&e);//时间,申请,教室,道路
int i,j,l,w;//边权
for(i=1;i<=n;i++)
scanf("%d",&c[i]);//安排
for(i=1;i<=n;i++)
scanf("%d",&d[i]);//另外
for(i=1;i<=n;i++)
scanf("%lf",&k[i]);//概率
memset(dis,127/2,sizeof(dis));
for(i=1;i<=e;i++)
{
scanf("%d%d%d",&a,&b,&w);
dis[a][b]=min(dis[a][b],w);
dis[b][a]=min(dis[b][a],w);
}
for(i=1;i<=v;i++)dis[i][i]=0;
for(l=1;l<=v;l++)
for(i=1;i<=v;i++)
for(j=1;j<=v;j++)
if(dis[i][j]>dis[i][l]+dis[l][j])//处理最短路
dis[i][j]=dis[i][l]+dis[l][j];
for(i=1;i<=v;i++)dis[0][i]=0;
for(i=0;i<=2009;i++)
for(j=0;j<=2009;j++)
for(l=0;l<=1;l++)
dp[i][j][l]=INF;
//memset(dp,0x3f,sizeof(dp));
dp[0][0][0]=0;
for(i=0;i<=n-1;i++)
for(j=0;j<=m;j++)
{
dp[i+1][j][0]=min(dp[i][j][0]+dis[c[i]][c[i+1]],dp[i][j][1]+dis[d[i]][c[i+1]]*k[i]+dis[c[i]][c[i+1]]*(1.0-k[i]));
if(j!=m)dp[i+1][j+1][1]=min(dp[i][j][0]+dis[c[i]][d[i+1]]*k[i+1]+dis[c[i]][c[i+1]]*(1.0-k[i+1]),dp[i][j][1]+dis[c[i]][c[i+1]]*(1.0-k[i+1])*(1.0-k[i])+dis[d[i]][c[i+1]]*(1.0-k[i+1])*k[i]+dis[c[i]][d[i+1]]*k[i+1]*(1.0-k[i])+dis[d[i]][d[i+1]]*k[i]*k[i+1]);
}
double ans=INF;
for(j=0;j<=m;j++)
{
ans=min(ans,min(dp[n][j][0],dp[n][j][1]));
}
printf("%.2lf\n",ans);
return 0;
}
总结
无
组合数问题 组合计数+前缀和
线性递推组合数,然后求一个前缀和。
对着大佬题解改改改一直没对,然后发现m可以大于n……
#include<cstdio>
const int N=2e3+10;
int t,k,n,m;
long long C[N][N],sum[N][N];
int main()
{
//freopen("in.txt","r",stdin);
scanf("%d%d",&t,&k);
for(int i=0;i<=2000;i++)C[i][0]=1;
C[1][1]=1;
for(int i=2;i<=2000;i++)
for(int j=1;j<=i;j++)
C[i][j]=(C[i-1][j]+C[i-1][j-1])%k;
for(int i=0;i<=2000;i++)
{
for(int j=1;j<=i;j++)
{
sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];
if(C[i][j]==0)sum[i][j]++;
}
sum[i][i+1]=sum[i][i];
}
while(t--)
{
scanf("%d%d",&n,&m);if(m>n)m=n;
printf("%lld\n",sum[n][m]);
}
return 0;
}
总结
考察组合数的递推
蚯蚓 优先队列+模拟
最开始敲了个85的模拟……想想也没那么简单。
学习了大佬题解,其中单调性很关键。
#include<cstdio>
#include<queue>
#include<cmath>
#include<algorithm>
#define re register
using namespace std;
template<class T>inline void read(T&x)
{
x=0;re int f=0;re char ch=getchar();
while(ch<'0'||ch>'9')f|=ch=='-',ch=getchar();
while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
if(f)x=-x;
}
template<class M>inline void write(M x)
{
re int buf[40],p=0;
if(x<0)putchar('-'),x=-x;
do{
buf[p++]=x%10;x/=10;
}while(x);
for(re int i=p-1;i>=0;i--)putchar(buf[i]+48);
putchar(' ');
}
const int N=7e6+10;
int n,m,q,u,v,t,sigma,t0,t1,t2,h0,h1,h2,cut1[N],cut2[N],now[N];
double p;
priority_queue<int>ans;
inline bool cmp(const int &a,const int &b){return a>b;}
int main()
{
//freopen("in.txt","r",stdin);
read(n);read(m);read(q);read(u);read(v);read(t);
p=(double)u/v;
for(re int i=1;i<=n;i++)
read(now[i]);
t0=n;t1=t2=0;h0=h1=h2=1;
sort(now+1,now+n+1,cmp);
re int top;
for(re int i=1,a1,a2;i<=m;i++)
{
if(h0>t0)
{
if(cut1[h1]>cut2[h2])top=cut1[h1++];
else top=cut2[h2++];
}
else if(now[h0]>=cut1[h1]&&now[h0]>=cut2[h2])top=now[h0++];
else if(now[h0]<=cut1[h1]&&cut2[h2]<=cut1[h1])top=cut1[h1++];
else top=cut2[h2++];
top+=sigma;
a1=floor(p*(double)top);a2=top-a1;
a1-=sigma+q;a2-=sigma+q;
cut1[++t1]=a1;cut2[++t2]=a2;
if(i%t==0)write(top);sigma+=q;
}
putchar('\n');
for(re int i=h0;i<=t0;i++)ans.push(now[i]);
for(re int i=h1;i<=t1;i++)ans.push(cut1[i]);
for(re int i=h2;i<=t2;i++)ans.push(cut2[i]);
for(re int i=1;ans.size();i++)
{
if(i%t==0)write(ans.top()+sigma);
ans.pop();
}
return 0;
}
总结
挺好的一个大模拟。
愤怒的小鸟 状压DP
dp[i]表示i状态时所需要的最少的小鸟数,state[i]表示第i条抛物线所打掉的小猪状态,
dp[i|state[j]]=min(dp[i|state[j]],dp[i]+1)
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define eps 1e-6
int n,m;
const int N=18;
int dp[1<<N];//dp[i]表示i状态时所需最少小鸟数
double x[N],y[N];
int state[510];//记录每条抛物线
int cnt;
int main()
{
//freopen("in.txt","r",stdin);
int t;
double a,b;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&n,&m);
cnt=0;
memset(dp,0x3f,sizeof(dp));
dp[0]=0;
for(int i=0;i<n;i++)
scanf("%lf%lf",&x[i],&y[i]);
for(int i=0;i<n;i++)
{
state[cnt++]=(1<<i);
for(int j=i+1,k=0;j<n;j++)
{
if((k>>j)&1)continue;
a=(y[j]*x[i]-y[i]*x[j])/(x[i]*x[j]*(x[j]-x[i]));
b=(y[i]*x[j]*x[j]-y[j]*x[i]*x[i])/(x[i]*x[j]*(x[j]-x[i]));
if(a>=-eps)continue;
state[cnt]=(1<<i);
for(int o=j;o<n;o++)
if(fabs(a*x[o]*x[o]+b*x[o]-y[o])<=eps)
{
k|=(1<<o);
state[cnt]|=(1<<o);
}
cnt++;
}
}
for(int i=0;i<(1<<n);i++)
for(int j=0;j<cnt;j++)
dp[i|state[j]]=min(dp[i|state[j]],dp[i]+1);
printf("%d\n",dp[(1<<n)-1]);
}
return 0;
}