首页 > 代码库 > 切糕(bzoj 3144)
切糕(bzoj 3144)
Description
Input
第一行是三个正整数P,Q,R,表示切糕的长P、 宽Q、高R。第二行有一个非负整数D,表示光滑性要求。接下来是R个P行Q列的矩阵,第z个 矩阵的第x行第y列是v(x,y,z) (1≤x≤P, 1≤y≤Q, 1≤z≤R)。
100%的数据满足P,Q,R≤40,0≤D≤R,且给出的所有的不和谐值不超过1000。
Output
仅包含一个整数,表示在合法基础上最小的总不和谐值。
Sample Input
2 2 2
1
6 1
6 1
2 6
2 6
1
6 1
6 1
2 6
2 6
Sample Output
6
HINT
最佳切面的f为f(1,1)=f(2,1)=2,f(1,2)=f(2,2)=1
/* 个人感觉很神的一道题目。 如果有解的话,会有一个p满足:(10^x-1)/9*8=L*p => 10^x-1=9*L*p/8 设m=9*L/gcd(L,8) 则存在p1使得 10^x-1=m*p1 => 10^x=1(mod m) 根据欧拉定理 10^φ(m)=1(mod m) 所以x一定是φ(m)的因数(这好像是某个定理)。 */ #include<iostream> #include<cstdio> #define lon long long #define N 400010 #ifdef unix #define LL "%lld" #else #define LL "%I64d" #endif using namespace std; int prime[N],f[N],num,qlen; lon q[N]; void get_prime(){ for(int i=2;i<N;i++){ if(!f[i]) prime[++num]=i; for(int j=1;j<=num;j++){ if(i*prime[j]>=N) break; f[i*prime[j]]=1; if(i%prime[j]==0) break; } } } lon gcd(lon a,lon b){ if(!b) return a; return gcd(b,a%b); } lon euler(lon x){ lon res=x; for(lon i=2;i*i<=x;i++) if(x%i==0){ res-=res/i; while(x%i==0) x/=i; } if(x>1) res-=res/x; return res; } void get_q(lon n){ qlen=0; for(int i=1;i<=num&&n>1;i++){ while(n%(lon)prime[i]==0){ n/=prime[i]; q[++qlen]=prime[i]; } } if(n>1) q[++qlen]=n; } lon mul(lon a,lon b,lon mod){ lon base=a,r=0; while(b){ if(b&1) r+=base;r%=mod; base+=base;base%=mod; b>>=1; } return r; } lon poww(lon a,lon b,lon mod){ lon base=a,r=1; while(b){ if(b&1) r=mul(r,base,mod); base=mul(base,base,mod); b>>=1; } return r; } int main(){ int cas=0;lon L;get_prime(); while(scanf(LL,&L)&&L){ lon m=9*L/gcd(L,8); if(gcd(m,10)>1){ printf("Case %d: 0\n",++cas); continue; } lon x=euler(m);get_q(x); for(int i=1;i<=qlen;i++){ if(poww(10,x/q[i],m)==1){ x/=q[i]; } } printf("Case %d: ",++cas); printf(LL,x);printf("\n"); } return 0; }
切糕(bzoj 3144)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。