T2 GMOJ2127. 电子表格
(File IO): input:excel.in output:excel.out
题目描述
也许你用过Microsoft Excel之类的电子制表软件，这类软件最令人称道的就是强大的公式计算功能。现在希望你也来实现一个具有最基本功能的电子制表软件。表格共有m列(0 < m • 26)，从左到右依次用A到Z的大写英文字母表示；有n行(0 < n <100)，从上到下依次用1到100的整数表示。这样，每一个单元格的位置就可以唯一地用它所在的列和行表示出来，例如从左到右第3列，从上到下第5行的单元格就可以用“ C5”来表示（注意，这里字母和数字中间没有空格）。
现在对表格进行了一系列的操作，这些操作主要就是赋值和查询。定义操作的输入规则
如下：
- 每个操作占一行，根据操作类型的不同，每行中可能有二至四个用空格隔开的“单词”；
- 每行的第一个单词指定了该操作涉及的单元格的位置；
- 每行的第二个单词指定了相应的操作，可能是： input,output,sum,avg
- (1).如果第二个单词是input，表示接下来的一个整数是要赋予该单元格的值，这个值是不超过1000的正整数
- (2).如果第二个单词是output，表示你需要在输出文件中输出这个单元格当前的值
- (3).如果第二个单词是sum，表示接下来输入的两个单词定义了一个矩形区域，该单元格的值就应该恒为这个矩形区域中所包含的单元格的值的和，直到该单元格被重新定义
- (4).如果第二个单词是avg，表示接下来输入的两个单词定义了一个矩形区域，该单元格的值就应该恒为这个矩形区域中所包含的单元格的值的算术平均数，直到该单元格被重新定义；
- “输入的两个单词定义了一个矩形区域”是指输入一个矩形区域的左上角和右下角的单元格的位置，这样就唯一确定了这个矩形区域；
- 所有时刻，每个单元格的值均为整数，如果不是，则向下取整；
- 如果某个单元格的值没有在上文定义，则它的值默认为0；
- 不会出现循环定义的情况；
- 在操作过程中所有单元格的值不超过231-1。
输入
第一行输入两个用空格隔开的正整数m和n，分别代表表格的列数和行数。
第二行输入一个正整数s，表示操作的总数。
以下s行每行输入一个操作，具体格式参见问题描述。
输出
对于输入数据的每一个“ output”操作输出一行结果。因此，输出文件的行数等于输入文
件中“ output”操作的个数。
样例输入
3 5
5
A1 input 100
B2 input 200
C3 sum A1 C2
C5 avg B2 C4
C5 output
样例输入
83
数据范围限制
对于30%的数据， m; n; s <= 10；
对于100%的数据， m <= 26， n < =100， s <=100。
Solution
当我看到上面的数据范围限制时，心中无比的开心……
n2暴力，是个地球人都不会TLE吧？
嗯，没错，我Wonderful Answer了了了了了……
本地过，样例太简单了吧……
每次有新的定义时直接计算该单元格的值是错误的，因为sum 和avg 的单元格必须随时更新（题目说明了这样的单元格必须恒等于要求的数值和或平均数）。
神奇的wrong！我也知道呀！
but我用了个栈来存储指令，每次操作后都让它执行一遍的……
大概是要把输入变代码吧。
没错，本地是过了，但是0分呐！
Code
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<cmath> 5 #include<map> 6 #include<set> 7 #include<vector> 8 #include<algorithm> 9 #define IL inline 10 using namespace std; 11 struct node{ 12 int x,y;//输出单元格 13 bool avg;//0->求和 1->求平均数 14 int x1,x2,y1,y2;//计算对象 15 }order[100]; 16 int tail=0; 17 string str; 18 int m,n,s; 19 int array[27][101]; 20 int main() 21 { 22 // freopen("excel.in","r",stdin); 23 // freopen("excel.out","w",stdout); 24 cin>>m>>n>>s; 25 int x,y; 26 while(s--) 27 { 28 x=getchar(); 29 while(x<'A'||x>'Z'){ 30 x=getchar(); 31 } 32 x-='A'; 33 y=getchar()-'0'; 34 cin>>str; 35 int i=0; 36 bool flag=0; 37 do{ 38 switch(str[i]){ 39 case 'i':{ 40 int value; 41 cin>>value; 42 array[x][y]=value; 43 break; 44 } 45 case 'o':{ 46 cout<<array[x][y]<<endl; 47 break; 48 } 49 case 's':{ 50 order[tail].x=x; 51 order[tail].y=y; 52 53 order[tail].avg=0; 54 55 order[tail].x1=getchar(); 56 while(order[tail].x1<'A'||order[tail].x1>'Z') order[tail].x1=getchar(); 57 order[tail].x1-='A'; 58 order[tail].y1=getchar()-'0'; 59 60 order[tail].x2=getchar(); 61 while(order[tail].x2<'A'||order[tail].x2>'Z') order[tail].x2=getchar(); 62 order[tail].x2-='A'; 63 order[tail].y2=getchar()-'0'; 64 65 tail++; 66 break; 67 } 68 case 'a':{ 69 order[tail].x=x; 70 order[tail].y=y; 71 72 order[tail].avg=1; 73 74 order[tail].x1=getchar(); 75 while(order[tail].x1<'A'||order[tail].x1>'Z') order[tail].x1=getchar(); 76 order[tail].x1-='A'; 77 order[tail].y1=getchar()-'0'; 78 79 order[tail].x2=getchar(); 80 while(order[tail].x2<'A'||order[tail].x2>'Z') order[tail].x2=getchar(); 81 order[tail].x2-='A'; 82 order[tail].y2=getchar()-'0'; 83 84 tail++; 85 break; 86 } 87 default:{ 88 flag=1; 89 break; 90 } 91 } 92 }while(flag&&i++<=3); 93 94 for(int i=0;i<tail;i++) 95 { 96 int tsa=0; 97 for(int tx=order[i].x1;tx<=order[i].x2;tx++) 98 for(int ty=order[i].y1;ty<=order[i].y2;ty++) 99 tsa+=array[tx][ty]; 100 if(order[i].avg) tsa/=(order[i].x2-order[i].x1+1)*(order[i].y2-order[i].y1+1); 101 array[order[i].x][order[i].y]=tsa; 102 } 103 } 104 return 0; 105 }
T4 GMOJ1574. X-因子链
(File IO): input:factor.in output:factor.out
时间限制: 1000 ms 空间限制: 131072 KB 具体限制
题目描述
给一个正整数X，一个长度为m的X-因子链是指这样一个序列：X0=1，X1，X2，。。。，Xm=X满足：Xi<Xi+1同时Xi|Xi+1(Xi+1能被Xi整除)
要求X-因子链的最大长度Len和长度为Len的X-因子链的数量。
输出
一行，两个整数，分别表示最大长度和该长度链的种数。
样例输入
100
样例输出
4 6
数据范围限制
（空）
Solution
稍作思考即可
P1问题转化
有一个数列，已知首项（为1）与尾项，每一项都是前一项的整数倍，且每一项都大于前一项。
如何使此数列最长？
设相邻两项的商为k
那么k一定是越小越好
但是所有k之积必须等于尾项
那么，k即为尾项的所有质因数。
所以第一问只需求x的质因数个数。
P1Code
int maxlen=0x3f3f3f3f,kind,x; int lian[10000000]; IL int prime_factor(int num) { int ans=0; if(num==1) return 1; if(num==2) return 1; int i=2; while(num!=1) { while(num%i==0) { num/=i; ans++; lian[i]++; } i++; } return ans; }
（使用IL（inline）加速调用）
填写lian[i]++是为了第二问的
P2看穿此题
排列组合类型题
因子链的数量，即这些prime factors可以排成多少种
以样例中的100为例
100=2*2*5*5=22*52
那么这串链有两种、四个元素
不去重的话，一共有4!种排列
但是这道题是要问的是种类（即去重的）数量……
现在只观察其中的两个2
2*2*？*？
这时两个2如果交换，会被算入全排列的数量中，但是不能被算入种类数量中！
两个2可以有两种摆放位置的顺序
那么4!=24除去两次2（两个2和两个5的）即可。
若有三个2呢？
三个2可以有3!=3*2*1=6种摆放位置的顺序
把n!除以6即可。
P2数学推导
设
其中ans为质因数的数量
所以
P2Code
若数据较小（不存在的），可使用dfs暴力枚举
1 IL void dfs(int depth) 2 { 3 if(depth==maxlen){ 4 kind++; 5 return; 6 } 7 for(int i=1;i<=x;i++) 8 { 9 if(lian[i]>0){ 10 lian[i]--; 11 dfs(depth+1); 12 lian[i]++; 13 } 14 } 15 }
于是20分……枉费了我辛辛苦苦思考呀！
阶乘运算
1 IL long long factor(int num) 2 { 3 if(num==0) return 1; 4 long long ans=1; 5 for(int i=1;i<=num;i++) 6 ans*=i; 7 return ans; 8 }
如果只是按照普通的阶乘的话，比较容易爆long long
解决方案有两种：
1、预处理（或者打表、记忆化……）出1~20的阶乘（极限应该是31！，但会超过unsigned long long，且数据也没有那么大哈哈），factor函数可用递归。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<cmath> 5 #include<map> 6 #include<set> 7 #include<vector> 8 #include<algorithm> 9 #define IL inline 10 using namespace std; 11 int lian[50000]; 12 unsigned long long kind,maxlen=0x3f3f,x; 13 unsigned long long fact[32]={1,1,2,6,24,120,720,5040,40320,362880,3628800,39916800,479001600,6227020800,87178291200,1307674368000,20922789888000,355687428096000,6402373705728000,121645100408832000,2432902008176640000}; 14 IL unsigned long long factor(int num) 15 { 16 if(fact[num]!=0) return fact[num]; 17 return fact[num]=factor(num-1)*num; 18 }
如果记忆化的话，其实就可以不用打表了（反正都是2、3毫秒的样子）。
2、把一个数的阶乘分解成
这样子
用结构体数组储存。
最后要相除则换成相减即可
——来自某某不让我去B组的老师的想法
Code
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<cmath> 5 #include<map> 6 #include<set> 7 #include<vector> 8 #include<algorithm> 9 #define IL inline 10 using namespace std; 11 int lian[50000]; 12 unsigned long long kind,maxlen=0x3f3f,x; 13 unsigned long long fact[32]={1,1,2,6,24,120,720,5040,40320,362880,3628800,39916800,479001600,6227020800,87178291200,1307674368000,20922789888000,355687428096000,6402373705728000,121645100408832000,2432902008176640000}; 14 IL int prime_factor(int num) 15 { 16 int ans=0; 17 if(num==1) return 1; 18 if(num==2) return 1; 19 int i=2; 20 while(num!=1) 21 { 22 while(num%i==0) 23 { 24 num/=i; 25 ans++; 26 lian[i]++; 27 } 28 i++; 29 } 30 return ans; 31 } 32 IL unsigned long long factor(int num) 33 { 34 if(fact[num]!=0) return fact[num]; 35 return fact[num]=factor(num-1)*num; 36 } 37 int main() 38 { 39 // freopen("factor.in","r",stdin); 40 // freopen("factor.out","w",stdout); 41 scanf("%lld",&x); 42 maxlen=prime_factor(x); 43 printf("%lld ",maxlen); 44 kind=factor(maxlen); 45 for(unsigned long long i=1;i<=sqrt(x);i++) 46 if(lian[i]>1) 47 kind/=factor(lian[i]); 48 printf("%lld",kind); 49 return 0; 50 }
这道题的代码在尝试着极限……0ms，512KB！
可是……真的有0ms……
