JZOJ2017.08.18 B组

T1

从前在一个美好的校园里，有一只（棵）可爱的弯枝理树。她内敛而羞涩，一副弱气的样子让人一看就想好好疼爱她。仅仅在她身边，就有许多女孩子想和她BH，比如铃，库特，等等。不过，除却巫山不是云，理树的心理只有那个帅气高大的男孩子——恭介，这让女孩子们不得不终日唉声叹气，以泪洗面。不过恭介是那样强大而完美，根本没有办法击败他，她们也只好咬牙忍痛度日，以待反击之时。
终于，她们获得了一次机会。机智的库特利用弹道学、密码学、宇宙学的知识设计出了一个密室，可以让进入的人无法从内部打开出口。库特设计密码的过程很奇葩，是由两个用整数坐标表示的n 维向量导出的。神奇的是，对于这两个向量中的任意一个，无论如何将它的坐标打乱(例如(a1,a2,a3)变成(a3,a1,a2)），打乱后的数量积都不会比原来的两个向量的数量积小。而库特就把原来的两个向量的数量积作为了密码。现在她们只用把恭介引入就可以了。但是，好事多磨，由于她们的粗心大意，在测试密室的时候不小心把自己给关了进去，而且还带走了密码纸。在外面的铃只找到了库特写着两个打乱后的向量的草稿。哇呼~能不能解救这些萌妹子，就看你了。

Input

三行。第一行一个整数N，表示N 维。
第2~3 行每行N 个整数，表示打乱后的两个向量(a1,a2,a3,a4…an),(b1,b2,b3,b4…bn).

Output

如题目要求，输出库特设计的密码

思路：

其实题意是两个序列相乘，乘积的最小值。
贪心，不解析。。。。。。

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
bool cmp1(long a,long b){return a<b;}
bool cmp2(long a,long b){return a>b;}
int main()
{
    long long a[1011],b[1011],n=0;
    scanf("%lld",&n);
    for (int i=1; i<=n; i++) scanf("%lld",&a[i]);
    for (int i=1; i<=n; i++) scanf("%lld",&b[i]);
    sort(a+1,a+n+1,cmp1);
    sort(b+1,b+n+1,cmp2);
    long long ans=0;
    for (int i=1; i<=n; i++)
    {

        ans+=(a[i]*b[i]);

    }
    printf("%lld",ans);
}

T2

Description

终于，在众亲们的奋斗下，最终boss 恭介被关进了库特设计的密室。正当她们松了一口气时，这个世界却发生了天翻覆地的变化：地面开始下沉，天空开始变成血红色，海水沸腾……一幅世界末日的图景。美鱼从她手中的古籍《若山牧水诗歌集》中发现了原因：白鸟は かなしからずや 空の青 海のあをにも 染まずただよふ 。大（xia）意（shuo）就是狡猾的恭介在创造这个世界的时候就篡改了法则。而这个法则的起源，就是一只生死之间的猫。这个猫被关在一个黑盒子里，盒子里有两个毒气罐，如果有任意一个毒气罐被打开那么猫将会被杀死，法则也能得到纠正。然而外界能控制的仅仅是这两个毒气罐被打开的概率。假设第一个毒气罐被打开的概率为1/x，第二个毒气罐为1/y（x,y 为正整数），那么当两个概率和为1/(n!)时，猫将会被莫名其妙地杀死。现在美鱼想知道，有多少对(x,y)可以让猫被莫名其妙杀死。

Input

一行，一个正整数n

Output

一行，满足题意的(x,y)对数。

思路：

设 m=n! ，由等式知x,y必定大于n!，所以再设 x=n!+k=m+k 带入 1/m=1/x+1/y 中化简得到y=m*m/k+m，因y为整数，所以要求k整除m*m，即k为m*m的因子，问题便转化为求n!n!的因子个数， 设n!=p1^e1 p2^e2 * p3^e3 …*pk^ek，则 n!*n!= p1^(2*e1) p2^(2*e2) …*pk^(2*ek) 。 则因子个数sum=(2*e1+1)(2*e2+1)…(2*ek+1)。
对于每个质数p，它的出现次数等于[N/p]+[N/p^2]+…，这样对于每个p是O(log_p(N))的，而1～N之间的质数有O(N/lnN)个，均摊一下就是O(n)的。
至于怎么算每个质数p出现的次数，p的次数等于[n/p]+[n/p^2]+[n/p^3]+….这个可以看成是把1~n里p的倍数的数字都除掉一个p，这样的数有[n/p]个，除完之后这[n/p]个数会变成1~[n/p]，然后再除掉一个p，就是[[n/p]/p]个，然后变成1~[n/p^2]，一直做下去。

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#define mo 100000000
using namespace std;
int len,n;
long long ans[30000],k,x,x1;
bool bz[700010];
void mul(int x)
{
    for(int i=1;i<=len;i++) ans[i]*=x;
    for(int i=1;i<=len;i++)
    {   
        ans[i+1]+=ans[i]/mo;
        ans[i]%=mo;
    }
    while(ans[len+1]>0)
    {
        len++;
        ans[len+1]+=ans[len]/mo;
        ans[len]%=mo;
    }
}
int main()
{
    k=len=1;
    scanf("%d",&n);
    ans[1]=1;
    for(int i=2;i<=n;i++)
    if(!bz[i])
    {   
        x=n,x1=0;
        while(x>=i) x/=i,x1+=x;
        x1=x1*2+1;
        if(x1*k>mo)
        {
            mul(k);
            k=1;
        }
        k*=x1;
        for(int j=i;j<=n;j+=i) bz[j]=1;
    }
    if(k>1) mul(k);
    printf("%lld",ans[len]);
    for(int i=len-1;i>0;i--) printf("%08lld",ans[i]);
}

T3

Description

在美鱼和理树后援团拯救世界的同时，外表柔弱的理树也开始坚强起来，思考着离开这个世界的办法。误打误撞地，她遇上了正在教室破坏课桌打开迷宫入口的沙耶。沙耶告诉理树，这个世界的出口就是这个迷宫的出口。于是理树毫不犹豫地跟沙耶一起跳进了迷宫。在迷宫里，两个女孩子互帮互助，一会儿割绳子，一会儿泡温泉，一会儿雕冰块，跌跌撞撞地走到了终点。不出所料，终点也有一个机关在等着她们。
终点的机关是一个立着的m*n 的方格棋盘，在有些格子上放了一个玩偶，而有些地方直接挖了个大坑。只有取走所有玩偶才能打开出口。但是，由于奇怪的设定，理树和沙耶不能直接触碰玩偶，他们需要操纵机器人来收集它。机器人的走法很奇怪，和国际象棋的马有点像，只不过马可以走任意方向的1*2 路线，它们只会由上往下走r*c(或c*r)的路线，不能回头。而机器人一旦经过一个有玩偶的格子，那个格子上的玩偶将被回收，并且在机器人离开时，那个格子会变成一个坑。理树可以把机器人放在任何一个有玩偶的格子上作为起点，也可以在任何一个有玩偶的格子回收机器人。机器人行走可以视为瞬移，只不过每一次设置新起点都会消耗1 时间。并且，有坑的格子不能落脚。

就在这个紧要关头，玩偶狂热爱好者的沙耶却流着口水智商归0。理树不得不转而求助你，帮忙计算出最少多少时间就能收集到所有玩偶。

Input

第一行包含4 个整数M、N、R、C，意义见问题描述。接下来M 行每行一个长度为N 的
字符串。如果某个字符是’.’，表示这个地方有一个玩偶；如果这个字符是’x’，表示这个地

方是坑。

Output

输出一个整数，表示最短时间。

思路：

其实这题也有用贪心的思想，答案=在每个实心点放的时间-最多可以瞬移的个数
那么我们可以将每一个实心点作为起点做一遍，求出他能到达的第一个没有做过的点
如果当前找到了一个做过的点，那么向上求出他的父亲能否换一个点，如果可以就返回true
以上可以用递归实现

代码：

var m,n,r,c,ans,max,i,j,doll:longint;
    w:array[1..50,1..50]of boolean;
    map:array[1..50,1..50,1..2]of longint;
    f:array[-50..100,-50..100]of boolean;
    dx,dy:array[1..4]of longint;
    a:char;
procedure init;
begin
  fillchar(f,sizeof(f),false);
  readln(n,m,r,c);
  for i:=1 to n do
  begin
    for j:=1 to m do
    begin
      read(a);
      if a='.' then begin inc(doll); f[i,j]:=true; end
      else f[i,j]:=false;
    end;
    readln;
  end;
  if r<>c then
  begin
    dx[1]:=r; dy[1]:=c;
    dx[2]:=r; dy[2]:=-c;
    dx[3]:=c; dy[3]:=r;
    dx[4]:=c; dy[4]:=-r;
  end else
  begin
    dx[1]:=r; dy[1]:=c;
    dx[3]:=r; dy[3]:=-c;
    dx[2]:=50; dy[2]:=50;
    dx[4]:=50; dy[4]:=50;
  end;
  fillchar(map,sizeof(map),0);
end;
function try(x,y:longint):boolean;
var  x1,y1,i:longint;
begin
  for i:=1 to 4 do
  begin
    x1:=x+dx[i];
    y1:=y+dy[i];
    if (f[x1,y1]=true)and(w[x1,y1]) then
    begin
      w[x1,y1]:=false;
      if (map[x1,y1,1]=0)or(try(map[x1,y1,1],map[x1,y1,2])) then
      begin
        map[x1,y1,1]:=x;
        map[x1,y1,2]:=y;
        exit(true);
      end;
    end;
  end;
  exit(false);
end;
procedure main;
begin
  init;
  for i:=1 to n do
    for j:=1 to m do
      if f[i,j]=true then
      begin
        fillchar(w,sizeof(w),true);
        if try(i,j) then inc(ans);
      end;
  write(doll-ans);
end;
begin
  main;
end.