Link:

Codeforces #174 传送门

A:

求原根的个数,有一条性质是原根个数为$\phi(\phi(n))$,多了一个不会证的性质

如果要确定哪些是原根的话还是要枚举,不过对于每个数不用枚举$p$次了

由于$\delta_p(x) | \phi(x)$,只要对欧拉函数值的约数枚举即可

不过此题好像直接$O(p^2)$枚举就行了……

#include <bits/stdc++.h>

using namespace std;
#define X first
#define Y second
typedef long long ll;
typedef pair<int,int> P;
typedef double db;
int p,res; int main()
{
scanf("%d",&p);
for(int i=;i<p;i++)
{
bool f=;int t=;
for(int j=;j<=p-;j++)
{
t=(t*i)%p;
if(t==) f=;
}
if(t*i%p!=) f=;
res+=f;
}
printf("%d",res);
return ;
}

Problem A

B:

简单分类

#include <bits/stdc++.h>

using namespace std;
int n,A,I;char s[]; int main()
{
scanf("%d%s",&n,s+);
for(int i=;i<=n;i++)
if(s[i]=='A') A++;
else if(s[i]=='I') I++;
if(!I) printf("%d",A);
else if(I==) printf("");
else printf("");
return ;
}

Problem B

C:

明明写起来很简单的一道题目被我弄复杂了

考虑维护序列中原来的值和增量

由于只要求当前最后一个数的增量,因此可以每次仅在增加的最后一个数打上标记

这样在删除数时将当前末尾的增量向前推就能满足要求

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef pair<int,int> P;
const int MAXN=2e5+;
int n,det[MAXN],dat[MAXN],lst;double sum; int main()
{
scanf("%d",&n);
dat[]=;lst=;
while(n--)
{
int op,x,y;
scanf("%d",&op);
if(op==) scanf("%d%d",&x,&y),sum+=x*y,det[x]+=y;
else if(op==) scanf("%d",&x),sum+=x,dat[++lst]=x;
else sum-=dat[lst]+det[lst],det[lst-]+=det[lst],det[lst--]=;
printf("%.6lf\n",sum/lst);
}
return ;
}

Problem C

结果我这个zz用树状数组维护了这个东西,还将一个增量写成了赋值WA了好几次……

#include <bits/stdc++.h>

using namespace std;
#define X first
#define Y second
typedef long long ll;
typedef pair<int,int> P;
const int MAXN=2e5+;
int n;ll bit[MAXN<<],dat[MAXN<<],lst,sum; void Update(int pos,int x)
{while(pos<=n) bit[pos]+=x,pos+=pos&(-pos);}
ll Query(int pos)
{ll ret=;while(pos) ret+=bit[pos],pos-=pos&(-pos);return ret;} int main()
{
scanf("%d",&n);
lst=;sum=;
for(int i=;i<=n;i++)
{
int op,x,y;
scanf("%d",&op);
if(op==)
{
scanf("%d%d",&x,&y);
Update(,y);Update(x+,-y);
dat[]+=y;dat[x+]-=y;
sum+=x*y;
}
else if(op==)
{
scanf("%d",&x);
Update(lst+,x);Update(lst+,-x);
dat[++lst]+=x;dat[lst+]=-x;
sum+=x;
}
else
{
sum-=Query(lst);
Update(lst+,-dat[lst+]);
Update(lst,dat[lst+]);
dat[lst]+=dat[lst+];dat[lst+]=;lst--;
}
printf("%.6lf\n",1.0*sum/lst);
}
return ;
}

Problem C

D:

可以发现真正的状态数只有$2*n$个,对于每一位只有加/减两种状态

发现$y$的值就是走过的总距离,这样用$dp[n][2]$记忆化搜索在当前状态走到边界有多远

每个状态在入栈时打上标记,如果再次走到形成环就是循环了

这样对于每个$i$就是对应于$dp[i+1][0]$时的解

注意将边界$dp[1][0/1]$设为-1

#include <bits/stdc++.h>

using namespace std;
#define X first
#define Y second
typedef long long ll;
typedef pair<int,int> P;
const int MAXN=2e5+;
int n,dat[MAXN],vis[MAXN][];ll dp[MAXN][]; ll dfs(int x,int d)
{
if(x<=||x>n) return ;
if(vis[x][d]==) return -;
if(vis[x][d]==) return dp[x][d];
vis[x][d]=;
ll dist=d?dfs(x+dat[x],):dfs(x-dat[x],);
vis[x][d]=;
if(dist==-) return dp[x][d]=-;
else return dp[x][d]=dist+dat[x];
} int main()
{
scanf("%d",&n);
for(int i=;i<=n;i++) scanf("%d",&dat[i]);
vis[][]=vis[][]=;
for(int i=;i<=n;i++)
{
if(!vis[i][]) dfs(i,);
printf("%I64d\n",dp[i][]==-?-:i+dp[i][]-);
}
return ;
}

Problem D

E:

由于每个点只有一个后继,因此最终图只有环和链

将环的情况排除后就只剩链了

对于一条链,一个可行解为$a_1*w_1+a_2*w_2...+a_n*w_n(a_1>a_2>...>a_n)$

但由于后面的限制条件不能直接背包,考虑将式子变换后消除限制:

设$d_i=a_i-a_{i-1},suf_i=w_i+w_{i+1}...+w_n$,

这样式子变为$d_1*suf_1+d_2*suf_2...+d_n*suf_n$

发现就是一个完全背包,且可以将各个链合并在一起考虑

但注意,这个背包要求$d_1,d_2...d_{n-1}>0$,否则无法保证条件成立

也就是原式每个物品的个数要达到${n-1,n-2...2,1,0}$的下限!

这样就先都取一个,并将其从$t$中减去再跑背包即可

#include <bits/stdc++.h>

using namespace std;
#define X first
#define Y second
typedef long long ll;
typedef pair<int,int> P;
const int MAXN=1e5+,MOD=1e9+;
int n,q,t,x,y,a[MAXN],nxt[MAXN],d[MAXN],dp[MAXN],cnt; int main()
{
scanf("%d%d%d",&n,&q,&t);
for(int i=;i<=n;i++) scanf("%d",&a[i]);
for(int i=;i<=q;i++)
scanf("%d%d",&x,&y),nxt[x]=y,d[y]++;
for(int i=;i<=n&&t>=;i++)
{
if(d[i]) continue;
int sum=;//t<0一定要及时退出!
for(int j=i;j&&t>=;j=nxt[j])
{
cnt++;
sum+=a[j];a[j]=sum;
if(nxt[j]) t-=a[j];
}
}
if(cnt<n||t<) return puts(""),; dp[]=;
for(int i=;i<=n;i++)
for(int j=a[i];j<=t;j++)
(dp[j]+=dp[j-a[i]])%=MOD;
printf("%d",dp[t]);
return ;
}

Problem E

同时还有一个注意点:$t<0$时要及时退出否则会爆$int$

最新文章

  1. [原创]AD9212采样方法
  2. IntelliJ IDEA License Server本地搭建教程
  3. Hbase之遍历超时处理
  4. vs2013 设置为中文版
  5. MongoDB性能优化指南
  6. .NET的微型Web框架 Nancy
  7. 微信小程序之注册和入门
  8. 连表查询都用Left Join吧
  9. Pushgateway 介绍
  10. rocketmq双主模式
  11. SpringSecurityOAuth使用JWT Token实现SSO单点登录
  12. Ex 6_20 最优二叉搜索树..._第六次作业
  13. 代理本地局域网其他PC上网的设置
  14. CSS渐变字体、镂空字体、input框提示信息颜色、给图片加上内阴影、3/4圆
  15. linux常用命令:date 命令
  16. 【Python】【元编程】【二】【描述符】
  17. AE(ArcEngine)定制工具Tool工具箱
  18. Git之代码合并及分支管理
  19. python 统计MySQL大于100万的表
  20. Openstack实现共有云多flat网络

热门文章

  1. JS之document例题讲解1(两张表之间数据转移、日期时间选择、子菜单下拉、用div做下拉菜单、事件总结)
  2. HH实习 acm算法部 1689
  3. Centos 7 安装jdk1.7
  4. LCD实验学习笔记(十):TFT LCD
  5. nmap导出处理脚本
  6. 【EverydaySport】健身笔记——背部训练
  7. 【Python学习笔记】使用Python进行主成分分析
  8. Linux 入门记录:十、Linux 下获取帮助
  9. 004_ssh连接慢的问题的解决?
  10. MyEclipse/Eclipse安装插件的几种方式