沉迷T3,不能自拔。
T1:
题意:求2的n次方modp,p为质数;
分析:
如果n小,就是快速幂。
如果n大,又是质数,就是费马小定理先把nmod了再直接暴力算。
有人写高精度快速幂,过了70,辛苦了。
#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<ctime>
#define ll long long
#define inf 2e8
#define modd 1e9+7
#define clr(x) memset(x,0,sizeof(x))
#define maxen(x) memset(x,127,sizeof(x))
#define maxer(x) memset(x,31,sizeof(x))
#define minus(x) memset(x,-1,sizeof(x))
#define each(i,n) for(int i=1;i<=n;i++)
#define minn(a,b,c) min(a,min(b,c))
#define maxx(a,b,c) max(a,max(b,c))
#ifdef WIN32
#define lld "%I64d"
#else
#define lld "%lld"
#endif
#define PROC "cell"
//for(int i=1;i<=n;i++)
//(double) (ll) LL (int)
//(double)clock()/CLOCKS_PER_SEC
using namespace std;
const int Maxn=3e5+5;
int len,num,p,ans;
char a[Maxn];
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void init()
{
ans=1;
scanf("%s",a);
len=strlen(a);
scanf("%d",&p);
if(len>4){
for(int i=0;i<len;i++)
num=(num*10+a[i]-'0')%(p-1);}
else{
for(int i=0;i<len;i++)
num=num*10+a[i]-'0';
}
}
void work()
{
for(int i=1;i<=num;i++)
ans=ans*2%p;
printf("%d",ans);
}
void debug()
{
//
}
int main()
{
freopen(PROC".in","r",stdin);
freopen(PROC".out","w",stdout);
init();
work();
//debug();
return 0;
}
T2:
题意:
在立方体里装管道,给定出发点和到达点,求最少使用管道数。
分析:
暴力,三维加上一个管道形态,加上一个流出方向。
推增量数组需要耐心。
#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<ctime>
#define ll long long
#define inf 2e8
#define modd 1e9+7
#define clr(x) memset(x,0,sizeof(x))
#define maxen(x) memset(x,127,sizeof(x))
#define maxer(x) memset(x,31,sizeof(x))
#define minus(x) memset(x,-1,sizeof(x))
#define each(i,n,m) for(int i=n;i<m;i++)
#define minn(a,b,c) min(a,min(b,c))
#define maxx(a,b,c) max(a,max(b,c))
#ifdef WIN32
#define lld "%I64d"
#else
#define lld "%lld"
#endif
#define PROC "tube"
//for(int i=1;i<=n;i++)
//(double) (ll) LL (int)
//(double)clock()/CLOCKS_PER_SEC
using namespace std;
const int Maxn=25;
int zlx[10][10]={{2,2,3,3,2,2,3,3},
{0,2,0,1,0,-2,0,-1},
{2,0,1,0,-2,0,-1,0},
{-2,-2,-3,-3,-2,-2,-3,-3},
{0,2,0,1,0,-2,0,-1},
{2,0,1,0,-2,0,-1,0}};
int zly[10][10]={{0,2,0,1,0,-2,0,-1},
{2,2,3,3,2,2,3,3},
{0,2,0,1,0,-2,0,-1},
{0,2,0,1,0,-2,0,-1},
{-2,-2,-3,-3,-2,-2,-3,-3},
{0,2,0,1,0,-2,0,-1},
};
int zlz[10][10]={{2,0,1,0,-2,0,-1,0},
{2,0,1,0,-2,0,-1,0},
{2,2,3,3,2,2,3,3},
{2,0,1,0,-2,0,-1,0},
{2,0,1,0,-2,0,-1,0},
{-2,-2,-3,-3,-2,-2,-3,-3}};
int zls[10][10]={{3,2,3,2,6,5,6,5},
{3,1,3,1,6,4,6,4},
{1,2,1,2,4,5,4,5},
{3,2,3,2,6,5,6,5},
{3,1,3,1,6,4,6,4},
{1,2,1,2,4,5,4,5}} ;
const int lim=Maxn*Maxn*Maxn*Maxn;
int stackx[lim],stacky[lim],stackz[lim],stacks[lim];
int xmax,ymax,zmax,stax,stay,staz,stas,finx,finy,finz,fins;
int vis[Maxn][Maxn][Maxn][Maxn],d[Maxn][Maxn][Maxn][Maxn];
int x,y,z,head,tail;
int qx[lim],qy[lim],qz[lim],qs[lim],xx,yy,zz,ss;
char sym,ano;
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void calc(int& s1)
{
sym=getchar();
ano=getchar();
if(ano=='x')s1=1;
if(ano=='y')s1=2;
else if(ano=='z')s1=3;
if(sym=='-')s1+=3;
}
void init()
{
xmax=read();ymax=read();zmax=read();
stax=read();stay=read();staz=read();
calc(stas);
if(stas==1)stax--;
else if(stas==2)stay--;
else if(stas==3)staz--;
else if(stas==4)stax++;
else if(stas==5)stay++;
else if(stas==6)staz++;
finx=read();finy=read();finz=read();
calc(fins);maxer(d);
}
void work()
{
vis[stax][stay][staz][stas]=1;
d[stax][stay][staz][stas]=0;
stackx[++tail]=stax;
stacky[tail]=stay;
stackz[tail]=staz;
stacks[tail]=stas;
head=1;
while(head<=tail){
int x=stackx[head];
int y=stacky[head];
int z=stackz[head];
int s=stacks[head++];
vis[x][y][z][s]=0;
each(i,0,8){
int tox=x+zlx[s-1][i];
int toy=y+zly[s-1][i];
int toz=z+zlz[s-1][i];
int tos=zls[s-1][i];
if(tox<=xmax&&tox>0&&toy<=ymax&&toy>0&&toz<=zmax&&toz>0&&
d[tox][toy][toz][tos]>d[x][y][z][s]+1){
d[tox][toy][toz][tos]=d[x][y][z][s]+1;
if(!vis[tox][toy][toz][tos]){
vis[tox][toy][toz][tos]=1;
stackx[++tail]=tox;
stacky[tail]=toy;
stackz[tail]=toz;
stacks[tail]=tos;
}
}
}
}
if(d[finx][finy][finz][fins]<=6)
printf("%d",d[finx][finy][finz][fins]);
else printf("Impossible");
}
void debug()
{
//
}
int main()
{
freopen(PROC".in","r",stdin);
freopen(PROC".out","w",stdout);
init();
work();
//debug();
return 0;
}
T3:
题意:
在路上走,转弯前后和出发到达都要时间加倍,求最短时间。
分析:
我当成记忆化搜索来做的,刘汝佳拆点,其实表示方法和我一模一样。就是存下来的方向和是否加倍。
不过有了拆点这个概念以后,我们就可以直接SPFA了,比我搜索要快很多还不容易错(在无向图上用记忆化搜索走单源最短真的很艰难啊!)
易错点就是到了终点可能还要走,但更新答案是不再走的(看程序意会)
#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<ctime>
#define ll long long
#define inf 2e8
#define modd 1e9+7
#define clr(x) memset(x,0,sizeof(x))
#define maxen(x) memset(x,127,sizeof(x))
#define maxer(x) memset(x,31,sizeof(x))
#define minus(x) memset(x,-1,sizeof(x))
#define each(i,n) for(int i=1;i<=n;i++)
#define minn(a,b,c) min(a,min(b,c))
#define maxx(a,b,c) max(a,max(b,c))
#ifdef WIN32
#define lld "%I64d"
#else
#define lld "%lld"
#endif
#define PROC "roller"
//for(int i=1;i<=n;i++)
//(double) (ll) LL (int)
//(double)clock()/CLOCKS_PER_SEC
using namespace std;
const int Maxn=105;
const int zl[4][2]={{0,1},{0,-1},{1,0},{-1,0}};
int tov[Maxn][Maxn][4],f[Maxn][Maxn][5][2];
int n,m,stan,stam,finn,finm,ans=2e8;
queue<int>qx,qy,qd,qp,qf;
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void init()
{
n=read();m=read();
stan=read();stam=read();
finn=read();finm=read();
each(i,2*n-1){
if(i%2){
each(j,m-1)
tov[i/2+1][j][0]=tov[i/2+1][j+1][1]=read();
}
else{
each(j,m)
tov[i/2][j][2]=tov[i/2+1][j][3]=read();
}
}
maxer(f);
}
void work()
{
qx.push(stan);
qy.push(stam);
qd.push(4);
qp.push(1);
qf.push(0);
f[stan][stam][4][1]=0;
while(!qx.empty()){
int x=qx.front();qx.pop();
int y=qy.front();qy.pop();
int d=qd.front();qd.pop();
int p=qp.front();qp.pop();
int fr=qf.front();qf.pop();
for(int i=0;i<4;i++)
if(tov[x][y][i]){
int tox=x+zl[i][0];
int toy=y+zl[i][1];
int pdd;
if(i==d)pdd=0;
else pdd=1;
if(f[tox][toy][i][pdd]>
f[x][y][d][p]+(pdd+1)*tov[x][y][i]+((pdd==1&&p==0)?fr:0)){
f[tox][toy][i][pdd]=
f[x][y][d][p]+(pdd+1)*tov[x][y][i]+((pdd==1&&p==0)?fr:0);
qx.push(tox);
qy.push(toy);
qd.push(i);
qp.push(pdd);
qf.push(tov[x][y][i]);
}
if(tox==finn&&toy==finm){
f[tox][toy][i][pdd]=min(f[tox][toy][i][pdd],
f[x][y][d][p]+2*tov[x][y][i]+((p==0&&i!=d)?fr:0));
ans=min(ans,f[x][y][d][p]+2*tov[x][y][i]+
((p==0&&i!=d)?fr:0));
}
}
}
if(ans==2e8)
printf("Impossible");
else printf("%d",ans);
}
void debug()
{
//
}
int main()
{
freopen(PROC".in","r",stdin);
freopen(PROC".out","w",stdout);
init();
work();
//debug();
return 0;
}