HDU 1271 数论，有点巧妙，改天再研究

来源：互联网发布：verizon的网络制式全网编辑：程序博客网时间：2024/06/04 19:05

假设A中去掉的数在第k+1位，可以把A分成三部分，低位，k，和高位。

A == a + b * 10^k + c * 10^(k+1)
B == a + c * 10^k

N == A + B == 2 * a + b * 10^k + c * 10^k * 11

其中b是一位数，b * 10^k不会进位，用10^k除N取整就可以得到b + 11c，再用11除，商和余数就分别是c和b了。但是这里有个问题a是一个小于10^k的数没错，但是2a有可能产生进位，这样就污染了刚才求出来的b + 11c。但是没有关系，因为进位最多为1，也就是b可能实际上是b+1，b本来最大是9，那现在即使是10，也不会影响到除11求得的c。因此c的值是可信的。然后根据2a进位和不进位两种情况，分别考虑b要不要-1，再求a，验算，就可以了。

迭代k从最低位到最高位做一遍，就可以找出所有可能的A。

# include<stdio.h>#include <stdlib.h>int cmp(const void *a,const void *b){return *(int *)a - *(int *)b;}int main(){int n,a,b,c,count,k,s[100],i;while(scanf("%d",&n)!=EOF && n){count=0;for(k=1;k<=n;k*=10){c=(n/k)/11;b=n/k-c*11;if((b!=0 || c!=0) && b<10){a=(n-b*k-c*11*k)/2;if(2*a+b*k+c*11*k==n){count++;s[count]=a+b*k+c*10*k;}}b--;if((b!=0 || c!=0) && b>=0){a=(n-b*k-c*11*k)/2;if(2*a+b*k+c*11*k==n){count++;s[count]=a+b*k+c*10*k;}}}if(count==0) printf("No solution.\n");else{qsort(s+1,count,sizeof(s[1]),cmp);printf("%d",s[1]);for(i=2;i<=count;i++){if(s[i]!=s[i-1])printf(" %d",s[i]);}printf("\n");}}return 0;}

#include<iostream>#include<cstdio>#include<algorithm>#include<set>using namespace std;void solve(int x,set<int> &result){int k,high,b2,a;int c,b1;for(k=1;k<=x;k*=10){high=x/k;c=high/11;c*=k;b1=high%11;if((b1!=0||c!=0)&&b1<10){b1*=k;a=(x-b1-11*c)/2;if(2*a+b1+11*c==x)result.insert(a+b1+c*10);   }b2=high%11-1;if((b2!=0||c!=0)&&b2>=0){b2*=k;int a2=(x-b2-11*c)/2;if(x==2*a2+b2+11*c)result.insert(a2+b2+10*c);}}}int main(){int x;while(cin>>x,x){set<int>result;//set结构不仅可以排序，而且可以去除重复的solve(x,result);if(result.empty()){printf("No solution.\n");}else{set<int>::iterator it=result.begin();printf("%d",*it);while(++it!=result.end())printf(" %d",*it);printf("\n");}}return 0;}