hdu5752 Sqrt Bo
题意:给一个数n,n < 10^100,操作f开根取下整,问几次可以得到1,5次以上直接“TAT”
题解:肯定是大数,但JavaBigInteger没有直接开方的函数,倒过来想,平方5次最大得到多少即可
import java.util.Scanner;
import java.math.BigInteger;
public class Main {
public static void main(String[] args){
Scanner cin = new Scanner(System.in);
BigInteger a1 = BigInteger.valueOf(4);
BigInteger a2 = a1.multiply(a1);
BigInteger a3 = a2.multiply(a2);
BigInteger a4 = a3.multiply(a3);
BigInteger a5 = a4.multiply(a4);
while(cin.hasNext()){
BigInteger a = cin.nextBigInteger();
if(a.compareTo(BigInteger.ZERO) == 0){
System.out.println("TAT");
}
else if(a.compareTo(a1) < 0){
System.out.println("1");
}
else if(a.compareTo(a2) < 0){
System.out.println("2");
}
else if(a.compareTo(a3) < 0){
System.out.println("3");
}
else if(a.compareTo(a4) < 0){
System.out.println("4");
}
else if(a.compareTo(a5) < 0){
System.out.println("5");
}
else{
System.out.println("TAT");
}
}
}
}
hdu5753 Permutation Bo
题意:c数组和h数组,如果h[i]>h[i-1] && h[i]>h[i+1],答案就加上c[i],h[1]-h[n]是1-n的全排列,h[0] = 0,h[n+1] = 0
题解:考虑连续的三位,由于互不相等,排列共有六种,其中有两种是算作贡献的(中间那位),所以贡献为h[i]/3
再考虑两边,贡献为h[i]/2
需要注意的n == 1的情况
int main(){
int n,temp;
while(cin>>n){
double ans=0;
if(n==1){
cin>>temp;
ans+=(double)temp;
printf("%.6f\n",ans);
}
else {
for(int i=0;i<n;i++){
cin>>temp;
if(i==0 || i==n-1){
ans += temp*1.0/2;
}
else {
ans += temp*1.0/3;
}
}
printf("%.6f\n",ans);
}
}
return 0;
}
hdu5754 Life Winner Bo
题意:分各种棋子,给一个初始(n,m),谁能走到(1,1)谁赢
题解:http://www.matrix67.com/blog/archives/6784
直接上男神博客,讲的特别清楚,尤其是皇后
其中没有讲马,不过马也是很好理解的
int main(){
int t,x,a,b;
cin>>t;
while(t--){
cin>>x>>a>>b;
if(x==1){
if(a&1 && b&1) printf("G\n");
else printf("B\n");
}
else if(x==2){
if(a==b) printf("G\n");
else printf("B\n");
}
else if(x==3){
///这些点可以经过偶数步,必定到达(1,1)
if(a==b && a%3==1) printf("G\n");
///这些点一步后,可以到达必胜点
else if(abs(a-b)==1 && (a+b)%3==2) printf("B\n");
///其他点均为不可到达点
else printf("D\n");
}
else {
a--;
b--;
if(a<b)
swap(a,b);
int k=a-b;
a=(int)(k*(1+sqrt(5))/2.0);
if(a==b) printf("G\n");
else printf("B\n");
}
}
return 0;
}
hdu5755 Gambler Bo
题意:给一个30*30 的矩阵,每个元素均为0,1,2,对(x,y)进行一个操作,可以将(x,y)+2,周围四个位置均+1,问经过哪些操作后,可以变成一个全0的矩阵
题解:将每个位置都弄成一个未知数,然后就成了一个900*900的方程,高斯消元的话,是O(n^3),对于这个题来说,每一个横行或竖行,最多就五个未知数,剩余的全为0,因此,可以搞
const int mod = 3;
int exgcd(int a,int b,int &x,int &y){
if(!b){x = 1; y = 0; return a;}
else{
int r = exgcd(b,a%b,y,x);
y -= x * (a/b);
return r;
}
}
int lcm(int a,int b){
int x = 0, y =0;
return a / exgcd(a,b,x,y) * b;
}
int A[MAX][MAX],free_x[MAX],x[MAX];
void Gauss(int n,int m){
int r,c;
for(r=0,c=0;r<n && c<m;c++){
int maxr = r;
for(int i=r+1;i<n;i++) if(abs(A[i][c]) > abs(A[maxr][c])) maxr = i;
if(maxr != r) for(int i=c;i<=m;i++) swap(A[r][i],A[maxr][i]);
if(!A[r][c]) continue;
for(int i=r+1;i<n;i++) if(A[i][c]){
///这里要保证运算都是整数,所以要求最小公倍数
///范围不大可以用int,int运算更快
int d = lcm(A[i][c],A[r][c]);
int t1 = d / A[i][c], t2 = d / A[r][c];
for(int j=c;j<=m;j++)
A[i][j] = ((A[i][j] * t1 - A[r][j] * t2) % mod + mod) % mod;
}
r++;
}
for(int i=r;i<n;i++) if(A[i][m]) return ;
///这里保证是没有自由变元的情况下。
///有自由变元的时候,不一定是x[i] 对应 A[i][m]应该找那一行最开始的一列不为0的那个
for(int i=r-1;i>=0;i--){
x[i] = A[i][m];
for(int j=i+1;j<m;j++){
x[i] = ((x[i] - A[i][j] * x[j]) % mod + mod) % mod;
}
int x1 = 0,y1 = 0;
///这里是用exgcd求逆元,也可以用费马小定理求,如果mod是素数
int d = exgcd(A[i][i],mod,x1,y1);
x1 = ((x1 % mod) + mod) % mod;
x[i] = x[i] * x1 % mod;
}
}
void Gauss_init(){
mem(A,0); mem(free_x,0); mem(x,0);
}
int dx[4]={0,1,0,-1};
int dy[4]={1,0,-1,0};
int a[34][34];
int main(){
//freopen("in.txt","r",stdin);
int t;
cin>>t;
while(t--){
int n,m;
cin>>n>>m;
for(int i=0;i<n;i++){
for(int j=0;j<m;j++)
scanf("%d",&a[i][j]);
}
Gauss_init();
for(int i=0;i<n*m;i++){
int r=i/m;
int c=i%m;
A[i][i]=2;
for(int j=0;j<4;j++){
int x=r+dx[j];
int y=c+dy[j];
if(x>=0 && y>=0 && x<n && y<m){
A[i][x*m+y]=1;
}
}
A[i][n*m]=(3-a[r][c])%3;
}
Gauss(n*m,n*m);
vector<int> v;
for(int i=0;i<n*m;i++){
while(x[i]){
v.push_back(i);
x[i]--;
}
}
printf("%d\n",v.size());
for(int i=0;i<v.size();i++){
int r=v[i]/m;
int c=v[i]%m;
printf("%d %d\n",r+1,c+1);
}
}
return 0;
}
hdu5761 Rower Bo
题意:从(0,a)出发,到(0,0)去,v1船速始终指向(0,0),v2水速始终指向x正方向
题解:当时做的时候,列了微分方程,但是没有解出来,数学太差,哎
int main(){
//freopen("in.txt","r",stdin);
double a,v1,v2;
while(cin>>a>>v1>>v2){
if(a==0)
cout<<"0.000000"<<endl;
else{
if(v1<=v2)
cout<<"Infinity"<<endl;
else{
double ans=(a*v1)/(v1*v1-v2*v2);
printf("%.10lf\n",ans);
}
}
}
}
hdu5762 Teacher Bo
题意:给10^5个点,问有没有相同的曼哈顿距离
题解:当时做的时候,没有太大的思路,现在看来,好水
由于给了距离m,所以不同点之间,最小1,最大2*m,而且本题只问有没有,直接做标记就好
int vis[200010];
int a[100010][3];
int main(){
int t,n,m;
cin>>t;
while(t--){
cin>>n>>m;
for(int i=0;i<n;i++){
scanf("%d%d",&a[i][0],&a[i][1]);
}
bool bool1=false;
memset(vis,0,sizeof(vis));
for(int i=0;i<n;i++){
for(int j=i+1;j<n;j++){
int temp=abs(a[i][0]-a[j][0])+abs(a[i][1]-a[j][1]);
if(vis[temp]){
bool1=true;
break;
}
else {
vis[temp]=1;
}
}
if(bool1){
break;
}
}
if(bool1) printf("YES\n");
else printf("NO\n");
}
}